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1. 


FOREWORD. 


The  objective  of  this  study  is  to  transplant  some 
new  understanding  of  friction  phenomena  into  the  science  of  explosives 
Initiation  and  safety,  and  to  suggest  fruitful  new  research  programs 
which  will  go  beyond  the  empirical  and  enhance  basic  understanding. 

It  is  essentially  a  planning  study. 

Mechanical  engineers  have  been  studying  friction  since  the 
fourteenth  century,  and  they  have  gained  a  degree  of  competence  in 
dealing  with  it  by  lubrication  and  by  modification  of  surfaces.  Their 
actual  understanding  of  friction  is  still  fragmentary;  but  some  of  the 
mechanical  and  lubrication  art  can  be  applied  to  the  problems  of 
explosives,  and  this  report  gives  an  introduction  to  that  art. 

The  art  selected  for  presentation  is  the  generation  of  heat 
and  hot-spots  by  mechanical  friction.  In  recent  years,  a  number  of 
investigators  have  addressed  themselves  to  quantitative,  mathematical 
models  of  this  phenomenon  in  metals;  and  it  appears  feasible  to  adapt 
their  treatments  to  explosives.  Successful  adaptation  would  lead  to 
a  quantitative  understanding  of  frictional  initiation  and  to  a 
concomitant  delimitation  of  what  can  and  cannot  be  done  about  it.  This 
report  includes  a  suggested  Research  Plan  to  do  that. 

This  study  has  searched  the  general  literature  on  lubrication 
and  micromechanics  through  the  year  1970,  and  has  Included  personal 
consultation  with  selected  authorities  in  the  fields  of  friction  and 
lubrication.  So  much  material  was  found  that  it  was  necessary  continually 
to  narrow  the  scope  of  this  review  in  order  to  preserve  adequate  depth 
in  the  selected  topic;  the  literature  is  enormous,  and  one  could  literally 
spend  a  career  studying  it  and  adapting  it  to  explosives.  Future  workers 
wishing  to  range  wider  than  the  scope  of  this  report  will  find  an  excellent 
starting  point  in  the  general  references  cited. 

This  report  is  the  product  of  Plcatlnny  Arsenal  Contract 
DAAA21-69-C-0558,  carried  out  by  the  author  on  a  non-profit,  spare¬ 
time  basis.  The  courtesy  of  Esso  Research  and  Engineering  Company 
in  approving  this  outside  activity  by  a  full-time  employee  and  in 
providing  library  access  is  gratefully  acknowledged,  as  is  the 
invaluable  counsel  of  the  author's  colleague,  Mr.  Alan  Beerbower. 

The  contract  was  sponsored  by  the  Explosives  Laboratory  and 
monitored  by  Dr.  Joseph  Hershkowits  end  Dr.  Bernard  Pollock 
of  the  Applied  Physics  Branch  thereof. 

JOHN  A.  BROWN 
Principal  Investigator 


2.  SUMMARY 


This  study  introduces  the  knowledge  of  asperity  and 
flash  temperature  statistics  which  has  been  contributed  by  the 
engineering  and  lubrication  fields  and  gives  an  anthology  of  primary 
and  summary  Jources  to  support  a  sound  theoretical  investigation 
of  its  appiiaation  to  the  explosives  field.  It  outlines  a  Research 
Plan  to  derive  an  integrated  mathematical  model  of  frictional  ignition 
and  to  test  its  validity  with  empirical  data.  It  suggests  how  a 
valia..ted  model  will  lead  to  better  safety  data  and  to  better  and 
safer  explosives  handling  techniques. 

Solids  do  not  contact  each  other  over  their  entire  touching 
surfaces,  as  they  appear  to  do,  but  rather  only  on  the  tips  of  the 
higher  surface  asperities.  The  total  of  real  contact  is  such  a 

siall  percentage  of  the  total  surface  area  that  the  touching  asperities 
are  loaded  to  their  yield  point  and  deform  until  their  aggregate  cross 
section  is  just  adequate  to  support  the  load. 

Most  of  the  surface  is  not  even  touching;  but  the  tiny  areas 
which  do  touch  -  the  functions  -  are  in  such  intimate  contact  that 
van  der  Waals,  electrostatic  or  even  Interatomic  bonds  form  and  lead  to 
strong  adhesions,  ^en  one  solid  body  is  forced  to  slide  over  another, 
the  adhesions  are  forcibly  sheared;  and  the  work  required  to  shear  them 
is  the  major  component  of  the  frictional  resistance. 

The  heat  generated  between  two  sliding  bodies,  which  is  equal 
to  the  work  expended  In  making  them  slide,  necessarily  appears  in  and 
only  in  the  junctions;  and  since  the  total  area  of  the  junctions  is 
quite  small,  the  resulting  temperatures  can  be  quite  high  -  hundreds 
or  even  thous'^nds  of  degrees.  If  one  knew  the  size  of  each  junction 
and  the  fraction  of  the  total  frictional  power  it  consumed,  one  could 
readily  calculate  the  resulting  temperature  in  that  particular  hot-spot. 

The  sl»o  of  individual  junctions  is  usually  not  known;  but 
the  average  size  can  be  estimated,  and  the  distribution  of  sizes 
appears  to  be  Gaussian  or  nearly  so  in  a  broad  range  of  cases,  so  that 
at  least  a  atatietiaal  deeoription  of  the  flash  teripcraturas  on  the 
asperity  tips  can  be  derived.  This  function,  coupled  with  published 
models  of  the  growth  of  hot-spots  to  explosion,  offers  a  way  to 
calculate  the  probability  of  ignition  of  a  given  explosive  in  a  given 
frictional  situation,  entirely  from  first  principles.  Comparison  with 
empirical  data  will  then  serve  to  test  the  validity  of  the  otodel. 

A  number  of  alternate  models  can  be  postulated,  involving 
different  assumptions  as  to  physical  mechanisms  of  heat  generecion 
and  as  to  what  surfaces  art  and  are  not  critical;  and  there  is  a 
laboratory  friction  sensitivity  teeter  cither  existing  or  suggested 
for  each.  Calculation  of  ignition  probabilities  as  above  offers  a 
new  and  powerful  way  to  taat  the  validity  both  of  the  assumptions  and 
of  the  testers,  and  thus  to  gain  a  deeper  understanding  of  toe  entire 
frictional  initiation  process. 


INTENTIONALLY 


BLANK 
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3.  THE  NATURE  OF  MECHANICAL  FRICTION 


The  classical  concepts  of  friction  coupled  with  modern 
concepts  of  the  nature  of  solid-solid  contact  open  up  a  new  approach 
to  understanding  the  basic  nature  of  the  frictional  initiation  of 
explosives. 


THE  CLASSICAL  UNDERSTANDING  OF  FRICTION 


Friction  is  the  resistance  to  motion  which  exists  in  varying 
degree  whnever  one  solid  object  is  caused  to  slide  over  the  surface  of 
another  one.  It  is  a  universal  attribute  of  matter,  and  its  application 
has  its  roots  in  prehistory.  The  use  of  sleds,  rollers  and  wheels,  often 
supplemented  by  liquid  lubricants,  dates  back  more  than  3000  years; 
although  the  scientific  study  of  frictional  phenomena  is  much  more  recent. 

Leonardo  Da  Vinci  was  probably  the  first  purposeful  observer 
of  friction  phenomena.  He  noted  that  the  friction  of  two  sliding  bodies 
is  independent  of  their  contacting  areas  and  that  the  total  frictional 
force  is  proportional  to  the  normal  load  between  the  contacting  surfaces 
(113).  These  two  principles,  still  valid  today,  are  usually  called 
"Amontons*  Laws",  after  the  French  scientist  Amontons  who  rediscovered 
them  In  1699  (114).  (There  are  limits  to  their  validity,  but  mainly  they 
are  true.) 


Amontons  (114),  Coulomb  (115)  and  Euler  (116)  further  quantified 
the  studies  of  friction;  and  by  1785  most  of  today's  general  understanding 
was  established.  These  scientists  hypothesized  that  friction  is  due  to 
the  interlocking  of  mechanical  protuberances  or  asperities  on  the  surfaces 
of  the  contacting  material  like  two  pieces  of  sandpaper  face  to  face. 
This  "Roughness  Hypothesis"  explained  Amontons'  laws  and  remained  the 
majority  view  right  through  the  nineteenth  century  and  into  the  twentieth. 

Beginning  about  1920,  however.  Interest  began  to  revive  in  the 
"Adhesion  Hypothesis"  as  growing  precision  In  measurement  revealed  more 
and  more  weaknesses  In  the  Roughness  Hypothesis.  "Adhesion"  hypothesizes 
that  actual  'Srelding"  occure  wherever  two  surfaces  touch  and  that  the 
frictional  resistance  to  sliding  is  the  force  required  to  break  the  welds. 
Interestingly  enough.  Adhesion  had  been  considered  by  the  original  workers 
and  rejected  because  they  could  not  reconcile  it  with  the  observed  fact 
that  the  frictional  force  was  independent  of  the  contacting  area. 

Today,  it  is  almost  universally  accepted  that  frictional  force 
is  due  to  a  combination  of  the  two  effects:  "Adhesion  plus  Plastic 
Deformation",  wherein  only  a  few  of  the  higher  asperities  of  each  surface 
actually  contact  a  few  asperities  of  the  other  surface.  As  the  nonsal 


load  Is  Increased,  the  high  asperities  are  squashed  down  so  that  the 
next  smaller  asperities  begin  to  contact  each  other;  and  this  process 
continues  until  there  is  enough  contact  area  to  support  the  load.  This 
actual  avea  of  contact  is  very  much  smaller  than  the  apparent  area  of 
contact  (the  gross  size  of  the  object's  face);  In  fact.  It  has  nothing 
to  do  with  the  size  of  the  face,  which  neatly  explains  why  the  frictional 
force  has  nothing  to  do  with  the  size  of  the  face  but  is  proportional  to 
the  normal  load. 


As  the  asperities  contact  one  another,  they  adhere;  and  the 
nature  of  the  adhesion  varies  widely.  Objects  of  like  composition  can 
often  form  actual  welds  or  interatomic  crystal  bonds,  but  there  are 
smaller  adhesive  forces  which  act  between  any  two  materials.  Bowden  ^  2  ) 
describes  them  as  Coulombic  and  Van  der  Waals.  Closer  than  about  20  A,  ^ 
the  force  is  mainly  coulombic  and  is  relatively  strong.  Beyond  about  20  A, 
the  force  is  mainly  Van  der  Waals  and  is  weaker  but  longer  range. 


It  can  thcrafora  ba  aaan  chat  forced  sliding  requires  the 
rupture  of  aaperity  adhesions  and  that  there  will  indeed  be  some  asperity 
interlocking  Just  as  Awontona  had  thought.  Quantitatively  it  tuma  out 
that  most  of  the  obaerved  frictional  force  is  due  to  adhesion,  and  only 
a  small  amount  of  it  is  due  to  asperity  interlocking  or  deformation. 
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A  complication  resides  In  the  fact  that  adhesion  Is  usually 
greatly  Inhibited  by  the  presence  of  an  Interfering  film  of  something 
such  as  a  surface  oxide,  or  moisture,  or  a  foreign  material,  or  even 
adsorbed  air.  Indeed,  were  It  not  for  these  Interfering  films,  things 
could  scarcely  be  slid  at  all.  Atomlcally-clean  metals  placed  in  contact 
In  a  good  vacuum  adhere  so  strongly  that  they  sieze  and  are  p^^rmanently 
welded.  So  do  bearings  which  are  run  under  such  heavy  loads  that  the 
protective  films  are  stripped  off.  Almost  everything  In  familiar 
experience  Is  protected  by  some  kind  of  film,  so  that  the  observed 
adhesion  Is  very  much  less  than  welding  -  it  Is  usually  more  akin  to 
Van  der  Uaals  forces  than  to  Interatomic  bonding  -  nevertheless,  adhesion 
still  dominates  over  Interlocking  of  asperities. 

Lubrication,  of  course,  is  simply  the  deliberate  interposition 
of  an  Interfering  film  of  oil  or  some  other  liquid.  The  resulting 
reduction  In  friction  owes  something  to  reduced  Interlocking  of  asperities, 
since  the  oil  takes  up  some  space;  but  it  la  mainly  due  to  reduced 
adhesion,  particularly  In  slowly-moving  or  essentially  stationary  systems. 

The  most  useful  quantitative  concept  Is  that  of  the  cn*>fHH.ent 
of  fnatiorit  the  proportionality  constant  between  the  frictional  force 
required  to  slide  a  given  object  and  the  normal  load: 

F  -  /L 

The  coefficient  of  friction  is  a  property  of  the  particular  materials 
involved  but  not  of  the  physical  size  of  the  system.  It  says,  for  example, 
that  the  same  force  is  required  to  slide  a  brick  across  a  table  regardless 
of  whether  the  brick  is  lying  flat,  resting  on  one  edge  or  standing  on 
end.  It  takes  twice  as  much  force  to  slide  a  stack  of  two  bricks  and  ten 
times  as  much  force  to  slide  a  stack  of  ten  bricks.  It  takes  the  same 
force  to  slide  two  stacks  of  five  bricks  each  as  to  slide  one  stack  of  ten 
bricks.  The  coefficient  s,  in  principle,  a  true  constant  and  is 
independent  of  either  the  load  or  the  contacting  surface  area.  Thus,  large 
and  small  objects  have  the  same  coefficients  of  friction. 
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It  is  frequently  stated  that  the  coefficient  of  friction  Is 
also  Independent  of  the  sliding  velocity,  but  this  Is  only  approximately 
true.  It  Is  well  known  that  the  friction  force  required  to  start  sliding 
Is  usually  greater  than  the  force  required  to  maintain  sliding,  and 
this  has  given  rise  to  the  notion  that  the'‘a  are  two  coefficients  of 
friction  -  static  (for  surfaces  at  rest)  and  kinetic  (for  surfaces  in 
motion).  These  are  normally  showr  separately  In  tables  of  friction 
coefficients.  Recent  work  has  shown  that  this  Is  a  gross  oversimplification 
(  1  ).  The  static  coefficient  is  a  function  of  the  time  of  contact,  and 
it  increases  as  the  contacting  asperities  slowly  yield  and  the  real  area 
of  contact  Increases.  The  kinetic  coefficient  varies  with  sliding  velocity, 
but  only  slightly,  usually  by  Just  a  few  percent  as  the  sliding  speed  Is 
raised  by  a  factor  of  ten.  For  most  practical  purposes  of  Interest  here, 
the  kinetic  friction  coefficient  may  be  considered  to  be  a  constant 
Independent  of  the  sliding  velocity. 

To  a  good  approximation,  the  friction  coefficient  Is  Independent 
of  the  roughness  of  the  sliding  surfaces  -  at  least  It  Is  in  the  roughness 
range  normally  encountered  In  engineering  practice.  Very  smooth  surfaces 
give  abnormally  high  friction  coefficients  because  the  area  of  real  contact 
Is  abnormally  high.  Very  rough  surfaces  give  abnormally  high  coefficients 
because  asperity  Interlocking  becomes  excessive.  But  in  the  roughness 
range  In  which  we  actually  find  most  surfaces,  the  friction  coefficient  is 
at  a  minimum  and  almost  Independent  of  roughness. 


(  I  ) 

Another  exception  occurs  when  e  rough  herd  body  slides  on  a  much  softer 
one.  Here  the  asperities  of  the  rough  eurfece  dig  Into  tic  softer 
material  and  the  ordinary  ''friction'*  is  augmented  by  gouging  and  plowing. 
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All  elementary  texts  on  friction  point  out  that  the  coefficient 
of  friction  of  a  given  material  is  approximately  equal  to  the  ratio  of 
the:  material's  shear  strength  and  penetration  hardness 

f  -  a/p 

because  the  sliding  process  involves  the  shearing  of  adhesions  which 
are  limited  in  size  by  the  hardness  of  the  material  (1,2) •  This  says 
that  friction  can  be  reduced  by  decreasing  the  shear  strength  of  the 
interface,  and  of  course  one  of  the  effects  of  lubrication  is  to  interpose 
a  low  shear  strength  interface  between  two  relatively  hard  materials. 

It  also  says  that,  if  one  has  two  different  materials  sliding  on  one 
another;  it  is  only  the  strength  and/or  hardness  of  the  softer  of  the  two 
which  matters. 

The  s/p  relationship  is  neither  exact  nor  fundamental.  It  is 
vitiated  by  involving  assumptions  which  are  only  approximately  true  and 
by  ignoring  effects  such  as  surface  energy  and  asperity  interlocking 
wh'ch  become  important  just  when  ones  gets  into  unusual  regimes  and 
needs  theoretical  guidance  the  most.  Moreover,  s  and  p  are  not  really 
Independent  quantities;  they  are  very  similar  quantities  which  depend  in 
almost  the  same  way  upon  such  factors  as  bond  strength,  nature  of 
dislocations,  etc.,  so  that  one  goes  up  with  the  other  and  the  ratio  is 
quite  similar  for  a  wide  range  of  materials.  For  example,  lead  and  low 
carbon  steel  vary  by  nearly  a  factor  of  100  in  shear  strength  and 
penetration  hardness,  but  /  is  nearly  the  same  for  steel  (1.0)  as  for 
lead  (1.2).  This  point  is  shown  graphically  by  the  following  figure 
which  plots  penetration  hardness  vs.  yield  stress  for  pure  metals  (  1  ). 


Coefficients  of  friction  do  not  vary  g:^atly  from  material 
to  material  as  might  be  expected.  Most  common  metals  exhibit  coefficients 
of  0.1  -  0.3  when  tested  as  one  normally  sees  them.  Most  non-metals 
run  slightly  higher  at  0.3  -  0.4.  Wood-on-steel  iZ  about  0.5,  as  Is 
wood-on-stone,  Iron-on-stone  and  leather-on-lron.  Most  solid  explosives 
run  about  0.6  -  1.1.  Carefully  cleaned  metals  run  about  0.7  ^o  1.3  In 
air'.  This  Is  not  a  large  range  of  variation;  many  physical  properties 
vary  by  orders  of  magnitude  from  material  to  material. 

Wider  variations  are  known.  Teflon,  for  example.  Is  anomalously 
low  at  0.04  due  to  Its  exceptionally  low  surface  energy;  and  freshly- 
prepared  copper  surfaces  In  high- vacuum  can  run  from  5  tc'  ^200  due  to 
cold  welding  on  contact.  Nevertheless,  these  cases  are  wildly  exceptional; 
and  the  fact  remains  that  most  materials  have  about  the  same  coefficients 
of  friction  -  about  0.5. 

To  some  extent,  the  striking  uniformity  of  coefficients  of 
friction  Is  due  to  the  fact  that  most  materials  are  handled  and  tested  In 
air,  and  rather  contaminated  air  at  that.  The  visible  "surface"  of  a  bar 
of  Iron,  for  example.  Is  not  Iron  at  all;  It  Is  Iron  oxide,  covered  with 
a  layer  of  adsorbed  moisture  and  atmospheric  gases  and  most  likely  a  layer 
of  oily  materials  from  adjacent  machinery  and/or  human  beings.  Non-metals 
may  lack  the  oxide  layer  (many  of  then  are  oxides  themselves);  but  they 
will  have  the  moisture,  gassy  and  oily  layers.  Consequently,  the  Interface 
Is  the  same,  or  nearly  the  same,  in  most  cases;  and  one  should  expect  the 
friction  to  be  similar  even  on  dissimilar  materials. 

For  most  ordinary  materials  of  experience  -  probably  Including 
explosives  -  most  of  the  frictional  force  Is  due  to  adhesion.  The 
"roughness  component"  is  small,  approximately  0.05  in  an  overall  coefficient 
of  0.5.  Other  conceivable  effects,  such  as  plowing  and  electrostatic 
attraction,  are  usually  negligible. 

It  should  be  noted  for  the  record  that  there  are  Important 
difficulties  with  the  adhesion  explanation  of  friction;  and  most  modern 
workers  feel  that  while  it  is  essentially  correct,  it  is  also  a  gross 
overslmpllfirntlon  of  the  friction  process.  Nevertheless,  it  serves  so 
v^ll  for  all  but  the  most  exacting  purposes  that  most  workers  in  the  friction 
field  would  rather  amend  it  than  abandon  it. 

The  foregoing  discussion  is  intentionally  simplified  to  serve 
as  a  first  Introduction  for  personnel  who  are  scientifically  trained  but 
unfamiliar  with  this  field.  Those  who  wish  to  delve  deeper  will  find 
excellent  r? sources  listed  in  the  general  bibliography.  The  books  by 
Rabinovlcs  (  I  )  and  Bowden  and  Tabor  (2,3)  are  particularly  recommended, 
as  is  the  NASA  publication  Tntardieoiplimry  Approach  to  Friation  and 
Wear  (  ^  ).  other,  mrre  specific,  references  will  be  found  in  the 
specific  citatluns. 


THE  NATURE  OF  SOLID-SOLID  CONTACT 


To  the  human  eye,  It  appears  that  two  flat  surfaces  placed 
in  contact  actually  touch  each  other  over  their  entire  areas;  but  this 
is  far  from  the  truth.  Visually  smooth  surfaces,  even  optically  polished 
ones,  are  quite  rough  on  a  microscopic  scale.  General  workshop  surfaces 
or.  metals,  for  example,  have  asperities  on  the  order  of  15  microns  high; 
and  even  the  smoothest  metallic  surfaces  have  asperities  100  -  1000 
Angstroms  high.  If  such  a  surface  were  slowly  lowered  onto  one  which  was 
truly  flat  on  the  atomic  scale,  it  would  first  make  contact  only  on  the 
peaks  of  the  highest  asperities,  and  the  remainder  of  the  surface  would 
not  be  in  contact  at  all. 


N \\\\\  \  ^  ^  V  \  \\\V  s  \v  ^  \^\\\\\\\\ 


If  the  material  were  ideally  rigid,  only  point  contact  would 
be  made  on  only  the  three  highest  asperities,  and  the  upper  object  would 
be  supported  on  a  tripod;  bit  in  actuality,  the  asperities  are  so  small 
that  the  weight  of  the  object  far  exceeds  the  yield  strength  of  the 
material,  and  the  asperities  squash  down  and  permit  other,  lower  asperities 
to  contact  also.  These  in  turn  squash  down  and  permit  still  more  contacts 
until  finally  there  is  enough  total  area  of  contact  to  support  the  load 
without  further  yielding. 


The  total  area  of  real  contact  is  only  a  small  percentage  of 
the  total  area  of  the  surface;  most  of  the  surface  is  not  touching  at  all 
and  never  will  touch,  even  if  the  load  is  Increased  enormously.  The  total 
area  of  contact  Is  approximately  equal  to  the  normal  load  L  divided  by 
the  penetration  hardness  p  of  the  material: 

Ay  -  L/p 

As  an  example,  consider  a  smooth,  100-gram  cube  of  copper  resting  on  a 
smooth  steel  plate.  The  loud  is  0.1  Kg,  and  the  penetration  hardness  of. 
cogper  (the  softer  material)  is  6000  Kg/Cm  .  This  gives2Ay,  ■  1.67  X  lO”^ 
Cm‘,  compared  to  a  total  area  of  the  cube  face  of  5.0  Cm^. 


Actual  cases,  of  course,  are  not  as  simple  as  pictured 
above.  The  formula  A  *  L/p  Is  for  purely  plastic  deformation,  and 
Is  an  approximation  anj^iray.  Real  asperities  no  doubt  deform  elartlcally 
over  part  of  their  travel,  and  elastic  deformation  Is  described  by 
the  more  complex  Hertzian  equations  (117),  of  which 


Where: 


A  Real  area  of  contact 
L  “  The  normal  load 

R  =  The  relative  radius  of  curvature  of 
the.two  surfaces 

E'  « 

E  «  Young's  modulus 
V  *  Poisson's  ratio 

Is  a  simplified  example.  Also,  in  real  cases,  both  surfaces  are  rough, 
and  the  contact  pattern  Is  much  more  complex  as  asperities  contact  other 
asperities  on  their  shoulders  or  even  In  the  opposing  valleys. 


Nevertheless,  the  total  area  of  real  oontaot  must  remain  that  given 
by  the  appropriate  deformation  equation,  since  it  is  a  function  only 
of  the  load  and  the  strength  of  the  material,  not  Its  shape.  The  total 
area  of  real  contact  Is  distributed  among  all  the  asperity  junctions. 

The  Individual  contacts  will  range  all  the  way  from  some  Just  barely  touching 
to  others  widely  squashed  out,  and  the  resulting  Junction  sizes  will 
range  all  the  way  from  mere  points  up  to  a  circle  proportional  to  the 
size  of  the  largest  and  highest  asperity. 
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lt  must  not  be  thought  that  asperities  are  anything  like  as 
high  and  jagged  as  the  sketches  above  would  suggest.  The  vertical 
scale  of  the  sketches  is  greatly  exaggerated  for  clarity.  Actual 
asperities  tend  to  be  low  and  broad,  with  contact  points  widely 
spaced.  Early  metallographlc  studies  of  metal  surfaces  showed  that 
solid  metal  surfaces  have  asperities  which  are  typically  10  to  300  uin. 
high;  that  their  slopes  are  shallow,  so  that  their  bases  can  be  50  to 
3000  yin.  across;  and  that  individual  contact  areas  are  about  100  to 
1000  yin.  wide.  The  following  electron  micrographs  of  a  gold  surface 
show  dome-like  asperities.  The  largest  are  roughly  400  yin.  across  and 
50  yin.  high;  the  smallest  are  less  than  20  yin.  across  and  5  yin.  high. 
These  data  are  from  Williamson  (  4  );  other  authors  have  reported  similar 
findings  on  a  wide  variety  of  materials. 


noun  1.— Electron  microcnph  of 
plated  gold  rarface.  <a)  Oentml 
▼lew  at  X  4000.  (b)  A  atofle  con¬ 
tact  arm  X  IBOOO.  t  l  \ 


One  would  like  to  know  the  size  of  each  individual  junction 
and  the  number  of  junctions  of  each  particular  size,  and  much  effort  has 
gone  into  ftudylng  these  parameters  in  real  systems.  People  have  pres.:ied 
opaque  solids  against  glass  plates  and  observed  the  junctions  through 
the  glass  (  3  ),  and  other  people  have  prepared  contour  maps  from  electron 
micrographs  of  aluminum  surfaces  (  A  );  but  the  most  promising  technique 
for  our  purposes  appears  to  be  statistical. 
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Greenwood  and  Williamson  have  shown  (  8  )  that  a  wide 
variety  of  surfaces  have  asperity  height  distributions  which  are 
very  nearly  Gaussian,  and  that  even  in  surfaces  whxch  are  strongly 
non-Gausslan  as  a  whole,  the  upper  50%  usually  Is  accurately  Gaussian 
For  example,  they  report  the  following  data  for  bead-blasted 
aluminum  (  8  ) : 


nmm  4.— OmmIsIIt*  hiigiit  dlsMbattw  pt  baad- 
biuM  aluaUnuk.  Bath  tha  diaMhoUaaa  at  all 
htighia  (.Z)  aad  of  fwak  htit»u  (•)  tfo  Qaaa' 
ilaa.  Tbp  proflla  of  tha  taiw  auiffece  la  thowa 
la  tha  apptr  dlafiam ;  tha  ?trtlcal  nagaMcatlaa 
la  00  tlBMa  the  horlaoatal  aagalOcatlOB.  ^  4  ) 


and  the  data  on  the  following  page  for  abraded  steel  (  8  ): 
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I  «*-  . 

_ I  tTST*  ■ 


•  1  ii  ^ii  w  ii“ 

HE«HT  MOVE  ARMTMRr  OMUM  (nICfohieM 

•••■  ••  IMM  Wl| 


VuitM  5.— Onanlattr*  height  dIatribntlOB  vt  atttd 
■iMl  tpMlmu.  IMitributltm  ti  >11  helghta.  X. 
Dtagflbatloa  of  poke,  •.  Ihlt  qnctmen  wtt 
abndad  on  400  grade  cerbonuidam  piper,  then 
■lid  agilmt  1  copper  block  dooded  with  oleic  ndd. 
It  ipproiUnately  10  kg.  ISO  cin/>  for  SOi  Al¬ 
though  the  distribution  it  Orst  sight  ippeira 
hlidtif  non-OiuMlnu.  In  tact  neirly  00  percent  of 
the  surfice  Is  OnusNlsn:  the  suitacc,  with  in 
srtusl  NtsndHnl  di'vluituu  of  SOn  In.,  would  be* 
hive  In  contact  is  If  Oiuiwlin  with  >  standird 
doTlitlon  of  bilf  tills.  The  protlle  of  the  ssbm 
■urtace  is  shown  In  the  niipcr  dlngram;  the  rertl- 
nl  nragniOratlon  Is  SOO  times  the  horlsontal 
nwgiilOcstlou. 


Willlnmnon  alno  argues  that  moat  surfaces  indeed  ought  to 
be  Gaussian  because  of  the  statistical  principle  that  any  quantity 
which  is  the  result  of  a  large  number  of  random  increments  and 
decrements  will  tend  to  follow  a  Gaussian  distribution  (the  Central 
Limit  Theorum  of  statistical  theory)  (  4  ).  Most  engineering  surfaces 
are  Indeed  the  end  result  of  a  large  number  of  independent  deforming 
processes  (grinding,  machining,  sandblasting,  polishing,  etc.):  ind 
many  other  surfaces  are  made  by  such  processes  as  casting  against 
engineering  surfaces  or  processing  in  equipment  made  from  materials 
having  engineering  surfaces.  Moreover,  the  Central  Limit  Theorum 
argues  that  in  almost  emy  attempt,  human  or  natural,  to  produce 
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unlformlty,  the  residual  errors  will  be  Gaussian  provided  the  attempt 
Is  good  enough  to  exclude  systematic  or  periodic  errors.  And  Williamson's 
findings  seem  to  say  that  even  systematic  errors  leave  a  Gaussian 
distribution  In  the  higher  asperities,  which  are  the  ones  that 
participate  In  solid-solid  contacts. 

So  It  seems  likely  that  we  can  arrive  at  a  useful  atatiatiaal 
description  of  asperity  and  Junction  sizes  even  If  we  cannot  determine 
the  sizes  of  Individual  Junctions  during  the  sliding  process.  Inasmuch 
as  the  frictional  Initiation  of  explosives  Is  Itself  a  statistical 
phenomenon,  a  statistical  description  of  the  Junctions  may  well  be 
quite  good  enough. 


A  number  of  Investigators  have  attacked  the  problem  of 
deducing  the  statistics  of  Junction  sizes  from  statistical  descriptions 
of  asperity  heights.  Greenwood,  for  example,  has  shown  (  9  )  that 
the  average  Junction  size  does  not  change  with  the  applied  load  If 
the  asperity  size  distribution  follows  any  kind  of  exponential  law. 

The  original  Junctions  grow,  of  course;  but  new,  tiny  ones  form  at 
Just  the  rate  needed  to  keep  the  average  constant  and  the  distribution 
essentially  normal.  Tabor  has  pointed  out  in  a  comment  on  Greenwood's 
paper  (  9  )  that  Gaussian  distributions,  which  as  we  have  seen  we 
expect  to  find  on  real  surfaces,  are  close  enough  approximations  to 
exponential  distributions  to  make  the  arguements  valid  for  them  as  well. 

We  already  know  the  total  junction  area  -  it  is  the  area  of 
real  contact.  If  we  could  determine  the  total  number  of  Junctions, 
we  would  then  have  the  average  Junction  area;  and  Greenwood's  statistics 
would  tell  us  the  entire  distribution  of  Junction  sizes. 


Rabinowlcz  (  1  )  has  investigated  the  number  of  Junctions 
existing  during  sliding  by  several  techniques,  including  some  statistical 
ones.  He  argues  that  the  real  area  of  contact  (the  sum  of  the  asperity 
Junctions)  is  constant  as  long  as  the  load  remains  constant,  but  that 
individual  Junctions  are  Csintinually  made  and  broken  in  a  random 

fashion.  This  gives  rise  to  random  variations  in  the  instantaneous 
(Coefficient  of  friction  over  a  sliding  distance  which  is  a  function  of 
the  average  Junction  size.  Rabinovicz  derives  the  diameter  of  the 
average  Junction  from  an  auCocorralatlon  analysis  in  which  the  friction 
coefficients  •’f®  measured  at  Intervals  off  a  friction- 

distance  pilot,  and  Ihe  autocorrelation  coefficient 


n 


I?  Vf  -  -  7) 


n  -  ic 


tlifi  -  7)^ 


is  calculated  (  22  }.  The  distance  k  at  which  the  autocorrelation  drops 
tc  zero  is  equal  to  twice  the  average  Junction  dimMter.  For  copper  on 
steel,  he  finds  an  average  Junction  diameter  of  9  X  10*^  cm.,  which  is  in 
reasonable  agreement  with  values  derived  from  measurements  of  wear  debris. 


Fig.  3.13.  AutooomUtion  ooefGoient  of  friction  trncea  m  n  function  of  diaUnco 
niong  the  iliding  tmck.  Since  re  drope  to  loro  in  18  X  10~*  cm,  the  average 
junction  diameter  ia  9  X  10~*  cm.  /  i  \ 
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Rablnowlcz  also  estimates  the  total  number  of  junctions  from 
a  different  treatment  of  the  same  statistical  data  (  1  ).  He  points  out 
that  If  an  assumption  Is  made  as  to  the  range  In  strength  between  the 
weakest  and  strongest  junctions,  and  a  factor  of  2  seems  plausible  In 
many  cases.  It  Is  possible  to  estimate  from  the  amplitude  fluctuations 
how  many  junctions  must  have  been  present  at  any  time.  He  derives  the 
expression 


— 

f  *  A/n 
Where : 

0  •  standard  deviation  of  friction  values 
f  •  the  mean  friction  coefficient 
n  •  the  number  of  junctions 

Rablnowlcz  also  deduces  the  size  distribution  of  junctions 
from  autocorrelation  data  by  using  all  the  autocorrelation  coefficients, 
not  just  the  one  at  which  the  correlation  reaches  zero;  and  he  finds 
size  distributions  In  good  agreement  with  those  deduced  from  wear 
particle  measurements  (  22). 

Ling  has  also  dene  calculations  In  this  area  (  13),  although 
his  work  has  been  more  concerned  with  the  deformations  of  asperities 
during  sliding  and  Is  better  considered  In  the  next  section  under  the 
generation  of  heat  between  sliding  solids.  In  general,  nis  results 
are  consistent  with  those  discussed  above. 

Jones,  et.al. ,  <  10)  and  Vhitehouse  and  Archard  (  20)  have 
also  recently  published  reviews  of  methods  of  obtaining  areas  of  real 
contact  and  of  the  properties  of  random  surfaces  In  contact.  Their 
conclusions  are  essentially  those  already  stated,  starting  from  slightly 
different  initial  assumptions. 

Tsukisoe  and  Hisakado,  In  a  recent  trilogy  of  papers  (16,  17,  18  ), 
present  detailed  calculations  of  the  separation,  the  real  area  of  contact, 
the  number  of  contact  points,  the  average  radius  of  the  contact  points, 
and  the  distribution  of  radii  of  the  contact  points  between  two  metal 
surfaces.  They  assume  that  the  surfaces  contain  a  large  number  of 
asperities  in  the  form  of  cones  of  equal  base  angles,  that  the  heights 
of  the  cones  are  Gaussian,  and  that  the  metal  deformation  is  plastic. 

They  ..iso  report  mplrical  measursnents  on  aluminum  surfaces,  with  good 
agreement  with  the  calculated  values.  These  papers  are  too  mathematical 
for  review  in  this  introductory  discussion,  but  they  deserve  the  most 
careful  study  by  anyone  who  may  wish  to  pursue  this  line. 
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In  summation:  although  there  are  still  a  number  of  unresolved 
problems,  there  is  general  agreement  that  real  surfaces  consist  of  arrays 
of  asperities  with  a  Gaussian  distribution  of  peak  heights,  that  such 
surfaces  contact  each  other  only  at  junctions  of  aggregate  area  much 
less  than  the  total  apparent  area  of  the  surfaces,  and  that  there  are 
techniques  for  calculating  the  number  and  the  size  distribution  of  the 
junctions. 


There  are  some  fairly  dubious  assumptions  built  into  some 
of  these  derivations,  for  example  the  assumption  of  purely  plastic 
deformation  of  asperities;  but  other  mathematical  analyses  have  shown 
that  It  really  does  not  matter  If  certain  other  assumptions,  such  as 
an  exponential  distribution  of  asperity  heights,  are  true  (  9  ). 

There  are  also  some  contradictory  assumptions  In  some  of  the  complementary 
arguments;  for  example  Rabinowlcz  assumes  all  junctions  the  same  size 
in  his  autocorrelation  treatment  (22  )  which  was  used  above  to  derive  a 
junction  size  distribution;  but  the  invocation  of  Gaussian  distributions 
can  rationalize  a  lot  of  that  kind  of  thing  too. 


All  the  models  have  been  criticized  on  the  basis  that  real 
surfaces  do  not  consist  of  Idealized  cones  or  wedges  or  sphere  sections; 
but  it  has  been  shorn  that  It  really  does  not  make  much  difference  what 
the  asperity  shape  is  or  what  the  exact  deformation  law  is  -  one  still 
gets  an  acceptably-Gausslan  distribution  of  junction  sizes,  and  one  gets 
reasonable  agreement  with  oitpirical  data  in  test  cases. 

Nothing  is  yet  known  about  the  topography  of  the  surfaces  of 
explosive  blocks,  and  some  measursMats  should  be  early  on  the  agenda 
of  any  progras  to  adapt  these  insights  to  the  study  of  explosives;  but 
It  would  be  surprising  to  find  that  explosive  surfaces  did  not  consist 
of  asperities  Gaussianly  distributed. 


the  generation  of  heat  between  sliding  solids 


The  generation  of  heat  between  sliding  solids  has  been 
familiar  since  the  first  caveman  slid  down  a  vine  and  burned  his 
hands.  Boy  Scouts  use  It  when  they  rub  two  sticks  together  and 
start  their  campfire.  We  contend  with  it  dally  in  the  bearings  of 
all  our  machines.  In  fact,  in  many  machines  such  as  auto*T!otive  or 
helicopter  transmissit.ns,  only  an  Infinitesimal  amount  of  the  oil 
used  is  actually  needed  to  lubricate  the  gears  -  all  the  rest  Is 
required  solely  as  a  coolant  to  carry  away  the  heat  generated  by 
the  friction  between  the  gears.  If  the  heat  is  not  removed,  the 
result  can  quickly  be  disastrous  for  the  gears.  In  the  case  of 
sliding  explosives,  failure  to  deal  with  frictional  heat  can  be 
even  more  dramatically  disastrous. 

It  Is  easy  to  determine  the  heat  generated  In  a  simple 
sliding  situation.  It  Is  simply  the  thermal  equivalent  of  the  %rork 
expended  In  making  the  object  slide  against  friction:  one  calorie 
for  each  4.18  X  10^  ergs.  This  datum  Is  seldom  easy  to  obtain  In 
a  hazard  situation  where  one  might  have  a  shower  of  complex  shapes 
cascading  down  under  gravity ,  but  It  Is  easy  to  obtain  in  a  laboratory 
situation  where  one  can  set  up  simple,  measureable  operations  and 
collect  quantitative  data. 

The  gross  amount  of  heat  Is  not  the  datum  of  interest  for 
an  explosive,  though;  one  wants  to  know  the  resulting  tgmperatxa^  In 
a  given  volume  of  explosive  and  the  eine  of  the  volume  which  reaches 
the  given  temperature. 

The  Interface  temperature  can  be  estimated  from  conventional 
heat  flow  equations  by  assuming  that  the  heat  Is  generated  in  the 
interface  and  conducted  away  into  each  of  the  two  sliding  objects 
according  to  their  respective  thermal  conductivities.  For  the  case 
of  an  infinitely  long  cylinder  rubbing  end-wise  on  a  flat  plate  at 
moderate  sliding  speeds: 

fLO 

’  4Jr(rr+  kil  (Equation  A) 

Where: 

AT  «  Kean  temperature  rise  above  the  rest  of  the  material 
f  ■  The  coefficient  of  friction 
L  •  Th..  aomal  load  (in  force  units) 

J  >  The  mechanical  equivalent  of  heat 
r  ■  The  radius  of  the  rubbing  cylinder 
and  k\  and  k2  are  the  thermal  conductivities  of  the 
two  contacting  surfaces. 
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The  simple  equation  above  needs  many  corrections  and 
elaborations  for  real  situations.  At  high  sliding  speeds,  far  more 
heat  flows  into  the  plate,  which  continuously  sends  new,  cool  material 
into  the  zone  where  heat  is  generated,  than  into  the  cylinder,  the 
sams  part  of  which  always  remains  at  the  interface.  This  problem  has 
been  attacked  by  m  any  investigators,  including  Blok  (  31  ) ,  Jaeger  (40), 
Bowden  and  Thomas  (33),  and  Archard  (32). 

Francis,  in  a  forthcoming  paper  (41),  derives  an  analytic 
expression  for  the  steady-state  interfacial  temperature  field  in  a 
sliding  circular  Hertzian  contact,  taking  into  account  the  difference 
in  bulk  temperatures  of  the  two  bodies  and  the  ellipsoidal  distribution 
of  the  frictional  power  In  the  contact  area: 


He  derlvea  a  formula  for  the  mntmm  interfacial  temperature  "  in 
terms  of  the  total  frictional  power  Q,  the  radius  of  the  contact*#, 
the  thermal  conductivity  k,  and  a  diaen:  lonlesa  parameter  B  which 
deacribaa  the  power  distribution. 
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Francis'  paper,  in  the  form  of  a  preprint  of  his  presentation  to  the 
25th  ASLE  Annual  Meeting  in  Chicago,  May  A-8,  1970,  is  inclvried  in  its 
entirety  in  the  appendix  because  it  '•■i  not  yet  readily  available  from 
standard  sources. 

For  hot-spot  calculations,  the  radius  r  and  the  load  L  of 
equation  (  A)  should  be  those  for  a  single  asperity  junction,  and  their 
sums  should  equal  the  area  of  real  contact  and  the  total  load,  respectively. 
If  one  knew  the  total  number  of  junctions  and  the  size  of  each,  one  could 
pro-rate  the  load  on  each  and  the  amount  of  heat  flowing  into  each,  and 
calculate  the  resulting  temperature  rise  at  each.  One  in  general  does 
not  knov  much  about  individual  junctions,  but  some  approximations  are 
possible. 

Rablnowlcz  (1)  shows  that  the  radius  of  a  junction  is  related  to 
the  surface  energy  y  and  the  hardness  p  of  the  softer  material  by  the 
order-of-magnitude  relation: 

p  =  12,000  ^ 

P 

and  that  the  load  L  an  asperity  will  catty  is  its  area  times  its  hardness: 


L  =  vp^p 

Substituting  ..n  Aquation  (A  ),  he  obtains  the  final  expression: 


AT  -  9400/yP 
J(ki  +  ^2) 


and  calculates  flash  temperatures  on  individual  junctions  as  follows: 


Material _ 

Brass  on  brass 
Steel  on  steel 
Bakellte/Bakellte 
Glass  on  glass 

The  approximations  Involved  are  rather  gross,  and  they  Include  the 
assumption  that  the  total  load  is  not  e  factor;  but  they  give  an  idea 
of  the  magnitudes  obtainable. 


JL-  _JL_ 
0.4  900  0.26 
0.5  1500  0.11 
0.3  100  0.0015 
0.9  500  0.0007 


[C/cm/sec) 


0.15 

0.75 

2.2 

70.6 
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One  could  also  plug  in  a  Gaussian  distribution  function,  as 
discussed  in  the  preceedlng  section,  and  obtain  a  statistical 
description  of  the  entire  array  of  junctions.  This  should  be  quite 
adequate  for  the  explosive  situation.,  since  explosive  initiation  is 
also  statistical  in  nature.  (In  fact,  the  statistical  nature  of 
explosions  is  one  of  the  strong  arguments  for  Gaussian  distributions 
of  asperities.) 

Ling  has  done  a  good  bit  of  this  (  36  ) »  and  has  derived 
expressions  for  the  instantaneous  temperature  of  an  asperity  junction 
using  a  stochastic  model  of  junction  formation  and  breaking  during 
the  sliding  process.  At  a  sliding  velocity  of  5000  ft/min,  he 
calculates  flash  temperatures  of  around  1800  °F  for  steel-on-steel, 
with  variations  from  approximately  500®F  to  2000®F. 

For  explosives,  still  another  factor  becomes  important  in 
attempting  to  calculate  the  temperature  of  a  frictional  hot-spot:  the 
explosive  will  begin  to  decompose  quite  exothermically  at  some  elevated 
temperature,  and  the  temperature  will  run  away,  perhaps  explosively. 
Fortunately,  this  process  at  least  has  been  well  studied,  and  expressions 
are  available  to  calculate  the  thermal  evolution  of  a  hot-spot  as  functions 
of  both  the  hot-spot  size  and  temperature.  For  example,  Zinn  ( 48  ) 
expresses  the  progress  of  an  exothermic  reaction  proceeding  by  first- 
order  kinetics  in  a  medium  of  thermal  diffusivity  k  by: 

3T/3t  =  +  (QZ/C)we-fi'/^ 

where  the  symbols  have  the  usual  meanings.  A  more  useful  expression 
is  the  Zinn-Mader  (49)  formulation  for  the  minimum  explosion  temperature: 

t  »  J^/(E/T  -  E/T.) 
exp .a  m  1 

where: 

t  >  time  to  explosion 

■  sample  thickness 
a  *  thermal  diffusivity  ■ 

E  *  activation  energy 

T  >  minimum  temperature 

for  explosion 
«  surface  temperature 

The  function  f  is  a  complex  one  best  handled  graphically: 


.\  »  thermal  conductivity 
p  *  density 
c  “  specific  heat 
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It  still  remains  for  someone  to. Integrate  all  these  approaches 
into  a  coherent,  quantitative  description  of  the  formation  of  hot-spots 
between  sliding  blocks  of  explosives.  Expressions  are  available  for  the 
amount  .>f  heat  generated,  for  the  ntuober,  size  and  distribution  of 
asperity  junctions,  for  the  temperature  profiles  in  and  between  the 
junctions  and  their  variation  with  time,  for  the  flow  of  heat  through 
the  junctions  and  into  the  bulk  of  the  explosive  blocks,  and  for  the 
thermal  evolution  of  hot-spots  in  the  explosive.  This  would  be  an 
ambitic.t-s  undertaking,  and  its  execution  is  far  beyond  the  scope  of 
this  study. 


Phnsiaal  mechanisms  for  heat  generation 

The  foregoing  discussion  has  dealt  with  the  fate  of  heat 
generated  between  sliding  solids  to  the  neglect  of  the  actual  generation 
of  that  heat.  One  does  not  need  to  know  the  origin  of  heat  in  order  to 
understand  its  flow  and  effects,  but  some  consideration  of  the  mechanisms 
of  heat  generation  is  useful  in  understanding  friction. 


Viscous  flow  can  take  place  in  either  a  liquid  or  a  solid,  and 
the  resultant  heating  is  directly  proportional  to  the  work  expended  in 
causing  the  flow.  For  example,  for  shear  in  a  thin  film: 

H  ■  uAU^/h 

where : 

H  •  work  per  unit  time  dissipated  into  heat 
y  ■  viscosity  of  the  flowing  material 
A  *  area  of  the  sheared  film 
U  ■  sliding  speed 
h  ■  thickness  of  the  sheared  film. 

For  a  typical  oil  in  a  reasonable  machine  bearing,  one  can  plug  in 
appropriate  values  of  3.3  X  10~^  for  y,  25  Sq.In.  for  A,  300  in/sec 
for  U,  and  lO'^  in.  for  h,  and  obtain  a  value  of  7430  in-lb/sec,  or 
1.125  HP,  or  840  watts  of  heating. 

If  all  the  heat  remained  in  the  oil,  and  the  oil  stayed  in 
the  clearance  space,  the  temperature  would  rise  at  the  rate  of  H/Vc, 
where  V  >  the  volume  of  the  film  and  c  •  the  heat  capacity  of  the 
oil  per  unit  volume.  Taking  c  *  140  ln-lb/cu.in.”F  for  a  typical 
petroleum  oil,  one  calculates  2120*F/sec. 


Such  drastic  heating  rates  do  not  occur  In  normal  practice 
because  heat  Is  lost.  Into  the  massive  heat  sink  of  the  metal  parts  and 
because  hot  oil  Is  continuously  extruded  and  fresh,  cool  oil  Is 
continuously  supplied.  In  cases  of  lubrication  failure,  the  reality 
of  the  above  numbers  Is  quickly  apparent  as  bearings  overheat  and 
sleze.  One  observes  actual  melting  of  even  hardened  steel  surfaces. 

Such  heating  rates  may  also  be  entirely  realistic  In  the  case 
of  a  sliding  explosive,  where  the  "lubricant"  Is  molten  explosive  or 
even  the  surface  layers  of  the  solid  Itself,  the  shearing  rate  may  be 
extremely  high  and  the  heat  conductivity  of  the  solid  Is  low. 

Calculation  of  heating  for  cases  other  than  that  of  a 
well-defined,  lubricated  bearing  Is  difficult,  because  It  Is  seldom 
possible  to  pin  down  the  properties  or  the  dimensions  of  the  material 
being  sheared.  In  fact,  even  the  calculation  for  the  bearing  Is  only 
grossly  correct.  It  gives  the  averaged  temperature  for  the  entire 
volume  of  the  oil  film;  but  the  oil  film  Is  not  sheared  uniformly. 

It  contains  a  shear  gradient  and  a  consequent  temperature  gradient. 

The  local  heating  In  zones  of  highest  shear  can  be  very  much  greater 
than  the  overall  heating. 

Non-unlformlty  of  heating  has  been  studied  by  Appeldoorn, 
et.aZ. ,  for  fluids  forced  through  a  capillary  (46  >.  For  adiabatic. 
Incompressible  flow,  the  bulk  heating  Is  a  function  of  the  driving 
pressure,  and  the  specific  heat  and  the  density  of  the  fluid: 


41^  -  P/C,pJ 
where : 

AT,  ■  the  bulk  temperature  rise 
P  "  -  the  driving  pressure 
C  *  the  specific  heat  of  the  fluid 
-  the  density  of  the  fluid 
J  -  the  mechanical  equivalent  of  heat 

61^  for  a  typical  hydrocarbon  oil  is  about  7*F  for  each  1000  psi. 

But  this  Is  only  an  average  value;  the  actual  temperature 
profile  is  very  non-uniform  across  the  capillary  because  the  shear 
is  very  non-uniform.  It  is  well  known  that  the  flow  profile  across 
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a  pipe  is  approximately  (not  exactly)  parabolic,  with  the  fluid  at 
the  wall  essentially  motionless  and  that  at  the  center  essentially 
not  in  shear. 


(Arrows  indicate  velocity  vectors) 


Appeldoorn,  et.al.  ^  have  shown  experimentally  (46,  47)  that  this  flow 
(and  shear)  profile  Is  reflected  in  a  temperature  profile  which  is 
similar: 


and  they  have  measured  temperature  rises  at  the  wall  as  high  as  140*F 
for  the  case  where  the  bulk  rise  was  only  7*F. 

For  the  case  of  one  solid  sliding  upon  another  without 
lubrication,  the  shear  gradient  may  be  very  steep  Indeed;  and  the  area 
of  shear  may  be  very  small,  i.e. :  the  area  of  real  contact,  which 
Is  vastly  smaller  than  the  bulk  area  of  contact;  so  that  the  local 
hot-spots  due  to  shear  may  be  very  hot  Indeed. 

Nothing  limits  the  above  phenomena  to  liquids.  Solids  can 
flow  too,  under  sufficient  pressure;  and  sufficient  pressure  specliically 
exists  in  an  asperity  Junction  already  loaded  to  the  plastic  flow  point. 
The  following  diagram  which  has  been  used  repeat.idly  by  Bowden  (  3  ) 
illustrates  the  point: 
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Two  solids  placed  in  contact  are  supported  on 
the  summits  of  surface  irregularities.  The  pressure  exceeds 
the  yield  pressure  />„,  of  the  solid,  u  hich  Jlows  plastically  until 
the  area  of  contact  is  sujFicient  to  support  the  load  W. 

Hence  .“I  =  The  material  around  the  regions  of  \ 

contact  will  be  elastically  deformed.  If  s  is  the  shear  strength 
of  the  junctions,  F  =  As. 


Solid  flow  also  exists  whenever  a  high  asperity  of  a  haifd 
surface  plows  a  furrow  in  a  softer  opposing  surface.  In  addition  to 
cutting  shavings,  such  plowing  usually  leaves  a  ridge  on  either 
side  of  the  furrow,  and  this  has  to  be  due  to  the  plastic  flow  of 
solid  material.  The  pressures  in  psi  are  tremendous;  so  are  the  shear 
gradients,  and  so  are  the  local  temperatures. 

Note  also  that  the  melting  point  of  a  rubbing  material  does 
not  limit  the  temperature  riae,  as  is  so  often  stated.  If  the  pressure 
and  the  sliding  speed  are  maintained,  the  molten  material  will  be  in 
severe  shear  and  will  be  heated  further  by  viscous  flow.  This  is  the 
reason  why  organic  explosives  which  catmot  be  ignited  by  slow  sliding 
pan  be  ignited  by  fast  sliding;  fast  sliding  forces  the  molten  layer 
Into  severe  shear.  It  is  also  the  reason  why  grit  Increases  the  friction 
sensitivity  of  an  explosive:  a  melting  asperity  collapses  and  reduces 
the  pressure  somewhat,  while  a  rigid  grit  particle  ?nalntalns  the  local 
pressure  and  forces  the  molten  material  Into  severe  shear  through  narrow 
clearances. 

The  phenosMna  of  viscous  and  plastic  flow  probably  account 
for  most  of  the  heat  produced  in  solid  wad  in  lubricated  sliding.  The 
breaking  of  asperity  adhesions  may  involve  other  phenomena  such  as 
surface  energy  effects  (see  Rabinowict,  chapter  6);  but  this  too  may 
be  accountable  by  the  heat  of  the  plastic  de/omation  of  the  asperity 
just  prior  to  rupture. 

Elastic  hysteresis  has  been  mentioned  as  a  possible  source  of 
heating  aa  an  asperity  deforms  elastically  under  load  and  then  recovers 
its  Initial  shape  with  release  of  Its  stored  energy.  In  most  cases,  most 
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of  the  stored  energy  probably  goes  back  into  the  deforming  asperity 
to  boost  it  on  its  way  as  it  departs,  but  at  least  one  case  has  been 
reported  (  45)  of  glass  particles  becoming  incandescent  apparently 
after  the  deforming  slider  had  departed.  In  viscoelastic  materials 
such  as  rubber  -  and  perhaps  some  PBX  -  there  is  a  bulk  hysteresis 
loss  which  is  significant.  In  fact,  it  is  the  major  mode  of  energy 
loss  in  a  rubber  tire  on  a  wet  or  icy  pavement,  and  it  accounts  for 
nearly  all  of  an  automobile's  traction  in  such  cases  (68,  69,  70  ). 

As  every  motorist  knows,  the  effect  is  very  small  compared  to  the 
normal  dry  adhesion  of  a  tire  to  pavement. 

It  should  also  be  noted  that  "elastic"  '*s  a  deceptive  word. 

Few,  if  any,  materials  are  perfectly  elastic.  A  repeatedly  flexed 
spring  will  get  hot,  not  from  elastically  stored  energy  which  is  given 
back  to  the  deforming  source,  but  from  hysteresis  losses  which  are  really 
flow  losses.  So  will  a  rubber  tire.  Nevertheless,  the  heating  is 
still  real;  and  hysteresis  losses  may  make  their  contribution  in  some 
cases.  Most  authorities,  though,  feel  that  it  is  a  minor  contribution 
compared  to  plastic  and  viscous  flow. 


There  are  two  other  Important  phenomena  associated  with 
sliding  and  grinding  processes:  fracture  and  electrostatic  discharges. 
They  are  both  outside  the  scope  of  this  study,  and  they  are  probably 
of  minor  consequence  in  most  frictional  initiations;  but  mention  of 
them  needs  to  be  made. 

Pox  and  Soria-Ruli,  of  the  Cavendish  Laboratory,  have 
studied  the  "apparent  temperatures"  generated  in  the  leading  edge  of 
a  crack  propagating  in  an  explosive,  and  have  found  them  to  reach 
values  as  high  as  several  thousand  degrees  (81  ) .  These  are  not  actual 
temperatures:  they  are  those  temperatures  one  calculates  from  the 
observed  rate  data  using  kinetic  parsMters  derived  fron  more  ordinary 
experiments,  and  they  may  or  may  not  c  rreepond  to  thermal  temperatures. 
Nevertheless,  they  are  moat  provocative;  and  they  suggest  that  surface 
energies  may  play  significant  roles  in  brittle  explosives. 

It  has  also  been  reported  (118)  that  the  work  required  to 
form  one  square  centimeter  of  new  surface  on  ^^2^3  highly  dependent 
upon  the  method  used: 


Method 

Eras 

Pulling 

500 

Zero  creep 

700 

Dry  crushing 

100,000 

Crushing  under 

liquid 

75,000 

Clearly,  crushing  requires  the  expenditure  of  much  more  work  than  is 
recovered  in  the  formation  of  new  surface.  The  rest  must  appear  as 
heat.  The  25%  reduction  by  wetting  is  also  provocative,  and  has 
similar  implications. 

At  this  writing.  Fox  has  an  active  program  on  the  subject,  and 
he  has  promised  a  paper  specifically  directed  toward  fracture  initiation 
of  explosives..  It  is  recommended  that  interested  investigators  await 
Fox's  paper.  This  study  contains  no  further  reference  to  fracture  except 
a  bibliographic  listing  of  publications  which  were  encountered  in  the 
field  in  the  course  of  the  literature  search. 

Electrostatic  discharges  are  frequent  concomitants  of  sliding, 
and  they  may  be  important  initiation  sources,  particularly  for  primary 
explosives:  but  they  are  not  "friction",  and  they  are  beyond  the  scope 
of  this  study.  Piezoelectricity  and  the  emission  of  exoelectrons  from 
freshly  broken  or  abraded  surfaces  are  closely  relate  a  phenomena  which 
may  also  be  important,  but  they  too  are  bmyond  the  scope  of  this  study. 
The  entire  subject  of  electrical  friction  phenomena  is  a  large  one,  and 
could  easily  absorb  a  study  larger  than  thlc  one.  The  bibliography 
contains  a  number  of  excellent  starting  places  for  suywne  interested  in 
pursuing  it. 
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4.  FRICTION  HAZARDS  TO  EXPLOSIVES 


This  study  has  found  no  magic  for  reducing  friction  or 
friction  hazards  -  except  for  the  ancient  magic  of  lubrication  - 
and  the  study  and  control  of  friction  hazards  remains  very  much  a 
Black  Art.  For  tne  time  being  at  least,  the  testing  and  the  safe 
handling  of  explosives  remains  highly  empirical;  and  innovations 
should  be  made  only  with  the  greatest  caution. 


LABORATORY  FRICTION  SENSITIVITY  TESTERS 

There  is  a  bewildering  variety  of  friction  sensitivity 
testers,  and  the  number  of  variations  is  approximately  equal  to  the 
number  of  active  laboratories.  The  testers  are  all  highly  empirical, 
and  their  value  is  usually  limited  to  a  particular  application  in  a 
particular  time  and  place;  they  do  not  yield  fundamental  data.  They 
are  only  partially  successful  even  in  yielding  empirical  data,  partly 
because  the  hazardous  stimulus  is  usually  quite  difficult  to  identify 
and  simulate,  but  mostly  because  the  essential  nature  of  friction  is 
only  dimly  understood.  In  these  circumstances,  it  is  neither  practical 
nor  orofitable  to  review  each  tester  in  detail:  it  seems  better  to 
review  the  principles  of  operation  and  critique  them.  Detailed 
descriptions  of  individual  testers  may  be  found  in  the  literature 
through  the  attached  bibliography. 

Host  friction  sensitivity  tasters  are  attempts  to  produce 
a  realistic,  usually  mall-scale,  simulation  of  actual  processing 
operations  or  accident  situations:  and  they  can  be  grouped  into  three 
categories: 

(1)  Those  which  shear  a  thin  layer  of  explosive  between  two 
rigid  plates  of  steel  or  other  material  of  construction. 

Some  of  these  are  lineal  and  single-pass,  and  sos«  are 
rotary  and  continuous. 

(2)  Those  in  which  a  block  of  explosive  is  rubbed  violently  on  a 
h.rd  or  abrasive  surface. 

(3)  Those  in  which  a  sample  of  explosive  is  subjected  to  extreme 
defonuition  in  an  imiMce  cr  extrusion  event. 


*  The  word  "empirical"  is  not  used  in  uerogation.  The  writer 
is  an  enthusiastic  e^iricist,  on  the  grounds  that  empiric im  is 
necessary  to  keep  theoreticians  honest.  But  it  is  also  true  that 
sound  theory  is  necessary  to  keep  smpAriclsts  relevant*  and  this  study 
attempts  to  emphasise  fuimiamentals. 


The  Shear  testa  are  the  most  popular,  and  they  are  represented 
by  the  (British)  Explosives  Research  and  Development  Establishment 
Mallet  Friction  Test  (  51),  the  E.R.D.E.  Emery  Paper  Friction  Test  for 
Sensitive  Explosives  (51  ),  the  E.R.D.E.  Sliding  Block  Friction  Test  (51  ), 
the  U.S.  Bureau  of  Mines  Pendulum  Frlctica  Test  (  51),  th^  USNOTS 
Pendulum  Friction  Test  (  51),  the  Picatlnny  Arsenal  Pendulum  Friction 
Test  (51  ),  the  Allegany  Ballistics  Laboratory  Sliding  Friction 
Machine  (56,  57,  58  ),  the  Esso  Friction  Screw  (66  ),  and  others. 

All  the  Shear  tests  squeeze  the  sample  between  two  tool  faces 
of  steel  or  other  material  of  construction  and  subject  It  to  shear  by 
means  of  a  sliding  motion  of  the  two  tools.  They  use  varying  degrees 
of  pressure  from  weights,  springs  or  hydraulic  cylinders;  and  they  use 
varying  sliding  speeds.  The  pinching  surfaces  can  be  steel,  fiber, 
aluminum,  glass,  or  any  other  material  the  investigator  thinks  realistic 
of  a  potential  real  event.  The  surfaces  can  be  smooth  or  rough,  and 
they  can  incorporate  added  grit  of  varying  hardnesses. 

Most  of  the  Shear  tests  are  Intended  to  be  simulative  of 
process  or  accident  hazards,  such  as  the  pinching  of  a  sample  between 
the  blade  and  the  wall  of  a  sigma  blade  mixer,  or  the  pinching  of  a 
sample  between  the  threads  of  a  bolt  and  Its  bolt-hole.  The  ABL 
Sliding  Friction  Test  (56,  57,  58  )  is  probably  the  best  of  these,  and 
it  has  been  widely  copied  in  this  country. 

The  common  weakness  in  all  the  Shear  tests  is  that  they  have 
little  to  do  with  the  friction  properties  of  the  explosive  itself;  the 
friction  is  mainly  between  the  tvo  plates,  or  the  mallet  and  the  anvil, 
with  the  explo8i'’e  playing  a  minor  role  as  a  lubricant  or  as  an  acceptor 
of  heat  generated  by  the  rubbing  of  the  two  plates.  The  explosive  is 
doubtless  sheared  like  a  viscous  liquid  to  some  extent,  but  it  is  never 
clear  to  what  extent*  The  tests  are  simulative  of  real  pinching  events, 
and  they  do  furnish  empirical  estimates  of  degree  of  hazard,  which  is 
their  intention;  but  they  are  not  susceptible  of  detailed  analysis,  and 
they  cannot  give  fundMentsl  insights.  They  coid>ine  too  many,  too 
complex,  phenomena. 

One  Shear  test,  the  Esso  Friction  Screw  (66  ),  was  developed 
ss  a  purely  empirical  test  to  correlate  with  accident  experience  in  the 
laboratory  handling  of  NF  compounds;  but  It  has  promise  of  yielding 
fundamental  data,  given  some  further  development.  It  subjects  the 
sample  to  slow  shear  between  two  plates  under  ever-increasing  pressure 
and  In  the  presence  of  added  grit*  The  main  friction  is  emphatically 
not  that  between  the  two  plates;  because  it  makes  no  difference  whether 
the  plates  are  steel,  aluminum  or  Teflon  In  most  cases,  and  it  makes 
little  difference  how  rapidly  the  screw  is  turned*  The  overriding 
independent  variable  is  the  hardness  (not  the  melting  point)  of  the 
added  grit. 
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On«  can  speculate  that  the  grit  particles  act  as  the  contacting 
asperities  between  the  two  plates,  and  that  the  explosive  acts  as  the 
lubricant  or  separating  film  between  asperities.  Pressures  are  such  as 
to  crush  the  grit  particles,  so  that  the  grit  hardness  (>  crush  strength) 
controls  the  asperity-asperity  contact  pressure.  Sliding  speeds  are  so 
low  as  to  preclude  gross  bulk  heating,  so  that  the  temperatures 
generated  must  be  simply  the  Flash  Temperatures  on  the  asperlt:  tips. 

The  flash  temperatures  are  probably  not  calculable  from  the  dry  friction 
models  discussed  in  Section  3,  since  the  contacts  are  lubricated  by  the 
sample;  but  they  may  be  calculable  from  models  currently  being  developed 
at  Esso  Research  and  Engineering  Company  for  Boundary  Lubrication.  This 
work  is  incomplete  at  this  writing;  but  it  will  be  published  during  1971 
(  80),  and  a  preliminary  version  is  given  in  the  Appendix  of  this  study. 

Violent  rubbing  of  a  block  of  explosive  on  a  hard  external 
surface  does  Involve  the  frictional  properties  of  the  explosive  itself; 
because  the  explosive  is  the  softer  material,  and  its  shear  strength 
and  penetration  hardness  will  be  the  ones  of  Importance.  Such  tests 
are  typified  by  the  Pante.'.  Skid  Test  (51  )  and  the  A.U.R.E.  Oblique 
Impact  Test  (  51).  The  Canadian  Friction- Impact  Test  is  also  an  exasiple. 

The  Skid  Test  and  the  Oblique  Impact  Test  hardly  qualify  as 
small-scale  laboratory  tests,  since  their  sample  Is  on  the  order  of 
20  pounds;  but  they  are  controllable  experiments,  and  they  could  yield 
fundamental  data  If  coupled  with  mathematical  modeling  such  as  is 
recommended  herein.  A  sosewhat  different  mathematical  analysis  of  the 
Pantex  Skid  Test  Is  currently  being  carried  out  at  Lawrence  Radiation 
Laboratory  (  62),  and  it  should  appear  soon. 

The  Canadian  Friction- Impact  Test,  in  which  a  torpedo  slides 
down  an  inclined  track  to  strike  a  50~mg  sample  a  glancing  blow  (  5i), 
probably  uses  too  small  a  sample  for  satisfactory  isolation  of  effects; 
but  20  pounds  is  certainly  not  necessary,  and  some  intermediate  sise 
ought  to  be  both  optlwai  and  convenient. 

Dyer  and  Taylor  have  recently  reportri  t  new  test  wherein 
a  25  Ml  cube  of  explosive  is  pressed  against  a  friction  surface  which 
Is  suddenly  jerked  away  with  a  sliding  motion  (61).  This  is  about 
the  right  slie,  although  it  may  e  just  a  bit  on  the  large  side  for 
laboratory  work,  “nils  test  ought  to  be  subjected  to  careful  a:.ialysls 
coordinated  with  a  mathematical  modeling  euch  as  is  suggested  in 
Section  S.  It  hes  great  premise. 

Dyer  and  Taylor  also  report  interesting  studies  mith  grit, 
and  their  data  aomewhat  parallel  those  from  the  Esso  Friction  Sersw. 

These  studies  also  ought  to  be  extended,  in  the  light  of  the  Esso 
results  and  in  Che  ?lght  of  the  principles  discussed  in  this  report. 


violent,  extreme  deformation  Is  employed  In  the  SUSAN  Test  (51), 
wherein  a  onc-pound  sample  of  explosive  contained  In  a  special  projectile 
is  fired  against  a  hard  target  at  controlled  velocities.  This  is,  of 
course,  an  "Impact’’  tsst;  and  LRL,  the  originator,  classifies  it  as  one; 
but  it  has  been  reported  that  the  sample  undergoes  extreme  deformation, 
expandli'^  to  several  times  its  original  diameter  with  a  corresponding 
reduction  in  thickness  before  exploding.  This  sounds  very  much  like 
plastic  or  viscous  flow  with  consequent  heating,  as  described  in  Section  3. 

For  that  matter,  any  Impact  test,  including  the  BuMlnes,  the 
NOL  and  the  Picatinny  Arsenal  Impact  tests,  look  to  this  writer  as  if 
they  had  a  large  component  of  violent  deformation  via  plastic  flow. 

This  may  be  one  reason  Impact  tests  are  so  notoriously  hard  to  analyze. 

There  are  several  other  friction  tests  currently  under 
development  by  various  groups,  but  none  of  them  show  much  promise 
of  breaking  new  ground.  The  AEG  has  one  in  which  an  abrasive  band 
pressed  against  an  explosive  sample  is  suddenly  pulled  away  with  a 
sliding  motion  while  the  pressure  is  maintained  (119).  The  "friction" 
is  greatly  confounded  with  pross  wear  of  the  explosive,  and  one  does 
not  know  whether  the  friction  is  between  the  abrasive  and  the  explosive 
or  between  -he  explosive  and  other  explosive  loaded  into  the  abrasive 
band. 


There  are  also  several  rotary  friction  machines  in  various 
stages  of  dev;-jlopme.nt.  There  was  one  at  Thlokol  a  number  of  years 
ago,  McDonnell  Douglas  Corporation  built  one  a  couple  of  years  ago 
(59),  E.R.D.E.  announced  one  in  February  1970,  (  60),  and  one  is  under 
development  at  Picatinny  Arsenal  (  120 « 

These  machines  in  general  shear  the  sample  between  rapidly 
rotating  plates  for  various  periods  of  time.  The  shear  is  something 
like  that  in  a  journal  bearing,  and  the  heat  generation  may  be  subject 
to  analysis  in  a  similar  manner.  The  treatment  is  discussed  briefly 
in  Section  3. 
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Some  suQcjested  new  approaahes 

The  planning  of  the  development  of  improved  friction 
sensitivity  testers  requires  the  identification  of  the  factors 
which  are  really  controlling  in  frictional  initiation  of  explosives, 
and  that  may  have  to  await  the  outcome  of  the  theoretical  studies 
recommended  in  Section  5.  .For  planning  purposes  here,  let  us  make 
several  alternate  assumptions  and  see  where  they  lead  us. 

Assumption  1:  The  important  fr^'otion  oharaoteristias  are 
those  of  the  solid  explosive  surface  itself.  This  assumption  is 
plausible  because  of  the  high  coefficients  of  friction  and  the  low 
thermal  conductivities  of  most  explosives  compared  with  those  of  most 
metallic  materials  of  construction.  If  this  be  the  case,  then  testers 
typified  by  the  ABL  Sliding  Friction  Machine  are  on  the  wrong  track; 
and  testers  such  as  the  Dyer  and  Taylor  machine  would  be  better.  The 
assumption  requires  that  container-on-container  friction  be  eliminated 
and  that  sliding  be  confined  to  explosive-on-explosive  or  explosive-on- 
container. 


If  the  assumption  is  verified,  then  further  development  of  the 
Dyer-Taylor  machine  is  Indicated.  It  needs  to  use  a  smaller  sample  for 
laboratory  safety  and  convenience,  and  it  certainly  could  since  most  of 
the  present  sample  is  far  removed  from  the  surface.  It  also  needs  to 
have  the  capability  of  higher  sliding  speeds  and  higher  p-essures  up  to 
the  strength  limit  of  the  explosive  block.  A  slightly  modified  apparatus 
could  also  measure  the  coefficient  of  friction  of  the  explosive  surface, 
and  that  should  be  done  by  adding  a  strain  gauge  to  the  slider. 

Assumption  2:  The  important  friction  characteristics  are 
those  of  the  container  surfaces ^  not  the  explosive  itself.  This 
assumption  is  also  plausible  because  of  the  generally  low  strength 
of  explosive  solids.  They  abrade  and  melt  under  friction,  and  the 
temperature  rise  is  strictly  limited.  If  this  be  the  case,  then  the 
ABL  machine  is  on  the  right  track  and  the  Dyer-Taylor  machine  is  wrong. 
Moreover,  the  hot-spots  must  be  generated  on  the  asperities  of  the 
test  tool  faces;  and  the  important  characteristic  of  the  explosive  is 
simply  its  thermal  sensitivity,  which  can  be  evaluated  more  directly  in 
other  ways  (for  example,  the  modified  Picatlnny  Arsenal  Autolgnitlon 
Test  developed  by  Coburn  and  Brown  (67  )). 

If  this  assumption  is  verified,  then  construction  of  an  ABL 
machine  is  recommended,  with  provision  to  interchange  tool  faces  of  any 
given  material  of  construction.  The  thickness  and  disposition  of  the 
explosive  sample  layer  probably  needs  to  be  optimized  too. 


The  Esso  Friction  Screw  provides  an  alternate,  or  perhaps  a 
complementary,  approach  to  this  assumption;  and  It  Is  a  lot  cheaper 
than  the  ABL  machine.  It  is  likely  that  the  Esso  Friction  Screw 
produces  its  hot-spots  by  grit-on-grit  friction,  with  the  explosive 
acting  as  a  heat  acceptor;  and  the  thickness  of  the  explosive  film  Is 
controlled  by  its  film  strength  and  viscosity.  The  grit  could  be 
metallic  particles,  and  the  problem  of  surface  finish  and  flatness  is 
obviated. 


Still  another  alternate  under  this  assumption  is  to  measure 
the  friction  characteristics  of  the  container  materials  without  any. 
explosive  present  and  measure  the  thermal  sensitivity  of  the  explosive 
separately.  This  requires  the  success  of  the  Section  5  theoretical 
studies,  and  it  is  probably  not  empirical  enough  to  inspire  confidence 
in  safety  decisions. 

Assumption  2:  The  important  factor  is  not  sliding  friction 
at  alii  but  rather  heating  from  viscous  or  plastic  flow  of  the  explosive. 
If  this  be  the  case,  then  none  of  the  existing  friction  testers  can  be 
very  good;  because  they  all  involve  a  great  deal  of  sliding  friction, 
and  the  flow  of  explosive  is  only  what  results  incidentally. 

A  machine  designed  to  test  this  factor  should  provide  a  fixed 
clearance  large  enough  to  preclude  container-container  contact,  and 
should  provide  a  controlled  pressure  to  cause  the  flow.  One  approach 
would  be  to  squirt  the  explosive  out  through  a  jet  like  a  hypodermic 
syringe,  with  the  test  flow  occuring  in  the  needle.  This  would  be  easy 
to  analyze  by  the  calculations  discussed  in  Section  3,  but  it  would  be 
difficult  to  construct  mechanically.  There  would  be  sliding  friction 
between  the  plunger  and  the  barrel,  and  quite  a  bit  of  plasticity  would 
be  required  of  the  explosive. 

A  better  approach  would  be  to  force  the  explosive  up  through 
an  annular  orifice  as  a  plunger  drives  down  into  a  closed  cup: 


liardened  steel 


Cast  or  packed 
explosive  sample 


The  plunger  needs  to  have  a  short  stroke  to  avoid  collision  with  the 
bottom  of  the  cup,  and  it  needs  to  move  with  a  controllable  high  speed. 
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It  also  needs  very  rugged  guides  In  order  to  maintain  close  clearance 
without  danger  of  accidental  contact  with  the  cup  wall. 

Such  a  test  is  obviously  workable  with  liquid  or  soft 
explosives;  it  is  not  so  obvious  how  well  it  would  work  with  granular 
or  brittle  explosives.  It  might  be  that  their  strengths  will  be  small 
compared  to  the  loads  imposed  by  the  test  fixture,  so  that  they  too 
will  flow.  There  might  also  be  a  large  amount  of  particle-particle 
surface  friction,  but  the  assumption  here  is  that  that  is  not  important 
anyway. 


Another  variation  of  this  concept  would  be  to  crush  a  pellet 
of  explosive  so  as  to  force  it  to  /low  out  into  a  thin  disk.  This 
might  be  done  in  a  hydraulic  or  pneumatic  press  with  a  short,  powerful 
stroke,  or  the  pellet  might  be  crushed  under  the  rim  of  a  wheel.  It 
might  also  be  done  in  an  impact  tester;  in  fact,  an  impact  tester  may 
actually  he  nothing  more  than  a  friction  tester  after  all.  The  writer 
has  yet  to  see  a  satisfying  analysis  of  its  action  on  solid  samples 
in  terms  of  compression  or  shock  waves. 

Crushing  has  certain  attractions  from  an  empirical  point  of 
view  because  it  is  a  realistic  process  and  because  it  automatically 
adjusts  itself  to  widely  different  physical  properties  of  the  explosive. 

A  liquid  or  soft  sample  will  flow,  with  heating  due  to  viscous  shear; 
a  brittle  sample  will  crush,  with  whatever  effects  the  formation  of  new 
surface  entail^  and  a  rubbery  sample  will  deform  and  recover,  with 
attendant  heating  due  to  hysteresis  losses.  The  same  considerations  make 
it  less  attractive  for  the  gathering  of  fundamental  data. 


HAZARD  REDUCTION  TECHNIQUES 


There  is  no  magic  to  eliminate  friction  or  its  effects. 
Friction  is  basic  to  the  very  nature  of  matter,  and  about  the  best  we 
can  do  is  to  moderate  its  effects  by  suitable  mechanical  design  and 
operating  limitations. 

Lubrication  is  nearly  magic,  though.  It  provides  a  low 
shear  strength  interface  and  thereby  drastically  reduces  the  coefficient 
of  friction,  and  it  also  provides  a  heat  sink  to  cool  the  hot-spots 
which  do  form.  The  best  lubricant  is  one  with  a  high  film  strength  so 
that  it  will  separate  the  solid  surfaces  and  keep  them  separated,  but 
the  materials  which  are  best  for  this  (fatty  or  petroleum  oils)  are 
difficult  to  get  rid  of  when  their  function  is  finished;  so  a  better 
choice  for  explosive  processing  is  water  or  a  volatile  hydrocarbon  or 
halocarbon,  or  a  lower  alcohol.  A  still  better  choice,  when  it  can  be 
done,  is  to  use  a  liquid  which  can  remain  in  the  finished  composition, 
such  as  the  plasticizer  in  a  composite  explosive. 


There  Is  nothing  new  about  lubrication,  of  course.  It  has 
always  been  the  practice  to  handle  primary  explosives  only  when  they 
are  wet  with  water  or  alcohol,  even  though  it  has  been  demonstrated 
that  wetting  does  not  reduce  their  sensitivity  to  Impact  or  severe 
friction  (  55).  It  does  reduce  their  sensitivity  to  gentle  friction, 
which  Is  what  one  has  In  normal  handling. 

Limitation  of  rubbing  pressures  and  velocities  is  also 
obvious  and  has  always  been  practiced.  Clearances  between  mixer 
blades  and  walls  should  be  as  wide  and  blade  speeds  should  be  as 
slow  as  Is  consistent  with  adequate  mixing.  Extruder  plungers  and 
screws  should  move  slowly.  On  the  other  hand,  bearings  and  compression 
rams  should  fit  as  closely  as  possible  In  order  to  keep  explosive 
materials  from  e  . taring  the  clearance  space. 

Grit  should  be  carefully  excluded,  because  it  can  introduce 
severe  friction  where  there  would  be  none  In  its  absence.  I*;  may  even 
be  desirable  to  substitute  softer  materials  for  Ingredients  of  the 
explosive  Itself.  For  example,  powdered  aluminum  Is  a  safer  choice 
than  powdered  boron  for  high-energy  rocket  propellants  because  of  the 
hardness  of  the  boron  (67  ).  Similarly,  the  use  of  aluminum  hydride 
gives  very  friction  sensitive  materials  because  of  the  hardness  of  the 
aluminum  hydride  (121). 

Low-fri.ction  tool  and  container  surfaces  are  desirable  on 
principle  because  of  the  danger  of  generating  hot  spots  by  their 
rubbing  on  one  another.  Brass-on-brass  is  better  than  steel-on-steel 
because  the  lower  coefficient  of  friction  and  the  higher  heat  conductivity 
both  act  to  reduce  the  flash  temperatures  attainable.  It  is  even 
better  to  make  one  surface  Nylon.  But  not  Bakelite  -  Bakellte's  very 
low  heat  conductivity  leads  to  very  high  flash  temperatures.  Surface 
finishes  should  preferrably  be  around  20  to  50  microinches,  because 
both  smoother  and  rougher  surfaces  give  higher  friction  coefficients. 

It  is  also  good  to  make  the  two  contacting  surfaces  of 
dissimilar  materials  which  have  little  or  no  mutual  solubility.  Steel- 
on-steel  can  sieze  if  the  parts  rub  hard  enough  to  cause  asperity-asperity 
welding,  but  lead-on-steel  cannot  because  lead  is  Insoluble  in  iron. 

Nylon  and  steel  are  similarly  Insoluble  and  cannot  sieze.  Most  plastics, 
however,  are  poor  choices  against  another  plastic  surface,  because  there 
is  often  mutual  solubility  with  high  coefficients  of  friction. 

Dissimilar  materials  are  good  from  another  point  of  view  as 
well.  They  give  an  opportunity  to  use  one  material  with  a  favorable 
surface  energy  and/or  yield  strength  and  another  with  a  favorable 
surface  finish  and/or  heat  conductivity,  and  have  the  best  of  both  worlds. 


Actually,  it  should  not  be  surprising  that  this  study  has 
not  produced  startling  new  insights.  After  all,  intelligent  people 
have  been  facing  these  problems  for  a  very  long  time;  and  they  have 
come  empirically  to  most  of  the  right  answers.  It  would  be  surprising 
if  this  study  had  found  them  wrong. 

What  this  study  can  do  is  to  illuminate  the  reasons  for  the 
empirical  rules  of  thumb  and  to  guide  the  choice  of  new  materials  such 
as  plastics  when  plants  are  modernized  and  the  old  rules  of  thumb  do 
not  cover  the  situation.  Plant  designers  should  get  up  to  date  on 
the  science  of  friction  by  studying  the  authorities  in  the  field,  and 
Rablnowlcz  (  1  )  is  an  excellent  place  to  start.  It  would  also  be  well 
to  retain  expert  consultants  to  review  the  friction  aspects  of  the 
design  of  new  facilities,  and  Rablnowlcz  himself  would  be  a  good  choice. 


INTENTIONALLY 


BLANK 
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5.  SUGGESTED  FRICTION  RESEARCH  TCx'ICS 


This  review  has  turned  up  a  number  of  places  where  the 
science  of  friction  and  friction  sensitivity  Is  not  yet  adequate 
but  where  there  are  clear  opportunities  for  advancement  and  where 
the  logical  next  steps  are  clear.  This  section  outlines  some  of 
the  more  attractive  and  better  defined  ones. 


A  MATHEMATICAL  MODEL  OF  FRICTION  IGNITION 


The  foregoing  discussion  has  mostly  been  rather  qualitative, 
partly  because  a  quantitative  analysis  of  the  various  processes  was 
beyond  the  objectives  of  this  study  and  partly  because  the  kind  of 
quantitative  analysis  which  Is  really  needed  has  not  yet  been  made. 
Nevertheless,  a  quantitative  study  is  needed  If  the  findings  are  to 
be  of  maximum  use.  It  is  necessary  to  determine  Just  how  large  the 
various  frictional  effects  discussed  really  are  In  situations  Involving 
explosives  and  to  compare  their  magnitude  with  those  of  empirically 
observed  explosion  hazards. 

A  quantitative  theoretical  program  is  recomnended  to  Integrate 
the  approaches  sketched  out  in  Section  3  into  a  complete,  quantitative 
model  of  the  formation  of  initiating  hot-spots.  Mathematical  expressions 
are  available  for  the  amount  of  heat  generated  In  sliding,  for  the 
number,  size  and  spatial  and  size  distribution  of  asperity  junctions, 
for  the  temperature  profiles  In  and  between  the  junctions  and  their 
variation  with  tine,  for  the  flow  of  heat  through  the  junctions  and 
Into  the  bulk  of  the  explosive  blocks,  and  for  the  thermal  evolution 
of  the  hot-spots  in  the  explosive.  This  is  essentially  a  job  for  a 
mathematician,  or  a  very  mathematically  inclined  Materials  Scientist. 

It  might  be  an  excellent  project  for  a  selected  Materials  Scientist  on 
a  post-doctoral  year. 

Host  of  the  mathematical  models  required  have  been  Introduced 
in  Section  3,  including  the  Zion  and  Mader  models  (48,  49  )  of  t'te 
growth  of  hot-spots  in  an  explosive.  The  only  additional  model  .  “oied 
is  one  for  the  critical  conditions  for  a  hot-spot  to  lead  to  explosion, 
and  this  is  furnished  by  the  work  of  Boddlngton  (SO  ).  He  has 
calculated  the  critical  sixes  and  tmeperatures  for  runaway  in  a  number 
of  coBBon  explosives  by  procx^duree  more  rigorous  than  those  of  Zinn 
and  Mader,  and  they  are  susmariied  in  the  following  figure: 
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Fiu.  3.  Critical  tvfn|icrature  aa  a  fuinlion  of  hut* 
•putaiae  for  iminc  miiimuii  rx|>luNiv.'a.  Str;ui:l>t  liiira 
ar«  tot  h  ^  0  amt  aluaiid  lie  tu  ihc  right 

liy  the  amouiita  0.301,  0.-t77  fur  It  I,  2 
tivfly.  A:  A|N»;  It:  IIMX,  I'KTN  (Krt-ab at  alo|ir) ; 
C:  UIIX;  li:  U-tryl;  K:  »‘Ui>h-m'»liiiitran>itH*; 
P:  ammonium  nilrate.  The  aht'it  curve  itim-emu 
docay.nf  n  thermal  aiiikc  rauard  hy  a  fiaami  fnm* 
moat  (oorteetml  to  k  c<iu^  I). 


Successful  execution  of  the  progrs*  envisioned  should  enable 
the  quantitative  prediction  of  the  probability  of  explosive  initiation 
by  any  given  frictional  process  without  recourse  to  eapirical  testinfi. 
knowing  only  the  fundamental  material  properties  of  the  materials 
involved  and  the  speeds,  pressures,  etc.,  of  the  mechanical  events. 
Comparison  with  empirical  reality  will  reveal  whether  the  chosen  frictional 
processes  are  indeed  the  critical  ones,  as  is  discussed  in  the  following 
paragraphs .  ^ 
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EXPERIMENTAL  TEST  OF  THE  MODEL  AND  APPLICATION  TO  HAZARD  CONTROL 


The  validity  of  the  model  developed  In  the  first  suggested 
project  will  depend  completely  upon  the  relative  importance  of  the 
processes  chosen  for  Incorporation*  and  this  furnishes  a  way  to 
determine  which  process  or  processes  actually  are  dominant  in  any 
given  hazard  situation. 

For  example,  if  a  model  assuming  the  generation  of  heat  by 
explosive-explosive  friction  is  found  not  to  correspond  to  reality; 
it  would  imply  that  that  assumption  is  not  valid  and  suggest  the 
substitution  of  a  different  assumption  such  as  the  generation  of 
heat  by  metal-metal  friction  or  by  viscous  shear.  The  model  which 
best  describes  empirical  reality  is  by  definition  the  one  which 
incorporates  the  frictional  processes  which  really  are  controlling. 

Once  the  critical  processes  are  identified,  the  lessons  for 
in-plant  hazard  reduction  will  be  obvious.  For  example,  if  the  study 
were  to  show  that  metal-metal  friction  was  not  important  and  that  only 
explosive-explosive  friction  would  lead  to  trouble,  we  could  quit 
worrying  about  materials  of  construction  and  concentrate  on  preventing 
the  sliding  of  one  block  of  explosive  on  another.  Similarly,  if  viscous 
shear  turned  out  not  to  be  important,  we  could  quit  irorrying  about 
maintaining  wide  clearances  and  concentrate  on  preventing  actual  metal- 
to<^etal  contact.  Nature  being  what  she  is,  conclusions  will  probably 
never  be  that  clear-cut:  but  there  should  still  be  lessons  for  hazard 
control. 


Modeling  will  also  reveal  tdtether  still  other  phenomena  such 
as  fracture  or  frictional  electrification  need  to  be  studied.  If  the 
models  account  adequately  for  the  empirical  observations,  "other  effects" 
can  pretty  well  be  forgotten;  however,  if  they  do  not,  it  will  be  a 
clear  indication  that  the  "other  effects"  need  to  be  explored  in  depth. 


DEVELOPMENT  OF  IMPROVED  FRICTION  SENSITIVITY  TESTERS 

As  discussed  in  Section  A,  there  is  at  present  a  wide  variety 
of  a^irical  friction  sensitivity  casters,  none  of  which  is  really  very 
much  good.  The  modelling  studies  propossd  id>ove  will  help  to  wsed  out 
some  of  them  by  rsvealing  which  ones*  prlnciplss  of  operstion  are  valid 
and  idiich  are  not.  They  will  also  help  to  guide  the  development  of  new 
and  better  ones.  This  subjsct  has  already  been  discussed  in  Section  4, 
and  it  will  only  be  recapp^  hers. 


It  is  suggested  that  three,  perhaps  four,  different  types 
of  testers  be  studied,  embodying  different  basic  assumptions  as  to 
the  critical  friction  process  or  processes.  These  would  be  (1)  A  Dyer- 
Taylor  type  machine  in  which  an  explosive  block  is  rubbed  on  a  hard 
surface  at  controlled  pressures  and  speeds,  (2)  an  ABL  machine  in  which 
the  explosive  is  a  thin  film  between  two  hard  rubbing  surfaces,  (3)  the 
Esso  Friction  Screw  in  which  the  explosive  is  a  thin  film  between  grit 
particles  rubbing  without  bulk  heating,  and  (4)  an  extrusion-type 
machine  in  which  viscous  and/or  plastic  flow  can  be  studied,  or  alternatively 
a  crushing  machine  c?  study  the  same  processes. 

This  study  is  complementary  to  the  modelling  study  irt  that 
each  machine  emphasises  some  basic  frlctlc  process  and  some  one  model 
ought  to  be  appropriate  to  each;  so  the  experimental  and  the  theoretical 
prograna  serve  as  tests  of  each  other. 

Again,  nature  being  what  she  is,  it  is  unlikely  that  results 
will  be  clear-cut  enough  to  discard  any  one  approach  completely;  and 
it  may  be  chat  all  the  testers  should  be  retained  and  used  where 
appropriate.  It  may  be  best  to  match  the  tester  to  the  particular  hazard 
problem,  with  the  guidance  of  the  theoretical  models. 


FRACTURE  AMD  ELECTROSTATICS 

This  report  is  believed  to  be  fairly  complete^inasregards 
friction  itself,  but  there  are  other  phenomena  associated  with  sliding 
and  grinding,  and  they  are  not  covered  herein.  Among  these  phenomena 
are  crushing  and  fracture,  and  frictional  electrification.  They  nay 
or  not  be  important  effects,  and  the  modelling  studies  above  are  designed 
partly  to  determine  whether  such  additional  processes  need  to  be  invoked 
to  explain  frictional  initiation  of  explosives.  If  they  are  needed, 
then  programs  arc  needed  to  investigate  them  in  depth;  and  each  program 
will  be  at  least  as  large  as  this  one. 

Fracture  and  electrostatics  are  beyond  the  scope  of  this  study, 
and  no  attempt  was  made  to  cover  the  literature  on  than.  Nevertheless, 
articles  on  them  were  inevitably  encountered;  and  the  ones  %rt)ich  were 
encountered  are  included  in  the  bibliography  to  serve  as  a  starting 
point  for  any  future  detailed  survey. 


*  Not  exhaus'.ive.  It  was  necessary  to  restrict  the  scope  of 
this  study  in  order  to  preserve  useful  depth;  but  it  is  believed  that  no 
major  Insights  into  friction  have  been  missed.  If  aiqr  have  been  missed, 
it  is  hoped  that  reviewers  will  bring  them  to  attention. 
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SUlfllARY  OUTLINE  OF  SUGGESTED  RESEARCH 

A  MATHEHATICAL  MODEL  OF  FRICTION  INITIATION 

The  amount  of  heat  available  from  sliding 
The  statistics  of  Junction  number  and  size 
The  statistics  of  Junction  temperature 
Critical  hot-spot  size  and  temperature 

EXPERIMENTAL  TEST  OF  THE  MODEL  AND  APPLICATION  TO  HAZARD  CONTROL 
Calculation  of  frictional  conditions  to  cause  ignition 
Comparison  with  oapirical  tests 

Conclusions  as  to  uhether  the  model  adequately  predicts 
events  or  idiether  additi )nal  phenomena  need  to  be 
taken  into  account. 

Conclusions  as  to  permissible  process  conditions 

DEVELOniENT  OF  IMPROVED  FRICTION  SENSITIVITY  TESTERS 

Construction  of  devices  embodying  different  basic  processes 
e  Dyer-Taylor  type  machine 
e  ABL  type  machine 
e  Esso  Friction  Screw 
e  Extrusion  or  crushing  type  machine 
Cross-checking  with  the  predictions  of  the  models  and 
comparison  %rith  known  handling  hazards 
Matching  of  the  taster  to  the  probloi 

POSSIBLE  NEW  PROJECTS 

Fracture  and  surface  energy 
Electrostatics 


INTENTIONALL 


BLANK 
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Interfacial  Temperature 
Distribution  Within  a 
Sliding  Hertzian  Contact 
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H.  A.  RAMOS 

Dt^.  «f  Ch>wical  Ewglwtring  A  OMmkai  Tachnology 
Inipafkil  Coliaga,  lemlon,  S.  W.  7 

analytic  aprasion Jar  the  tltoily-slaU  inter  facial  temperolure 
field  in  a  itidag  cmnlar  Hertzian  contact  ts  derived,  taking  into 
acemmt  the  tllifsoidul  distribution  of  the  fitictional  power  and  the 
difftreiKe  between  the  bulk  temperotiares  oj  the  two  bodies.  Jot  the 
case  where  one  act  fact  is  slatvm/uy  end  the  other  rapidh  wevii^ 
with  respect  to  the  contact  Other  caus  way  be  treated  in  a  smilar 
manner,  it  is  shown  that  the  temperature  at  cun  point  da  the 
interfioce  tan  be  approximated,  to  an  accuracy  tmprormg  with 
odocity,  by  ha^  the  harmonic  mean  of  the  two  swjact  tempeiatuits 
attained  at  that  pout  if  each  bode  were  to  receive  all  the  fiictnmal 

hightr  than  that  pten  by  Blok 's  width  used  Joemnln. 

IMTAOOUCnON 

To  char»c«ri»  compkirly  fhe  ini«rr4cc  be«v«fcen  t>»-o 
dkling  bodifL  one  mu»i  br  abia  lo  deierininv  ihc  icm- 
peranirt  disthbution  ovrr  the  ar*R  of  conuct  due  to 
frictioiul  heating  or.  tacking  tha.  the  maaimum  and  nr 
mean  tempentum.  Nf  uurement  of  il<ete  le'nperstunrt 

Mndtia  a)  Mm  YMi  AM  Aimm(  NtaaMaf 
M  CMtaoh  May  4^  tfTO 

NoMrtHXATi^u 

a  at  thermal  Aiiltttiviiv 
B  m  dimemioaltto  patametrr  vR  « 
fl  K  B  rtMnpIrtr  rtlipiK  mteara-t. 

G  at  duMrwMaalm  hulli  ttenpetaiMfr  diltrtrnrr 
I  B  thermal  twMiuiimiy 

L  M  heal  Muitte  haif-Mtdih  in  the  dtrenmn  nf  nwiKin 
f  m  Crklional  power  denui«  (heat  per  unit  time  per  tuui  atea< 
B  heat  Ilua  aerrm  the  iiurHate  due  to  the  bulk  temperature 
dilleremte 

<2  B  toul  frittmtal  power 

(h  M  power  ttawitntd  acroH  the  ituerface  aa  a  lewdi  ed  thr 
bulk  icM^Ktaiure  dwIrretuT 
e  V  dhtaace  irom  the  tenter  ol  the  heat  Miurtre 
B  m  radiua  et  thr  heat  wmne 

T  m  temperamw  at  the  iltdme  inietface  mmut  thr  Oulk  tetn- 
pent  we 


i$  a  notoriously  difficult  problem,  hence  titere  is  consider¬ 
able  demand  for  reliable  mathematical  techniques.  Sev¬ 
eral  methods  fur  calculating  interl'acial  temperatures  have 
been  devised  { /-  7i  for  the  case  in  which  the  power  density 
(i.c.,  the  rate  of  frictional  heat  production  per  unit  area) 
is  constant  over  the  contact  area.  However,  in  engineering 
systems  the  area  of  sliding  contact  is  uiten  formed  by 
elastic  deformation  of  curved  bodies  in  nominal  point  or 
line  contact,  in  which  cases  the  pressure  is  exactly  or 
approximately  cllipiicallv  distributed  over  the  contact 
area  according  to  the  classic  equations  of  Heru  Even 

sftHc  *c^d!  vh*  iNrm*«ikr# 

distribution  is  still  fairly  close  to  elliptical  lit).  I'hus,  given 
that  the  frklion  force  contnbuted  by  a  differential  ele¬ 
ment  of  area  is  proportional  to  the  normal  load  on  that 
element,  the  power  densitv  m  a  sliding  Heruian  contact 
will  be  ellipiieallv  distributed.  .\n  important  example, 
and  the  prinetpai  application  fur  the  trvulis  of  this  palter, 
is  the  4-ball  machine,  tor  a  pin-on-disc  machine  with  a 
hemisjdmical  slider),  widely  used  in  inbologv-  meareb 
and  lubricant  testing,  which  effei  is  sliding  contact  be- 
iween  two  tphem  such  that  tlut  contact  area  it  staiionarv 
vm  one  ol  them. 

r,  a  surface  tempcratuie  wnhin  a  uatMt-vars  hrai  source  on 
a  single  surCace  minus  the  bulk  irmperaiute 
r,  «  surface  temperature  wnhin  a  moving  heat  wHirce  on  a 
smgle  surface  minus  thr  bulk  lemperaiute 
•  ■  dimensionlevs  t«ordnuie  ’  a  I  a-  s  /.i 
e  ■  sekoetty  with  n-vpect  lo  the  heat  vniice 
a,  ,r  M  Cartesian  -  •inrdniaies  denotint  pnviiHtn  wnhin  the  heat 
wuree 

B  n  «  cstndinaies  of  the  maximum  inirrfaisal  tempera- 
twe 

a  a  true  heal  partition  cocfhcimi 
B  at  eAceiise  hrai  partninn  corfheimt 
hm  M  local  desiaiion  ftum  a 
t  ■  absolute  imiprratuie  at  the  siohne  inietface 

*  hulk  trmpr<aiuirs  nf  the  Maininarv  and  moving 
bndict.  frsprriivels 
m  bulk  irmprtalurc  dilferenre 
It  B  dunenMunli  m  lunriion  of  a 
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lii  ilii'  |>,i|>i'r,  .III  iiii.ilvtii  <•x[)I^s^illIl  for  the  sieudy-slatc 
II. o  il. II  i.ii  i(  iii|)riatiiri-  (lolriliiiiioii  in  a  radially  sym- 
I'li  iiii  I  Irri/Mr.  l  oni.ii  t  '.Iiiliii”  iiiulcr  unluliricatrd'or 
>1. .mill. ir\  lulu u. Ill'll  (oi'.diiims  is  (li-rivi'il  takin^  into 
.uiiMiiil  till'  rllipsoid.il  pinvi-r  ili'iisity  distriliiiiiiin.  The 
iir.iiMii'iit  is  |iii'sriiii'd  for  the  rasr  in  svliirh  one  surface 
i>  Ni.iiion.irv  wiili  ri-spi'ci  to  ilii-  coniart  area  hut  is  valid 
.o  well  I'll  lerlain  oilier  eases  svliieh  will  be  discussed. 
We  ssill  Ills!  derive  the  lemperaluie  disiribution  on  a 
..•v.Y  siirl.iie  for  (a)  a  stationary  beat  source  and  (b)  a 
iiiosiiu:  lie.tt  source.  The  teinper.iluie  distribution  for 
'luliiu;  lonlaci  will  ilieii  be  obtained  front  these  two 
(lismbiitioiis,  iniliallv  In.  the  case  where  the  two  bodies 
have  the  same  bulk  lenipeiatuie  and  subser|ueMlly  for  the 
more  general  case  of  une.|ual  bulk  temperatures. 

In  this  pa|)er,  the  svmbol  7' will  always  ri'iiresent  the 
surface  temperature  imiruir  above  the  temperature  in  the 
bodvisi  far  Iront  the  heal  souice.  I'lie  subscripts  .t  and 
I  '.sill  refer  to  the  smitle  surface  stationary  and  moving 
sounes,  respectively. 


fig.  1— |o)  Uniform  anil  (bj  nlliptoidol  powor  danti*y  ditIribuHant. 
Beth  diilribuliont  hnva  th*  loma  total  powor  0. 


the  elastic  deflections  tf(r')  for  the  corresponding  pressure 
distributions  (//).  For  the  ellipsoidal  power  density  (3), 


(4), 


STATIONARY  HEAT  SOURCE  ON  A 
SINGLE  SURFACE 


The  steady  stiiie  (i.e.,  after  an  infinitely  long  time) 
stirlace  temperature  due  to  a  stationary  point  heat  source 
on  a  plane  surface  is  given  by  ( 10) 


dr  ~ 


•hibk  ’ 


m 


where  b  is  the  distance  from  the  heat  source,  i/Q  is  its 
power,  and  k  is  the  thermal  conductivity.  An  identical 
relationship  holds  lor  the  elastic  normal  deflection  of  a 
pi. me  surface  tine  to  a  norntiil  point  force  dl\'  acting  on 
the  surface  (/  / ); 


dw  = 


dW 
2i:bE’ ' 
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where  E'  =  E/2[\  —  r'-’),  E  is  Young's  modulus,  and  »•  is 
Poisson's  ratio.  Therefore  the  normal  deflection  u<(x, 
of  a  plane  surface  due  to  a  pressure  distribution  p(x,  y), 
obtained  by  double  integration  of  Eq.  |'J|.  is  an  exact 
analogue  of  the  steady-state  surface  temperature  field 
7‘,(  V,  rl  resulting  from  a  stationary  heat  source  of  power 
dei.sity  i  )  distributed  over  a  finite  area  on  the  surface, 
which  is  given  by  double  integration  of  Eq.  (1). 

An  analytic  expression  is  required  for  T,(t)  resulting 
from  the  ellipsoidal  power  density  distribution 


-(r/W^  [3] 


where  O  is  the  total  power  produced  by  the  heat  source, 
R  is  its  radius,  and  r  is  the  distance  from  its  center.  This 
is  plotted  in  Fig.  1  in  the  dimensionless  form  (2ir7fV3QJ 
For  comparison,  the  uniform  power  density  dis¬ 
tribution  q  =  QJ'tR-  is  plotted  as  well.  Using  the  analogy 
presented  above,  the  surface  temperature  fields  r,(r)  for 
these  two  heal  distributions  can  be  written  directly  from 


and,  for  the  uniform  power  density, 

T,(r)  =  i£(, //?)-£-  (r  <  R),  [5] 

ir  irRk 

where  E(t/R)  is  the  complete  elliptic  integral  of  the 
second  kind  whose  modulus  is  r/R.  In  Fig.  2,  the  two 
temperature  distributions  are  plotted  in  the  dimensionless 
form  (tr/fl/Q)7',(r//f).  The  average  temperature  over  the 
source  area  (r  <  /?)  for  the  ellipsoidal  distribution  is 

Z2-nRk' 

which  is  4.1c;.  greater  than  that  for  the  commonly  em¬ 
ployed  uniform  heat  source. 

FAST  MOVING  HEAT  SOURCE  ON  A 
SINGLE  SURFACE 

Let  (x,  y)  be  a  Cartesian  coordinate  system  in  the 
surface  with  its  origin  at  the  center  of  the  heat  source. 


Fig.  3— Surfoc#  fo.nparolura  field  within  •  (tatfonary  hool  muk*  for 
(o)  uniform  on  f  |b)  •tliptoidol  powor  doiivity  ditlHhuHofU. 


SO 


+  t  the  .*»!h  »riocrf»  s  o# 

surface  »»i!h  rrtprt!  i<>  ibr  wHircr  Aftharri  (?)  has 
vh(n»n  ihai  when  ihe  dimen\H>nle«  pafametrr  r#-«  > 
!().  where  a  is  ihe  ihermai  dilfusiviiv.  ihe  heat  flow  in 
the  y  direction  may  be  neclecied.  and  the  temperature 
distribution  lin  a  heat  source  of  finite  area  can  be 
determi  ed  by  dividing  the  area  into  diflerential  stript 
parallel  to  the  sliding  direction;  the  temperature  profile 
along  any  strip  is  the  same  as  that  for  an  infinitely  long 
(in  they  direction)  band  source  (for  which  naturally  there 
is  no  heat  flow  in  the>  direction)  of  width  equal  to  the 
strip  iength  and  having  the  same  power  density  profile 
as  the  strip.  The  situation  is  illustrated  in  Fig.  3  which 
shows  the  ellipsoidal  power  density  distribution,  a  differ¬ 
ential  strip,  and  its  equivalent  semi-infinite  band  source. 

To  apply  this  method  to  our  fast  moving  ellipsoidal 
heat  source,  we  require  an  analytic  expression  for  e 
temperature  profile  across  a  band  source  having  an  ellip¬ 
tical  power  density  profile  (a  vertical  section  through  an 
ellipsoid  being  an  ellipse).  An  exact  analytic  solution  is 
not  possible  (2),  however  Blok  (/2)  has  derived  an  ap¬ 
proximate  expression  whose  accuracy  improves  asymptot¬ 
ically  with  increasing  vL/a,  where  the  band  width  is  2L. 
Blok’s  relationship  can  be  expressed  most  conveniently, 
after  a  change  of  variable  and  rearrangement,  as 

I’l 

where  u  =  +  ^/^)'  ^  replaced  by 

i[2M(l  -  r)  -  1),  Since  q(x)  =  -  {x/Lfy'^,  tq. 

(7)  is  an  elliptic  integral  which  reduces  to  {13) 

where  il{u)  =  (1  —  u)K{u^^^)  +  [2u  —  l)£(u*^*),  and 
K{u^^'^),  E{u^^^)  are  the  complete  elliptic  integrals  of  the 
first  and  second  kind,  respectively,  both  of  whose  moduli 
have  the  value  u‘^".  The  function  ^(u)  is  plotted  in  Fig. 
4. 

Cameron,  et  al.  (6)  showed  that  for  a  uniform  {q{x)  =  e 
constant)  band  source,  the  temperature  profile  r,(u)  is 


Fig.  4— Th*  diminuenltii  function  f2(u)  which  pretcriboi  Iht  shop* 
of  th#  lurfoco  timporoluro  prefilo  ocrou  on  infinitoly  fast  moving  bond 
MMirco  b';  ving  on  ollipticol  power  density  profile.  The  leading  edge  of 
Ihe  bond  is  on  Ihe  left. 


given  almost  exactly  for  10  <  nL/a  <  x  by  multiplying 
the  asymptotic  (I'L/a  =  oo)  expression  by  the  factor 
1  +  u"*^-j0.65(a/t'Z,)  -I-  0.44(a/t’Z,)^^“|,  Assuming  this 
correction  factor  to  be  valid  for  the  elliptical  case  as  well, 
Eq.  [8]  becomes 


TJu)  = 


3k 


m^r 


As  vL^S  d'?CrC2S'?S  th^  ror!*?r**0^  r^r'lor  in- 

creases  the  maximum  temperature  and  moves  it  toward 
the  origin. 

Following  Al  chard’s  method,  we  now  can  write  directly 
the  temperature  distribution  within  a  fast  moving  circular 
heat  source  having  an  ellipsoidal  power  density  distribu¬ 
tion  by  setting  L{j)  =  {R^  —  and  =  {3Q/27r 

/?2)(1  -  ~  13QZ,0’)/2ff£3)(secFrg.3).Thu.'=we 

have,  from  Eq.  [9] 


+  0.44(j^r*])L>/W(u).  [10) 


This  expression  holds  for  B  =  vR/a  >  ~10.  For  B  10, 
vL/a  for  the  outer  strips  will  be  less  than  10,  however  this 
would  not  be  expected  to  contribute  a  serious  error  be¬ 
cause  (a)  the  correction  factor  will  hold  approximately 
for  vL/a  <  10,  and  (b)  the  outer  strips  are  only  a  small 
fraction  of  the  area  of  the  source;  e.g.,  L  <  R/2  for  only 
6%  of  the  toul  area. 

The  maximum  temperature  and  its  location  (x,  j>)  are 

r  .  d  at  (0.652£,  0) 

^  ttRMB^^^  for  the  asymptotic  case, 


fig.  3— ImomMc  vIow  showing  Iho  fast  moving  oiligtoltlal  gowor  q  .  (0  633/2.  0) 

doniity  disihbutlon,  o  dillotonllal  strip  (wWHi  dy,  longih  21|i  and  its  T  =  2.013  ^ ' 

ogulvolont  somi^nAnilo  bond  sourc.  (width  ai|.  nRAB^'^  for  B  =  10. 
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For  =  10, 


7: 


1.148 


(14J 


TTius,  over  the  range  of  applicability  of  Eq  (10],  the  effect 
of  the  correction  factor  is  chiefly  amplification  of  the 
temperature;  th<-  translation  of  the  maximum  is  negligi¬ 
ble.  For  the  asymptotic  case,  the  maximum  and  ayerage 
temperatures  are  16. 1'i  and  4.3  <  higher,  respectively, 
than  those  (3)  for  a  circular  unifiem  heat  source,  where 
the  maximum  temperature  occurs  at  the  trailing  edge 

{R,  0). 


SLIDING  HERTZI/>N  CONTACT 


wjihio  the  CDRUci  area  she  irmpr.atiirr  fs.'lds  of  'he 
ftnomg  and  siaticnan  turfates  at  the  inifTfate  »»ju*d  be 
«  7,**,  rt  and  =  I  —  «» i7,i  i.  » ■.  respect  is  civ.  vhere  and 
r,  are  as  ilclined  b\  Eqs  jldl  and  [}j  However,  ii  was 
pointed  out  bv  Blok  ( I !  that  if  the  trinperasurc  fields  on 
the  two  surfaces  are  to  coincide  over  the  entire  contact, 
the  ratio  of  the  two  heat  fluxes  into  the  surfaces  cannot 
be  constant  over  the  contact:  assuming  that  the  interfa¬ 
cial  temperature  field  lies  between  oT',,(t,  r)  and  (1  —  a) 
T,(x,y),  then,  in  the  regions  where  n  7", ,  <  (1  —  a)T,,  the 
fraction  of  the  local  frictional  power  entering  the  mov¬ 
ing  body  will  have  values  greater  than  o,  and  conversely 
for  the  rest  of  the  contact  area.  Thus  the  distributions  over 
the  contact  of  oQand  (1  —  a)Q_  will  in  fact  be  (a  +  Aa) 
^(r)  and  [1  —  (o  Ao)|^(r),  where  Aa  is  some  unknown 
function  of  x,  j>  which  expresses  the  deviation  from  the 
ellipsoidal  heat  flux  distributions  ag(r)  and  (1  —  a)q(r). 
However,  a  reasonable  and  mathematically  practical 
approximation  to  the  resulting  interfacial  temperature 
field  Tlx,  y)  may  be  had  by  assuming  that,  since  at  a  given 
point  aT^,  and  (1  —  o)7',  are  linearly  dependent  on  oQ 
and  ( 1  —  n)(l,  the  deviations  in  the  two  surface  tempera¬ 
tures  caused  by  the  heat  flux  deviations  ^Aa(x,^)^(r)  are 
similarly  proportional  to  Aa(x,  y)Q.  Then 


Now  that  we  have  both  the  stationary  and  high  velocity 
surface  temperature  fields  within  an  ellipsoidal  heat 
source  on  a  single  surface,  the  task  remains  to  determine 
the  steady-state  interfacial  temperature  field  for  sliding 
Hertzian  contact  W'e  assume  here  that  ihe  isso  liodies 
have  the  same  bulk  temperature  (i.c.,  the  temperature  far 
from  the  contact);  the  more  general  case  will  be  treated 
in  the  next  section. 

Assuming  that  a  ithin  the  contact  area  the  gap  between 
the  surfaces  is  sufficiently  small,  which  obtains  for  “dry” 
sliding  or  boundary’  lubrication,'  there  will  be  no  tem¬ 
perature  discontinuity  across  the  interface.  The  moving* 
body  will  receive  a  greater  fraction  n  of  the  total  frictional 
power"  than  the  stationary  body,  becau.se  the  stationary 
side  of  the  interface  needs  only  enough  power  (1  -  o)Q 
to  mamtatn  the  steadv-state  temperature  field,  while  the 
moving  .side  of  the  interface  requires  the  greater  amount 
of  power  nQ  to  raise  its  temperature  from  the  ambient 
value  to  the  higher  values  at  the  interface.  As  c  increases, 
a  proportionately  greater  area  of  the  moving  surface  must 
he  “brought  up  to  temperature”  per  unit  time,  and  thus 
the  heat  panition  coelTiciciil  a  increases  to  accommodate 
the  increased  power  requirement;  for  the  asymptotic  case, 
a  =  1. 

Intarfacial  Tamptrertur* 

If  the  frictional  power  density  q(r)  were  to  divide  itself 
between  the  two  bodies  in  the  ratio  a/(  1  —  a)  everywhere 


'Thicker  elastohydrodynamic  lubricant  films  preteni  a  more  complex 
problem  (14,  15' 

n'he  terms  "moving"  and  "stationary"  are  defined  with  respect  lo  the 
'  ont  ict  area. 

'Given  here  by  the  product  of  velocity,  normal  load,  and  coefficieni 
of  Inction. 


(o  y-  Aa(x,>)]r^,(x,^)  =  T\x,y) 

=  (1  -  a  -  Aa(x,  v)]7',(x.^),  [15] 

and 

T{x,y)= - j - ? - j -  (5  >-10).  [16] 

T,{x,  y)  r^(x,  y) 

Thus  the  interfacial  temperature  field  in  the  circle  of 
contact  can  be  approximated  by  half  the  harmonic  mean 
of  the  two  “single  surface”  temperature  fields— that  is, 
the  temperature  fields  which  would  exist  if  each  body 
were  to  receive  all  of  the  frictional  power  Q.  Note  that, 
due  to  the  assumption  leading  to  Eq.  [15],  T\x,  y)  is 
independent  of  a. 

The  subsequent  presentation  will  be  restricted  for  sim¬ 
plicity  to  the  case  where  the  two  contacting  bodies  have 
the  same  thermal  constants  k,  a  but  can  easily  be  followed 
through  for  dissimilar  materials;  the  same  arguments 
hold. 

7',.(y,  y)  is  a  decreasing  function  of  B,  while  T,(x,  _y)  is 
independent  of  B.  Thus  from  Eq.  [15],  asB~*  ao,  IkaTJ 
AftT, -eO,  and  7\x,  )r)-»a7’„(x,  >),  while  aT^(x,  y)-* 
T^{x,y).  For  the  asymptotic  case,  1\x,  y)  -  T^{x,y).  This 
is  apparent  in  Fig.  5  which  shows  the  profiles  along  the 
X  axis  to  T,  oT,,,  (1  -  a)T,  plotted  in  the  dimension¬ 
less  form  (7rRk/Q)T{x,  0)  for  B  =  10,100.  The  value  of 
o  was  calculated  using  Eq.  |23]  derived  below.  Thus  for 
high  velocities  the  temperature  field  of  the  stationary 
surface  is  forced  to  cunf^orm  to  a7',,(x,  >),  and  the  form 
of  (1  —  a)T,[x,y)  is  of  little  consequence.  Stated  another 
way,  the  power  which  flows  across  the  interface  in  order 
to  equalize  the  two  temperature  fields  is,  for  large  values 
of  B,  a  negligible  fraction  of  the  power  oQ  received  by 
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Figt.  So,  Sb— Surfac*  tampcrahir*  preflUt  along  th«  nntorlina  «l  a 
ilMing  Hortiion  contact  for  two  different  eliding  velecitiei.  The  leading 
edge  ef  the  centect  circle  i(  on  the  left.  T  it  the  interfacial  temperature, 
T,.  the  temperature  ef  the  moving  body  if  it  were  to  receive  all  the 
frictional  heat,  and  aT^,  1 1  -  o)T,  are  the  two  tingle  turfoce  tempera- 
turet  taking  into  account  the  partition  ef  the  hoot.  At  ■  i'.creatoi, 
r-*aT,. 


the  moving  body,  and  thus  7Tx,  >)  does  not  significantly 
depart  from  aT^(x,  y).  Consequently  the  accuracy  of  the 
approximation  [16]  increases  with  B. 

Since  T",  and  both  lie  on  the  x  axis,  it  can 

be  seen  from  Eq.  [16]  that  the  maximum  interfacial 
temperature  T„„  will  occur  on  the  x  axis  as  well.  Putting 
^  =  0,  L  =  /?,  r  =  |xl  =  /?|2u  -  1 1,  Eqs.  [4]  [10]  [16] 
give 


irM 


nx.o) 


a/3ir 


1  -  %2u  -  1)» 


Q(ii) 


1  4- 


0.65  .  0.44 

B 

„l/2 


SrnHWf'^ .  ‘7,4 1.^.  0*  varxs  bv  and  can  »r  appro* 

inuiied  by 

0)  0.540  OSSOB  '  ”*  [19] 

Finally,  from  Eqs,  [16]  [18]  [19], 


7’»..=  7Tx.,  0)  = 

1.852(?/tr/tt 

1.996  -  l.OOlB-o"*'*^  o!537fl-“*”  + 

(B  >  10),  [20] 

which  provides  a  working  formula  for  computing 

Haol  Pgrtition  Cogfflcitnt 

Over  the  clement  dxdy,  the  power  entering  the  moving 
body  is  {T/T^)q  dx  dy  =.  qdx  dy/[\  +  (r^/r,)].  Thus  the 
heat  partition  coefficient  a  is  given  by 


[21] 


However,  since  this  expression  cannot  be  integrated 
analytically,  and,  besides,  it  is  only  as  good  an  approxi¬ 
mation  as  Eq.  [16],  we  offer  the  following  alternate 
method  for  determining  a. 

The  best  approximation  to  the  true  value  of  a  will  be 


perature  fields  of  the  two  surfaces  be  identical  over  the 
contact  area.  Thus  an  analytic  approximation  can  be 
obtained  by  equating  the  mean  temperatures  of  the  two 
surfaces,  as  did  Jaeger  (2),  giving 


=  (I  -  a)r,,. 


Substituting  Eqs.  [6]  [13]  yields 


[22] 


0.870B>^^ 

1.047  .  0.774 
B  '*"B^ 


(B  >  10)  [23] 


which  expresses  the  manner  in  which  a  approaches  unity 
as  B— f  00.  At  the  lower  limit  B  =  10,  a  =  0.709;  the 
moving  body  receives  roughly  twice  as  much  power  as 
the  stationary  body. 


[17] 

which  is  plotted  in  Fig.  5  for  B  =  10,100.  The  position 
at  which  occurs,  obtainable  from  Eq.  [17],  de¬ 
pends  only  on  B  and  is  equal  to  0.652B,  0.396B  for 
B  a  00, 10,  respectively.  Over  the  range  oo  >  B  >  10, 
varies  by  14%  and  can  be  approximated  by  the 
expression 

Ci/irB*r.(x„,  0)  =  1.078  -  0.290B“»”‘.  [18] 


UNEQUAL  BULK  TEMPBULTURB 

In  the  preceding  section,  the  bulk  temperatures  of  the 
two  sliding  bodies  were  taken  to  be  equal.  In  practice, 
however,  this  condition  seldom  occurs.  If  the  heat  parti¬ 
tion  coefficient  a  is  greater  than  ‘/j,  and/or  if  the  thermal 
masses  or  heat  loss  characteristics  of  the  two  bodies  are 
unequal,  the  difference  between  the  bulk  temperature  will 
increase  asymptotically  with  time  due  to  the  frictional 
heating.  In  addition,  of  course,  external  heating  can 
contribute  to  the  bulk  temperature  difference. 


nem  omwii>irr  thr  rtirci  o#  tirnn|tul  !)>«%  tf-mfirta- 
lurrt  ‘>n  the  Hfrfjctji  frmprraiiurf  held  In  addiium  fo 
iha*  fntiionai  pn»rr  ti  lo  both  btaiir^  fnHH  ih«‘ 

inirrfacr.  ihetr  «  now  <ui  amount  of  poM^rr  which  flow 

from  the  hotter  body  airtru  the  iiurrlace  into  the  tuider 
body.  The  siluaifion  ly  illutiratcd  in  Fii?  b  vktiere,  remem- 
berini;  that  the  symbol  7'  refers  to  the  temperature  in¬ 
crease  above  the  bulk  lemperaturr  and  desii(natini{  abso¬ 
lute  temperature  bv  the  symbol  6*.  fh  e,  y)  is  the  absolute 
interfacial  temperature  Held,  and  ^5,  are  the  two  bulk 
temperatures  (assumed  to  be  uniform  far  from  the  inter¬ 
face).  The  form  of  the  heat  flux  distribution  .t,  y)  is 
indeterminate.  However,  if  we  assume  (/jf.r,  y)  to  be  ellip¬ 
soidal,  then  (l^  and  (^become  additive,  greatly  simplifying 
the  mathematics;  effectively  0^  only  alters  the  heal  parti¬ 
tion  coefficient  a.  .As  long  as  is  sufficiently  less  than 
the  error  due  to  assuming  an  arbitrary  specific 
(x,y)  distribution  will  be  small.  Taking  Q_^  to  be  positive 
when  it  flows  from  the  moving  body  to  the  stationary 
body,  we  define  a  parameter  ft  such  that 

(i»  =  (a  -  m  [241 

Thus,  when  8^  -  6^,  =  Atfj  =  0,  Qj  =  0  and  ft  =  a. 
and  hence  a  —  fi  have  the  same  sign  as  Affj,.  The  net  heat 
flow  from  the  interface  into  the  moving  body  is 

=  [251 

and  the  net  heat  flow  into  the  stationary  body  is 

(I -«)Q+Q»  =  (l -^)Q.  (26) 

Thus  P  is  the  effective  heat  partition  coefficient.  Note 
that,  w  hercasn  ranges  only  from  */••  (when  B  =  0)  to  unity 
(when  B  =  00),  P  can  have  any  value,  positive  or  nega¬ 
tive.  depending  on 

Ihe  absolute  interfacial  temperature  field  6{x,  y)  can 
be  derived  as  in  toe  preceding  section;  the  same  argu¬ 
ments  apply.  Eq.  (15)  becomes 

\P  +  Aa(*.  jr)ir„(x.  y)  +  =  0{x,  y) 

=  |l  -  Ao(.r.  y)17;(r..y).).(lj,.  [271 

Note  that  the  single  surface  temperature  fields  PT,.{x,  y) 
or  (I  —  P)7\(x,  y)  will  be  negative  (cooler  than  the  bulk 
temperature)  when  heat  flows  out  of  the  respective  body — 


«a  Q, 


(i-adO  0» 

ey. 


Sthtnwiic  raprtMntsKon  of  tho  hoot  Aotw  at  llio  lIMinp  intor- 
fxo.  Q  it  iho  liicliona!  power,  and  O,  it  tho  hoot  flow  rotullinp  from 
Iho  bulk  lompotaluro  difloronco  .Ut,  -  (*»,. 
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eqtmakrnt  ic»  ait  eflppiaedkdl  !<«■«  «Mt  .m  tW  aaeists. 

I'fom  Eq  {27S, 

8<X.  yt  =  9^  +  — - U^lIL-.  (B  >  -10).  f28) 

T,{x,  y)  ^  r,(*,  y) 

As  before,  9{x,  y)  is  independent  of  p. 

AsB-r  ».#(*.  y)  -  ~  {•^»/7',(*.J'))17',.(.v,>). 

Thus  when  the  bulk  temperatures  are  unequal,  the  inter- 
facial  temperature  held  for  the  asymptotic  case  no  longer 
conforms  to  7’^(.t,  y).  In  addition,  T^{.x,  y)—*0,  hence 
8(x,  y)-*  9 

In  Fig.  7  the  relative  values  of  9^^,  9^,  arc  plotted 
as  functions  of  the  dimensionless  bulk  temperature  differ¬ 
ence  G  =  xrRklOf^ '()  over  the  range  —  1.3  <C  C  <[  1.3  for 
B  =  10,100.  The  horizontal  line  is  the  average  bulk  tem¬ 
perature  y.,(e^  +  ej.  As  B-»  oo,  and  this 

trend  may  be  seen  in  the  plot.  The  location  of  the 
maximum  interfacial  temperature  on  the  .v  axis  varies 
with  G as  shown  in  Fig.  8.  When  G  =  -f  1 . 1 78, 9^^  =  9/^ 
for  all  values  of  B,  and  9^^  is  situated  at  the  center  of 
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Fig.  7— Rolotivo  volwoi  of  tho  maximum  intorfociol  tomporohiro  ond 
Iho  two  bwlk  tomporolurot  ploHod  01  funclioni  of  Iho  dimoniionloit 
bulk  tomporoluro  difforonco  C  for  hwo  difforont  iliding  volocilioi. 


Fig.  I— location  of  Hio  moximum  inlorfocial  fomporohiro  on  Iho  con- 
toriino  m  a  function  of  Iho  bulk  fomporoturo  difloronco.  Tho  ordmoto 
'  .0  li  Iho  trailing  odgo,  and  0  it  Iho  conlor. 


-f.,  !  _ 


tfar  amtaet  t»'  AaG  t _ Rums  t»«rani  V 

o:  the  cDnu(.< 0). 

Tbe  efecmc  heat  partition  coeWkicBt  fi  »  given  by 


4f'- 


However,  as  in  the  preceding  section,  we  can  obtain  an 
anaiytic  expression  for  fi  by  equating  the  average  tem¬ 
peratures  of  the  two  surfaces  over  the  contact  area.  Thus, 


+  ».,  =  (  1  -  +  9u  1301 


iCi,  f)  =  '  v  .  ;36i 

Eq.  (35 1  is  also  applicable  to  the  one-dimemional  pr>birm 
of  a  sliding/ rolling  Hertzian  contact  band  such  as  would 
result  from  contact  between  gear  teeth  or  cylinden,  here, 
T^j(s)  are  of  the  form  of  Eq.  (9). 

The  maximum  intcrfacial  temperature  is  usually 
the  most  useful  single  quantity  and  the  one  most  often 
referred  to  as  the  flash  temperature.  In  his  landmark 
paper  on  flash  temperature.  Blok  ( /)  derived,  for  a  circu¬ 
lar  contact  area  with  a  uniform  power  density  distribu¬ 
tion,  an  expression  for  based  on  the  two  assumptions 


which,  together  with  Eqs.  [22]  (6],  gives 


aT,{R,  0)  =  (I  -  a)r,(0,  0) 


P 

i-P 


1  -  1.132C 

1.132C  ’ 
a 


(31)  and 


=  y2(«7'„(0,  0)  -(-  ( 1  -  a)7’.(0,  0)].  (38] 


where  (1  —  a)/o  is  as  given  by  Eq.  [23].  Rearrangi  g 
Eqs.  [24]  and  [31], 

Q,/Q  =  a -/3=  1.132aC,  [32] 

ft  =  a-  1.132aC;  [33] 

✓ 

thus  both  P  and  the  ratio  Qj/Q  are  linear  in  G,  ihe 
constants  depending  only  on  B. 

DISCUSSION 

A  number  of  cases  where  both  surfaces  move  with  re¬ 
spect  to  the  contact  can  be  treated  by  the  method  of  this 
paper.  Jaeger  (2)  showed  that  for  low  velocities  (B  < 
—0.2)  the  moving  single  surface  temperature  distribution 
may  be  approximated  by  that  for  a  stationary  source. 
Thus,  designating  the  velocities  of  the  two  bodies  with 
respect  to  the  contact  area  by  r,,  Uj,  we  can  determine 
9(x,y)  for  any  combination  of  values  of  B,.  Bj  greater 
than  + 10,  less  than  - 10,  and  from  —0.2  to  -f-0.2.  When 
Bj  >  10,  B,  can  range  from  -0.2  to  -1-0.2  without  in¬ 
validating  the  results  of  the  previous  sections.  When  both 
]Bj],  ]B2|  <0.2,  r,  is  replaced  by  7",,  and  Eq.  [28] 
becomes 


9{x,  y)  =  %{9^  +  9J  +  %r.(x,  y).  [34] 

When  both  ]B|1,  [B,]  >  10,  9(x,y)  will  be  given  by 

1 _ 

9(>t,y)^9^  4-  — 


y) 


1 


[35] 


^*t(A  y)  y) 


If  V|,  Vt  are  of  opposite  sign,  T^^ix,  y)  must  be  replaced 
by  r,2(— Jt,  y).  For  the  case  of  pure  rolling  contact^ 
(»,  =  Bj),  Eq.  [35]  becomes 


Although  the  heat  partition  criterion  [37]  is  convenient 
in  that  it  docs  not  require  knowledge  of  the  complete 
functions  r,(*,  y)  and  TJ,x,  y).  it  seems  a  rather  poor 
approximation  to  the  condition  that  the  temperature 
distributions  of  the  two  surfaces  must  coincide  over  the 
entire  contact  area.  The  assumption  [38]  neglects  the  fact 
that  at  high  velocities  T  approaches  7',..  as  shown  earlier, 
and,  in  addition,  it  is  likely  that  occurs  at  a  point 
nearer  to  the  trailing  edge  (R,  0)  than  to  the  center  (0. 0). 
However,  as  Blok’s  equation  is  the  one  commonly  used 
to  calculate  flash  temperatures,  it  is  the  logical  equation 
with  whicli  to  UUI  laCJ. 

The  equation  derived  by  Blok  (/)  for  B  >  —20  is 


l,362Q,/ffB* 
1.596  4-  B'''^  ■ 


(39) 


The  ratio  of  Eq.  [20]  to  Eq.  (39j  ts  l.iJ  for  B  =  100  and 
1.38  for  B  =  10,  the  lower  limit  of  applicabilitv  of 
Eq.  [20], 

Experimental  support  for  the  validity  of  Eq.  |20]  is 
contained  in  data  presented  by  Fein  (/7,  18k  who  found 
a  discrepancy  between  lubricant  transition  temperatures 
in  a  slow  speed  pin-on-disc  machine  and  those  determined 
by  calculating,  using  Blok’s  equation,  at  the  transi¬ 
tion  point  in  a  4-ball  machine  where  the  parameter  B 
was  in  the  vicinity  of  10.  He  further  showed  that  this 
discrepancy  could  be  eliminated  hv  multiplying  the  nu¬ 
merator  of  Eq.  ]3yi  by  1.33  which  comp.ires  favorably 
with  the  ratio  1.38  between  Ecjs.  |20).  |39|  at  B  =  10. 

It  should  be  noted  ii,  passing  that  a  lomntonly  over¬ 
looked  facUK'  in  the  calculation  of  flash  tempt  ature  is 
the  fact  that  the  value  of  the  thermal  ettnductivity  k  (to 
which  the  thermal  dilfusivity  ti  is  proportional)  can  vary 
markedly  with  composition,  thermal  history,  and  tem¬ 
perature  ’•'or  example,  the  v.ilues  of  X  for  various  tvpical 
iron  alloys  listed  in  Table  I  ( 19)  are  seen  tt»  var\-  through 


*lt  ihould  be  noted  that,  for  pure  rolling  in  the  absence  ul  a  hydro¬ 
dynamic  lubricant  him,  the  rnciional  power  cunti^ls  priiiiarilv  ut  the 
hyttereiii  in  the  elastic  defurmation  cycle  and  is  diMohuted  ntiiumi- 


formly  tiver  a  ro/iiw  Mhiwe  rlTnovr  radiiit  Irnm  the  rentri  ul  luni.u  i 
w  iM-rh.'in  J  or  '  liniri  H  >  l'<i  'I  hr  rlhpMiid.il  pl.iii.ir  hc.ii  o.i„,  t  ,„,n 
iherclorr  be  a  los  ai.urair  approniinjiiun  when  ihere  it  no  shihiig. 
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r.«i  t  IfILRMAL  rivrn  ttf  VaR!»»1>  iron  4LLr»VS 


Ailov  Tvri 

C«»ia«»*irnw  '.weight  ■:  i 

Tmaaaat 

CositH  c.ttvmr 
(calyctn  *ec  ’C) 

C 

Mn 

St 

Cr  Ni 

Onnui 

Pure  iron 

.178 

Cast  iron 

116 

57 

1  54 

112 

Carb  steel 

23 

64 

.124 

Carbon  steel 

1  22 

35 

.108 

Alloy  steel 

,34 

55 

.78  3.53 

Mo  39 

.079 

Cu  .05 

T1  tool  steel 

.70 

4.0 

W18.0 

038 

VI. 0 

Stainless  type  410 

.13 

1.0 

10 

1 1.5-13.5 

.057 

Stainless  type  304 

.08 

2.0 

to 

l8.0-'20.0  8.0-12.0 

.036 

a  factor  of  3.  For  EN  HI  or  AlSl  52100  steel,  generally 
used  in  the  fully  hardened  state  for  bearing  components, 
tempering  due  t(.  frictional  heating  could  result  in  as 
much  as  a  67‘r  increase  in  i  above  its  value  for  the 
as-received  fully  hardened  condition  (20). 

CONCLUSIONS 

An  analytic  expression  for  the  steady-state  intcrfacial 
temperature  field  within  an  unlubricated  (or  boundary- 
lubricated)  sliding  circular  Hertzian  contact  has  been 
derived  in  terms  (/i  the  lotai  irictionai  power,  veiocity, 
contact  radius,  the  thermal  properties  of  the  materials, 
and  the  bulk  temperature  of  each  body.  The  detailed 
analysis  has  been  presented  for  the  case  where  one  of  the 
surfaces  does  not  move  with  resfject  to  the  contact  area 
(f,  =  0),  as,  for  example,  in  the  4-ball  machine,  and 
where  v-j  >  lOa  although  the  relationships  derived 
would  hold,  to  a  lesser  ilegiec  of  accuiacy,  foi  Kiwet 
velocities.  The  basic  method  presented  in  this  paper, 
however,  can  be  applied  to  any  combination  of  c,. 
outside  of  the  range  O.'Ja,  /?  <  jej  <  Kid  /?. 

Although  not  an  exact  solution  to  the  problem,  the 
results  of  this  paper  are  more  reliable  than  previously 
reported  approximations  in  the  following  respects:  (a) 
Account  is  taken  of  the  fact  that  the  [viwer  density  is 
ellipsoidally  distributed  ttver  the  contact  circle,  (b)  .A 
factor  is  intniduced  into  the  .asymptotic  expression  for  the 
temperature  distribution  7’.(  t,  r)  within  a  moving  he.it 
source  on  a  single  surface  whiih  extends  the  range  ol 
validity  of  the  results  down  to  /i  =  Id.  (cl  It  was  shown 
that  the  interfacial  temfjerature  distribution  VU.  »•)  can 
be  approximated  bv  half  the  harmonic  mean  of  the  two 
single  surface  temperature  fields  7‘,(  t,  >•),  7',,( ».  n  and  that 
the  accuracy  of  this  appmximation  improves  with  in¬ 
creasing  veliK'ity  as  7h.  »)  approaches  /,(».  n.  The  ex- 
jiression  for  7„,,„  was  obtainetl  direiih  from  /is.  ri 

In  atldnion,  the  problem  and  its  analv'ie  solution  base 
been  extended  to  iiuhule  the  more  get  lal  cast  of  un- 
'■i|0.d  hulk  tempeiatuifs  of  the  two  diiii  .g  bodies.  This 


alters  the  heat  partition  coefficient  and  the  magnitude 
and  shape  of  the  interfacial  temperature  Held.  .As  /?  — *  oo, 
the  interfacial  temperature  approaches  the  bulk  temper¬ 
ature  of  the  moving  body. 

The  expression  derived  for  the  maximum  flash  tem¬ 
perature  7'„,„  gives  a  value  from  33%  to  38"^  greater 
(depending  on  the  value  of  B)  than  Blok’s  widely  used 
formula  (/). 
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A  Critical  Survey  of  Mathematical 
Models  for  Boundary  Lubrication 

A.  BEERBOWER 
Government  Research  laboratory 
Esso  Research  and  Engineering  Co.,  Linden,  New  Jersey  07036 


The  various  equations  representing  special  cases  of  boundary 
lubrication  are  examined  and  grouped  as  “models".  The  study 
indicates  that  there  are  no  serious  mconsistencies,  but  that  each 
model  IS  so  seriously  limited  in  its  coverage  that  large  gaps  exist 
w  coverage.  Methods  for  closure  of  these  gaps  ore  suggested  and 
numrncnlly  explored  to  a  limited  extent.  Hope  is  extended  that 
in  the  near  future  it  ivill  be  possible  to  compute  wear  rates  from 
basic  physical  pinperlies  without  actual  exfierimentation  on  the 
specific  lubricant  and  surfaces.  This  would  constitute  a  useful  new 
system  for  machine  element  design. 

INTRODUCTION 

Elastohydrodynainic  lubrication  models  fornr  a  rea¬ 
sonably  well  integrated  matrix,  but  those  related  to 
boundary  hrbrication  are  isolated  and  apparently  iiteott- 
sisteni.  However,  there  are  enough  of  them  to  justify  a 
serious  attempt  at  unification.  The  scope  of  this  project 
was  selected  to  exclude  dry  solid  film  lubricatioir  and  gas 

Prticnled  at  on  Amaricpn  Socialy  of  lybricotlon  Enginttrt  paper  ot 
the  ASM* '  ASLE  luhtlcolion  Conference  held  in  Cincinnati,  Ohio, 

October  13-15,  1970 


a  ~  iractiorial  film  delect  (for  additive  in  Model  If  O) 

.4,  cross-si  ctiotial  aie.i  of  wear  track  (cm"’) 

/),2  =  work  of  adhesion  (eres/ciii") 

fi  =  fractional  li.ise  fluid  film  defect 

H  —  electric  charging  energy  of  met,\l  (ergs/cni‘) 

Y  =  (1  +3  /»■)"  '■ 

Yy  =  Slirf.ice  free  energy  of  licjuid  extrapolated  to  0”K 
(ergs/ci.i-) 

Y|  -  sorfarr  free  etin  gs  of  solid  (ergs  cm*) 

Y;  --  snif.tce  free  energy  of  lifpiid  li  igs  etn‘) 

Yj.j  —  irtlerf.n  i.il  fri'e  <  irergv  (erg-  cni-r 
f  -  addilive  coim  tiit.irroii  in  rnoi  .ii  iioii 
r.;  cnristanrs  in  liniie  wear  e(|nations 
(j  r-.  oxygen  cometilratioti  in  Inijricaiil  (gnt  gm) 

M  rr:  l.ondott  tone  Miliihility  |iar.mieter  (eal  cm')"'* 


bearings,  and  to  include  all  liquid  lubrication  even  when 
the  liquid  is  thickened  to  a  grease  or  a  solid  additive  is 
present.  The  load  regime  is  defined  by  the  zone  A')’  in 
Fig.  1,  where  “relative  load”  may  be  taken  as  the  force, 
generalized  in  regard  to  geometry  required  to  produce 
the  wear  effects  illusiraied. 

As  stitdy  proceeded,  it  became  necessary  to  recognize 
additional  regimes.  The  first  of  the.se  is  AX.  in  which  the 
original  asperities  interfere  with  the  clastohydrodynarnic 
icgiiiu  in  OX,  Hi  Ic.iM  uiiiil  liiey  are  smootricd  out  by 
one  Of  more  “break-in”  processes.  'I'liis  “running  in" 
regime  is  transitory  in  time  but  of  vital  importance.  If 
the  proce,sses  include  ehemica!  reactions,  the  regime  of 
“corrosive  wear”  extends  from  A  to  li.  The  position  of 
B  is  highly  variable;  with  very  inert  systems  it  may  lie 
at  A',  while  witli  very  active  /:/’ (extreme  pressure)  agents 
it  shifts  all  the  way  to  Z'.  In  ordinary  systems,  it  lies 
between  .V  and  and  the  BY  regime  is  known  as  “adhe¬ 
sive  wear.”  Finally,  at  there  is  an  abrupt  transition  to 
“sculling  wear”  or  “scoring.”  Tlie  I'Z  regime  is  of  little 
interest  to  the  operator  of  successfully  lubricated  ma¬ 
chines,  but  the  exact  location  of  Y  is  of  great  interest  to 
designers. 

Sp  —  |)(«lar  solubility  piriamrlcr  (eal/rm'^l"  '* 

Sfi  Iisdrogen  boinliirg  solribilily  p.iramcirr  (eal  cm')"  ' 

d  =  sliding  dislaiice  (cm) 

rfy  “  zero  wear  .sliding  dislaiice  (cm) 

I)  =  siiifacc  rniiglmess,  peak  to  peak  (;i  in  or  ;»  in) 

<  =  dii’lecn  ic  constant  of  liibiicaiil 
K  =  heat  of  adsorption  (c.il  'mol  or  crgs/rtii') 

AF  —  dinercncc  in  heat  of  .idsot|)iion  of  addilive  (A),)  and  base 
fluid  (A.'j) 

h\,E.,  —  Young's  mndiilris  of  solids  (psi) 

G  =  nilowuble  fraction  of  for  zero  wear 
Gf,  ez  a  for  'JDIIO  passes 
A  =  Plank's  (onstani  ib  :)  x  Id  ■'  erg  ser ) 

//  Meriz  eoniael  diim  nsion  (in  oi  cm) 

X  =  Uoll/.m.m's  (onsiarrl  (1,1  x  Id  ergs,  deg  (',) 

"weal  eoeftuieni  '  for  melal-mel.il  loniael  ari'a 
X,.X._,  r-  tlierinal  eondra  livities  ((  .d  sei  -t  in-  ()) 


NoSlI  Ne.l..XTl'RF. 
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REUmvE  LOAD 

fig.  1— RaginiM  of  lubricoHon. 


.-Vs  a  useful  simplification,  it  was  decided  to  eliminate, 
at  least  for  the  present  analysis,  the  possibility  of  foreign 
abrasive  contaminants  reaching  the  wearing  surfaces.  It 
»».i\  also  as-sumed  that  the  lubricant  was  adequately  in¬ 
hibited  against  bulk  oxidation  and  that  changes  due  to 
ihrimal  decomposition  were  negligible.  However,  a  full 
f.iiige  of  atmospheres  was  accepted  as  a  necessary  part 
i  f  the  job,  with  those  more  corrosive  than  humid  air  being 
v-t  .iside  temporarily.  Also  temporarily,  the  ambient  at- 
t’l<l^phcric  pressure  was  set  at  one  atmosphere.  Thus,  the 
vojic  includes  break-in,  (/)  corrosive  w'car  due  to  lubri- 
1  ant  or  to  atmospheric  O2  +  H20,  (2)  adhesive  wear  and 
ir.-»nsitinr>  nnini  jocatiors  {3).  The  sfTect  Cl  iubricaticn  on 
fatigue,  except  as  it  marks  a  limiting  condition,  does  not 
vrem  to  be  ready  for  this  type  of  consideration. 

Mooeis 

•Strictly  speaking,  any  equation  relating  to  boundary 
lubricatioii  could  be  considered  a  “model”  and,  therefore, 
ili  sciving  t)f  separate  handling.  Fortunately,  many  of 

A',_,  =  iiiterfarial  wetting  factor 

V  =  mass  diffusivity  (cin®/sec) 

/  =  number  of  res-olutiuns  or  strokes 
=  1.  for  zero  wear  period 
/  =  I'oeHieienl  of  friction 
•  i  -.  rx|Minent  in  finite  wear  equation 
M  -•  molecular  sveight  (g/mole) 
r|.i  _.  =  I’oisson's  ratio  of  solids 
.T  minilKT  of  molecules  per  utiit  area  (cm”*) 

^  -  miiniM'r  of  p.isses 
_  N..  =  .V  for  zero  wear 

|i.irti,d  pifNsure  of  preferred  component 
/ =  lloss  pies.sure  of  metal  under  static  loading  (gni/cin*  or 
kg  mm'-) 

v  fr.ieiional  surface  coverage  by  additive 
y.  m.sximum  stress  in  Hertz  contact  area  (psi) 

V  e  be.ii  of  surfiice  formation  (tii>s/cin‘’) 

’  .  i.olms  Ilf  s|)licrical  member  (in  or  cm) 

I!  •  mill. It  g.is  constant  (cal/mole  °K) 

!•  ileiiMts  of  oil  (gm/cc) 

I'.,  ■-;  <>  at  tiO'F  tl5‘C) 

total  entropy  change  associated  with  adsorption  (cu) 


these  have  been  rounded  up  and  screened  into  sonic 
coherence  by  MacGregor  ct  al  (-f)  and  by  Rowe  {5 ).  These 
will  be  referred  to  as  Model  1  and  Model  II  rcs|}cctively. 
In  addition  to  these,  it  is  necessary  to  consider  a  third 
iT-”'dcl  which  has  not  yet  been  completely  formulated. 
Model  III  is  based  on  the  general  concept  of  irrevcrsable 
reactions.  These  include  the  formation  of  “friction  poly¬ 
mer”  a  material  generated  under  conditions  of  low  wear 
from  the  base  stocks,  apparently  by  catalytic  action  of 
freshly  sheared  metal  surfaces.  A  fourth  model  came  to 
the  writer’s  attention  too  late  for  complete  analysis  at 
present,  and  a  fifth  is  given  minor  attention  as  a  special- 
purpose  matter. 

BASIS  OF  MODEL  I 

This  model  is  the  result  of  a  great  deal  of  experimentai 
work  at  IBM,  an  account  of  which  is  given  by  Bayer  (6, 
7,  8).  It  may  be  described  as  essentially  geometric  and 
empirical,  though  a  definite  relationship  to  the  Archard 
wear  and  Palmgren  fatigue  equations  has  since  been 
established.  Only  four  ometinics  three)  “grades”  of 
lubricant  are  considered.  I'hcsc  were  selected  to  cover  the 
range  of  oils  commonly  encountered,  along  with  the  dr\’ 
situation.  Of  course,  this  constraint  is  very  irksome  to 
those  concerned  with  lubricant  development,  especially 
since  no  means  is  provided  for  interpolation  or  extrapo¬ 
lation  to  other  grades  of  oil  or  grease  beyond  assuming 
that  the  changi  will  mn  be  rui  ilie  beiici.  Anoiiier  serious 
limitation  is  that  it  is  not  possible  to  adjust  the  mode' 
for  environmental  conditions,  these:  being  preselected  as 
air  at  one  atinosiihcre,  30%  relative  humidity,  and  a 
temperature  of  about  22  C.  It  is  assumed  that  .speeds  will 
be  kept  low  enough  that  surface  temperatures  will  not 
be  significantly  above  ambient,  so  there  is  no  adjustability 
for  speed.  A  fourth  limitation  is  that  parameters  nuisi  be 
looked  up  for  specific  material  pairs  or  couples.  'I’liere  is 


S  =  distance  traveled  per  revolution  or  stroke  (in  or  cm) 
=  iiiaximuni  shear  stress  in  coiilart  area  (psi) 

—  yield  stress  in  shear  of  solid  (psi) 
t  =  time  (see) 

<0  =  vibrational  time  of  .adsorbed  molecule  (sec) 
y'p  =:  ratio  of  tg  values  for  additive  and  base  fluid 
7j  =  critical  temperatures  of  lubricant  (°K) 

7j  r=  bulk  temperature  of  lubricant  ("K) 

T„  =:  melting  point  of  lubricant  (°K) 

7’,  =  temperature  of  surface  ('K) 

7]  =  transition  temperature  at  scufTiiig  (°K) 

0  =  contact  angle  (deg) 

U  =  sliding  velocity  (cm/sec) 

{/,  =  lubricant  flow  velocity,  average,  into  gap  (cm/sec) 
y  =  wear  volumes  (cm'’) 

=  molar  volume  (em'’/mol) 

W  =  load  (,gms  or  lbs) 

III  =  width  of  wear  i.-:u  k  (cm  or  in) 

U'u  =  «'  for  end  of  zero  wear  period 

X  =  ratio  of  iitolnr  di.ameieis 

,V  =  diameter  of  adsorbed  molecule  (i  tn) 

{  =  dearamc  of  betitiiig  (cm) 
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a  wide  selection  of  alloys,  and  also  plastics  for  one  mem¬ 
ber  of  the  pair,  but  only  three  metals  arc  available  for 
the  other  member:  52100  steel,  502  stainless  steel  and  a 
65/35  brass.  It  is  possible  to  make  limited  adjustments 
in  this  part  of  the  model;  for  instance,  the  same  alloy 
heat-treated  to  a  dilferent  hardness,  or  plated  on  a  differ¬ 
ent  metal  substrate  can  be  handled.  However,  this  re¬ 
quires  the  assumption  that  the  coelTicient  of  friction  is 
the  same  as  for  the  material  listed,  which  is  hard  to  justify. 

'I’he  reasons  for  these  limitations  lie  pe.rtly  in  the  history 
of  the  model.  It  vvas  developed  primarily  for  use  in  de¬ 
signing  business  machines  to  opeiate  in  air-conditioned 
olfices.  Release  to  the  general  engineering  public  was  a 
later  consideration,  and  apparently  suliicient  warnings  of 
its  limited  applicability  were  not  included.  .-Xs  a  result, 
lubrication  engineers  have  teixled  to  criticize  Model  I  for 
not  serving  purposes  clearly  outside  its  scope,  rather  than 
looking  to  sec  m  what  directions  it  needs  to  be  extended. 
The  authors  (f>)  do  not  plan  to  extend  the  model  as  it 
already  serves  the  original  purpose  at  IBM.  However, 
they  suggested  that  anyone  interested  in  doing  so  could 
prepare  tables  for  other  temperatures  and  humidities. 
Plastics  in  particular  need  such  special  handling. 

The  limited  version  available  lias  two  levels.  Model  lA 


Use  of  the  model  is  by  stepwise  calculafion.  FiiNt,  the 
contact  stress  is  calculated  from  tbe  Hertz  formula 
appropriate  to  the  geometry,  the  loan,  the  radii  of  curva¬ 
ture,  the  Young’s  modulus  ?nd  Poisson’s  ratio  of  the 
materials  and  the  ccefficicnt  of  friction  of  the  system,  F. 
A  stress  concentration  factor  may  be  used  to  correct  the 
Hertz  formula  for  sharp  corners,  and  non-Hertz  cases  can 
also  be  handled.  Values  of  F  must  be  looked  up  from  the 
same  table  as  Gg,  obtained  experimentally  or  estimated 
by  analogy  from  listed  systems. 

The  number  of  passes  for  zero  wear  (Nq)  in  passes  may 
be  estimated  by  combining  Inequality  [1]  and  Equa¬ 
tion  (2). 


However,  the  lifetime  of  each  member  must  be  con¬ 
sidered  separately.  For  the  "fixed  spot”  member,  on  which 
the  contact  spot  does  not  move,  the  zero-wear  travel  dg 
in  inches  of  relative  travel  is  given  by 

do  =  SLo  =  HMo  [4) 


is  for  design  at  the  “zero  wear”  level,  defined  as  ilamage 
low  enough  not  to  change  the  original  surface  roughness. 
This  may  be  considered  to  be  a  means  of  locating  point 
A  in  Fig.  1.  Model  IB  is  for  finite  wear,  and  is  further 
subdivided  into  two  types  of  wear  apparently  corre¬ 
sponding  to  regimes  XB  and  BY. 

Both  models  are  tied  to  a  concept  of  unit  travel  called 
a  “pass”  and  defined  as  tbe  distance  of  sliding  ctiual  to 
the  Hertz  contact  area  dimension  in  the  direction  of 
motion. 


ANALYSIS  OF  MODEL  lA 


This  model  uses  the  criterion  for  zero  wear  that  the 
maximum  shear  stress  in  the  contact  area  must  be 
smaller  than  a  certain  traction,  6',  of  the  yield  stress  in 
shear  of  the  material,  t^.  Hence,  the  following  inequality 
defines  “zero  wear  conditions”: 


'him  ^  ^'^1/ 


ni 

y. 

When  the  number  of  passes  (iV)  is  2fl()00(i  has  the 
reference  value  Cg.  Gg  can  have  one  of  three  possible 
values,  depending  on  tbe  materials  and  lubricant,  for  full 
hydrodynamic  lubrication  Gg  -  1.00.  For  boundary 
lubrication  Gg  -  0  54  designates  systems  with  low  sus¬ 
ceptibility  to  tr.insfer  and  Gg  —  0.20  for  tliose  with  high 
susceptibility.  No  systems  showing  Gg  values  intermediate 
between  0.20  and  0.54  have  ever  been  delected,  indicating 
that  these  represent  diircrcnt  sub-regimes.  V'alues  of  Gg 
iire  tabulated  (•/)  for  over  500  systems,  and  for  unlisted 


where  //  is  the  length  of  the  Hertz  contact  spot  in  the 
direction  of  travel,  .S'  the  distance  traveled  per  revolution 
or  stroke  and  the  lifetime  in  revolutions  or  strokes. 

For  the  “moving  spot”  member,  on  which  the  contact 
spot  keeps  ntoving, 

do  =  SLo  =  SNo  (5) 

The  lower  of  these  d^  values  is  taken  as  the  lifetime  of 
the  system.  Failure  at  this  point  is  believed  to  be  by 
fatigue  (8). 

It  is  not  possible  to  do  a  great  deal  of  mathematical 
testing  on  Model  I.\.  as  the  only  inputs  arc  load,  geome¬ 
try,  properties  of  the  two  materials  and  F.  The  value  of 
Gg  is  automatically  selected  from  the  tables  (7)  or  esti¬ 
mated  to  be  0.20.  For  an  example,  the  ball-on-planc 
geometry  of  Rowe  (,'5)  and  Bayer  (6)  may  be  put  through 
some  com|)utations.  Unfortunately,  no  fixed  F  value  for 
Rowe’s  copper  (hcmi-)sphcre  is  available,  but  the  lange 
0.26  to  0.57  can  be  reconstituted  from  his  paper.  He  gives 
hardness  values  for  his  copper  and  steel  of  88  and  220 
kg/inm-  resijcciivcl) . 

The  equations  for  the  sphere  are  (13) 

t„m  =  -2*',)Vl6]-bF^  16] 

where  qo  is  the  maximum  stress  in  the  Hertz  contact 
area  and  »•,  is  Poisson’s  ratio  for  the  metal  of  the  sphere. 
For  this  geometry 
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Wherc  H'  is  ihc  load,  E,,  and  arc  Young’s  modulus 
for  sphere  and  plane,  and  r  is  the  radius  of  the  sphere 
(0.312  cm  or  0.123  in). 

By  definition 


zero-wear  condition  described  in  Equation  (3).  The  wear 
on  the  first  surface  is  predict  'd  by  a  differential  equation 
which  is  basically  a  statement  that  wear  is  a  function  of 
the  energy  dissipated  per  pass: 


[81 

ForCoppcr,£j  =  16.5  X  10®, i',  =  0.34, =  !2,800pse; 
for  Steel,  =  30  X  10®,  V2  =  0.30;  System  =  0.20 
(the  conservative  assumption). 

With  these  inserted,  the  minimum  zero  wear  travel  of 
this  system  in  inches  is 

.  8.68  X  10-*-'^  ,Q, 

®  (0.0064  +  ^  ' 

This  was  evaluated  for  the  loads  and  F  values  shown 
in  Table  I.  The  loads  arc  considerably  below  Rowe’s 
range,  to  avoid  the  criterion  <  //  ( IG).  This  is  not 
clearly  expressed  in  the  published  versions  {4,7)  though 
it  is  implicit  in  inequality  ( I  ].  The  /i  values  cover  the  entire 
-ange  from  MacGregor’s  lowest  to  Rowe’s  highest. 

We  may  conclude  that  Model  lA  is  useful  in  locating 
point  A  in  Fig.  1  but  needs  generalization.  The  extreme 
sensitivity  of  to  F  creates  a  problem,  and  a  more  so¬ 
phisticated  method  for  estimating  F  would  be  highly 
desirable. 


d[A,hr,^Wf'^]  =  CdN  [10] 


This  equation  can  be  integrated  only  for  specific  situa¬ 
tions,  of  which  only  the  Rowe  example  used  above  will 
be  illustrated.  Converting  d,,  the  cross-sectional  area  of 
the  wear  track,  to  the  more  useful  wear  track  width  (in 
this  case,  equal  to  ball  scar  diameter)  is  done  by  an 
approximation  based  on  a  triangular  area  formula 


""  16r 


[HI 


Wear  on  the  ball  results  in  a  flat-to-flat  geometry,  so 
that 


4  M'(0.25  +  F2)®'® 


TTW 


[12] 


Making  these  substitutions  and  integrating  for  ball 
wear. 


w  =  C(Lr5)0'  '®/’®®3''(0.25  +  F^)®--’®®  [  1 3] 


EVALUATION  OF  MODEL  IB 

This  model  (7)  requires  an  immediate  decision  as  to 
the  type  of  wear  anticipated,  not  by  an  analytic  or  com¬ 
puter  process,  but  based  on  the  experience  of  the  user  or 
on  a  quick  simulative  experiment.  There  arc  two  possible 
modes;  one  involve:-  high  transfer  of  material  from  one 
surface  to  the  other.  If  this  corresponde.l  to  regime  YZ 
in  Fig.  1,  it  would  be  outside  the  scope  of  this  study. 
However,  as  shown  under  Model  IlIA,  it  corresponds  to 
the  “Archard  wear  law”  under  conditions  obviously  re¬ 
lated  to  corrosive  wear.  The  other,  with  little  or  no  tran.s- 
fer,  can  be  analysed  by  using  many  of  the  factors  already 
considered  in  Model  lA.  It  must  be  recognized  at  the  start 
that  only  one  of  the  surfaces  is  expected  to  show  finite 
wear,  while  the  i  her  one  is  expected  to  remain  in  the 


where  C  includes  the  origiital  C'  limes  several  numerical 
constants.  Neither  C  nor  Ccan  be  developed  directly  from 
material  properties. 

By  eslabli.shing  a  formiilization  of  the  definition  of  zero 
wear,  it  is  possible  to  establish  a  useful  relationship  be¬ 
tween  Models  lA  and  111.  Ibis  is  done  by  defining  tin- 
depth  of  the  “zero”  wear  scar  as  half  the  surface  rough¬ 
ness,  D.  From  this  and  the  triangle  approximation,  the 
zero  wear  scar  width 

T„  =  {8,D)^'-^  1 14) 

Assuming  that  all  the  other  factors  in  Equation  (M) 
remain  constant 

(15) 


Taiu.f.  I 

-  PREDlCrri'  U  ZERO  WEAR 

rR.AVEE  (in)  FOR  COPPER 

IIAI.I.  ON 

STEEL  DESK 

CotrnciFNT 
or  FmcTioN 

(F); 

.08 

.16 

.24 

.32 

.40 

.48 

.56 

64 

Load  (If) 

Poll.NOS 

.«)! 

280. 

4.63 

0.205 

1,88  X  10  ’ 

2.78  X  10  J  1 

5  I18  X  10 

1  4(i  X  10  ' 

4.50  X  10  1 

.004 

7,10 

0.115 

2.08  X  10  ■'  1 

Mill  X  10  « 

6.80  X  to  » 

1.  41  X  10 

3  03  X  lO  '- 

1  Hi  X  10  « 

.015 

O.'.W) 

3.38  X  10  >  r 

1..50  X  10  ♦ 

1.37  X  10  > 

2.03  X  10  * 

1,15  X  10  ’ 

1  07  X  10  1 

t.’O  X  10  ’ 

.005 

4,1‘t  X  10 

2J  6.78  X  10  » 

3,00  X  10  •* 

2  75  X  10  » 

4,0t.  X  10  " 

8  31  X  10  !• 

2  14  X  10  » 

6,58  X  10  '>> 

.250 

1.15  X  10 

«  1.86  >;  10  •« 

8.24  X  10  " 

7.57  X  10  » 

1.12  X  10  ’> 

2.2'l  X  10 

5  00  X  10  " 

1.81  X  HI  0 

Note:  All  vslues  lirlnw  llir  line  are  unrrliulilr  liy  tlir  erilerioii  tfo  <  II  (Ki). 
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whcrc  m  is  0.1 18  for  ball  wear  and  0.  KIO  for  plane  wear. 
In  the  hitfh  transfer  rei'inie.  these  c.xponcnts  are  0.25  and 
0.33  respectively.  As  shown  under  Model  IIIA,  m  —  0.25 
is  equivalent  to  the  “.Archard  wear  law.” 

It  has  been  demonstrated  that  /.^  is  independent  of 
surface  finish  (/•/)  it  is  proper  to  substitute  the  value 
obtained  from  Equation  (!1)  into  Equation  [13]  as  this 
point  renresents  the  intersection  of  the  zero  wear  line  with 
the  finne  wear  curve.  This  is  a  discontinuity,  where  the 
ball  ceases  to  have  a  spherical  contact  on  the  plane. 
Another  discontinuity  comes  when  the  second  member 
starts  to  wear,  and  Equation  (13|  is  not  applicable  above 
this  point. 

It  is  evident  that  Ecpialion  [13]  shows  a  much  greater 
dcpetidcnce  of  wear  on  surface  finish  than  most  investi¬ 
gators  find  in  the  boundary  regimes  HO).  Bayer's  tying 
D  to  the  coarser  surface  (7)  has  been  reconsidered  and 
the  writer  believes  it  best  to  use  the  roughness  of  the 
surface  e.xpected  to  first  show  wear.  This  is  necessary  by 
the  definition  [I  f]  and  is  compatible  with  the  general 
nature  of  adhesive  wear.  It  is  important  to  use  peak-;  .■ 
peak  values  for  D.  If  only  root -mean-square  values  are 
available,  they  should  be  multiplied  by  2  to  con¬ 
vert  them. 

I'hcre  is  a  case  in  which  the  roughness  of  the  non¬ 
wearing  surface  is  significant,  but  it  is  so  special  as  to  be 
outside  the  scope  of  this  study.  It  is  usually  the  result  of 
very  poor  design  in  which  the  harder  member  acts  as  a 
flic  on  the  softer  one.  This  is  known  as  “two-body  abrasive 
w'ear  ,  aiul  is  set  aside  along  witli  the  similar  case  in 
which  hai  d  foreign  matter  embeds  in  one  surface  to  create 
the  “file.” 

The  result  of  computation  using  Model  IB  with  the 
sitme  conditions  as  for  zero  wear  in  Table  1  are  shown 


in  Table  2.  Rowe’s  copper  sphere  roughness  of  I  micro- 
inch  RMS  was  converted  to  2.8  microinches  peak-to- 
peak.  To  provide  results  comparable  to  him,  the  program  ^ 
was  set  up  to  read  out  in  ^'/d-  using  the  proper  geometri¬ 
cal  formula  to  convert  T  to  V.  The  d  values  were  selected 
to  cover  Rowe’s  ranges  of  3  to  20  hours  at  0.1  to  100 
cm/scc. 

The  results  compare  very:  favorably  with  Rowe’s  Fig. 

4,  where  the  lightest  load  gives  a  value  of  I'/rffry  of 
1.2  X  10  '  -’  cc/ern-g.  Transforming  Table  11  to  the  same 
units,  we  find  that  for  F  =  0.24  and  (7=10*  cm  (3900 
in),  the  values  of  this  function  vary  from  1.3  X  lO"'"*  at 
the  lightest  load  to  2.3  X  10'  *”  at  the  first  load  beyond 
the  t/g  <  //  criterion.  Similar  calculations  for  m  =  0.250 
gave  values  of  2  X  10'  "  and  3  X  10'**,  respectively, 
indicating  that  m  ■—  0.1 18  was  the  proper  choice  despite 
Rowe's  assumption  of  the  “Archard  law.” 


BASIS  OF  MODEL  II 

The  Rowe  model  was  presented  in  three  papers  of 
increasing  degrees  of  sophistication  II,  12).  It  is  based 
on  reversible  adsorption.  The  basic  equation  is  an  adapta¬ 
tion  ol  the  “Archard  wear  law” 


Pm 


[16] 


where  is  a  dimensionless  wear  cocflicieiit  descriptive 
of  the  surfaces,  W  the  load  and  the  flow  pressure  under 
static  loading.  The  other  factors  are  obtained  as  follows: 


Y  =  (1  -f  (17) 


Table  2-  PREDIC'l TU  WEAR  R.ATES  FOR  MODEL  IB  (cm\cni) 
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Tbc  conslant  3  is  from  the  von  Miscs  equation  for  two 
dimensional  plastic  flow.  Spurr  (13)  has  further  justified 
this  value. 

The  fractional  surface  film  defect  is 

\ -Xtxp(-E/RT,)-\ 

where  A'  is  the  diameter  of  the  spot  associated  with  an 
adsorbed  molecule,  E  the  heat  of  adsorption,  R  the  gas 
constant  and  U  the  sliding  velocity. 

The  surface  temperature  is 

(WP 

1.040  x  \Q-^FUy.  [19] 

where  Tj  is  the  bulk  labricant  temperature  and  and 
the  thermal  conductivities  of  both  surfaces.  This  is  the 
well-acceptcd  “surface  temperature  rise’’  equation. 

The  fundamental  time  of  vibration  of  an  adsorbed 
molecule  is 


/MV  2^3v0.ri 

lo  =  4.75  X  [20] 

where  Af  is  the  molecular  weight  of  the  lubricant,  its 
molecular  volume  and  its  melting  point. 


ANALYSIS  OF  MODEL  ilA 


It  must  be  mentioned  that  de  Boer  (14)  has  recently  taken 
considerable  exception  to  the  practice  of  regarding  tg  as 
the  vibration  time.  He  finds  that  for  adsorbed  gases, 
1^  =  h/kT  where  h  is  Plank’s  and  k  Boltzman’s  constant, 
liquids  show  lower  values  by  a  multiplier  which  is  very 
difficult  to  predict. 

Values  of  F  cannot  be  predicted  with  any  degree  of 
certainty.  Many  values  can  be  picked  up  from  Mac¬ 
Gregor’s  Table  III-3.  Fortunately,  values  of  y  are  not 
subject  to  even  as  much  variation  as  those  of  F,  which 
range  only  from  0.08  to  0.29  in  MacGregor's  lubricated 
systems.  (Values  up  to  0.39  on  Be-Cu  seem  to  be  a  rare 
exception.)  Hence,  y  =  1 .06  covers  all  his  data  with  ±5% 
error,  though  Rowe  found  1.1  <  y  <[  1.4  for  his  system. 

To  predict  the  effects  of  this  variation  on  T,,  we  may 
use  Rowe’s  example  of  copper  on  steel  at  10  cm/sec, 
4000  g  and  31 PK.  The  temperature  rise  would  be  only 
1.5  to  5.5°K  over  the  range  of  y  from  0.08  to  0.29.  At  his 
maximum  speed  of  100  cm/sec,  the  rise  would  be  15  to 
55°K,  so  that  T,  vyould  be  346  ±  17.3°K.  Thus,  again 
a  fixed  value  of  0.175  would  only  lesult  in  about 
±5%  error  Wc  cannot  be  complacent  about  this,  how¬ 
ever,  as  processing  these  values  through  Eq.  [18]  results 
in  a  seven-fold  range  in  n,  the  fraction  of  lubricant  film 
defects,  and  hence  in  wear  rate. 

Some  hope  for  a  more  exact  input  may  be  derived  from 
Spurr  (15),  who  arrived  at  an  equation 

,  It  rosfy  entfi  \ 

'  >  _ ii 


All  of  these  input  values  can  be  estimated  or  arc  known 
quite  pa-cisely.  Rowe  tested  this  equation  by  back- 
calculating  from  wear  experiments,  but  the  present  plan 
is  to  use  it  for  design  purposes  and  then  judge  the  results 
against  both  general  and  anomalous  experience.  For  this 
purpose  one  may  regard  II  ',  U,Af,  and  <0  as  preselected 
parameters.  /*„,  is  obtainable  from  the  Diamond  Pyramid 
Hardness,  Vickers  Hardness  Number,  etc.  or  from  Mac¬ 
Gregor’s  Table  111-4  (4),  k„,  has  a  theoretical  value 
for  hemispherical  wear  particles  of  0.33,  shown  by  Rowe 
to  be  of  the  right  order  of  magnitude. 

'Fhc  values  of  ,V  arc  not  so  generally  available,  and  it 
is  proposed  to  use  a  sitnpitr  parameter.  (GV,.  which 
haT pr(>vcSnsrlicT>f  grcaiTvaluc  i h’pretlTcnng  surface  free 
energies.  This  is  the  diameter  of  the  sphere  associated  with 
the  surface  molecule.  Hence, 

,  ^  /  -exp(-jf//?7.)  \ 

®  U(M/TJ'>'')  ^ 

Values  of  T„  ate  readily  available  for  pine  compounds, 
but  for  mixtures  such  as  lubricating  oils  they  simply  do 
mil  exist.  (The  pour  point  only  Isidk  .ties  that  enough  wax 
ha%  coagulated  to  gel  the  lii|uid.)  Heme,  in  such  cases 
a  generalized  meliino  point  ba.sed  on  the  critical  or  pseu- 
dm-riiicai  poin;  will  bt  usctl, 

:=  0.40  T,  (22) 


where  is  the  yield  stress  in  shear  of  the  softer  material, 
and  arc  the  contact  angles  of  water  on  the  rider  and 
specimen,  or  “fixed”  and  “moving  spot”  rcs|)ectively, 
when  lubricated  with  a  partial  film  of  stearic  acid.  His 
measurement  of  by  Vickers  hardness  number  is  pre¬ 
sumably  equivalent  to  Rowe’s  Diamond  Pyramid  Haid- 
ness  and  MacGregor's  Microhardness,  If  so,  is  ibiaincd 
from  by  MacGregor’s  Fig.  I1I-3. 

Spurr’s  reasoning  points  up  a  weak  spot  in  Rowe’s,  in 
that  the  latter  proceeded  as  if  E  were  the  same  for  copjicr 
and  iron,  and  that  only  the  film  on  one  side  of  the  contact 
need  be  considered.  To  fulfill  Rowe's  definition  of  o, 
which  is  the  same  as  Spurr’s  for  cos  0^  cos  0,,  it  can  be 
proved  by  elementary  probability  theory  that  it  must  be 
the  product  of  a,  and  o,  which  arc  equal  to  the  cos 
and  cos  8,  respectively.  'I'hcse  can  be  calculated  by 
Equation  |18]  or  [21]  by  using  the  appropriate  E,  and 
A',,  the  heats  of  adsorption  on  rider  and  speciinc'’  While 
this  change  in  definition  will  nut  affect  cases  near  n  =  I. 
at  the  level  n  <  0.01  which  Rowe  indicated  to  be  of  most 
interest  it  will  produce  a  drastic  change.  In  effect,  exp 
(-Aw/iir,)  is  replacet!  byexpl  -(E,  -♦■  A',),''A’7‘,)wliiifi  for 
similar  metals  iscxpi  ~‘'2E:R'/\).  Thus,  Rowe's  eminently 
stuisfactoiy  value  t»f  =  I  I7(KI  is  cut  in  half.  The  alter¬ 
native  is  Itt  accept  a  geometric  mean  definition  n  = 
\/n/n,  which  wnultl  pieseive  the  value  t»f  E  at  some 
ex|X'nse  o,  logic.  Both  alternatises  will  be  explored  nu- 
fiierically. 
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Obviously,  the  simultaneous  solution  of  Equations  [19, 
21  and  23]  is  beyond  an  algebraic  manipulation,  but  it 
is  certain  that  an  iterative  computer  program  can  be  set 
up  to  handle  it. 

Evaluation  of  E  is  the  most  difficult  task,  and  also  the 
most  important  in  testing  Model  II.  Very  few  values  arc 
available  in  the  literature,  and  Rowe  was  forced  to  com¬ 
pare  his  experimental  heats  (or  energies)  of  adsorption 
with  other  authors’  works  of  adhesion.  These  arc  not  the 
same,  though  the  40%  or  so  difference  may  be  trivial 
compared  to  other  uncertainties. 

The  heat  of  adsorption  is  equal  to  the  energy  of  ad¬ 
hesion,  and  according  to  PhilippofT  {!€) 


(1  - 

=  truvvvw 


(24) 


where  A 12  is  the  work  of  adhesion.  A',,  is  a  relationship 
defined  by 


K 


12 


£i2\  _  ,  _  (Yl  -  Y12) 

>2  /  >2 


[25] 


and  is  numerically  equal  to  cos  0  until  it  exceeds  1.900, 
Yl  and  y.^  the  surface  free  energies  of  the  solid  and  liquid, 
and  Yj2  their  interfacial  free  energy.  Other  definitions  arc: 


n  fa.- 


1261 


12 


I  12 


..  n 

12V 


'M2' 


I ‘>11 
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The  important  difference  between  this  treatment  and 
that  used  by  many  aut!4ors  is  that  A‘,2  is  able  to  assninc 
values  greater  than  1. 1'his  avoids  the  usual  trivial  solu¬ 
tion  /1,2  =  2y2.  which  meiely  states  that  the  joint  will 
part  between  lititiid  layers  rather  than  at  the  interface. 
It  should  be  noted  that  the  Fhili|)poir  e()uations  do  not 
provide  a  value  of  /•.'  from  the  values  of  y,  and  y-.*  <'ven 
when  full  wetting  (A‘  >  1.)  exists.  I'hus,  a  further  meas¬ 
urement  of  A’,..i  or  Y|.'  is  necessary, 

A  solution  to  the  y,2  problem  is  advcK-ated  by  Fowkes 
(/7),  who  sets  up  as  his  basic  equation 

fi2  =  Yl  +  Y2  --  2  Vy/y;.^  |28| 

where  y,**  and  y.A  are  the  surface-free  energies  of  solid 
and  litpiid  due  to  "(lis|)ersion"  forces  alone.  He  tested  this 
against  interfat  ial  free  energies  l>etween  mercury  and 
watei ,  and  found  need  ftir  a  correction  for  polarity.  TIicit 
is  some  question  of  the  validity  of  this  approach,  but  it 
will  be  usetl  as  one  alternative.  It  is  es()eeiallv  useful  in 
a  modified  form  using  Hansen's  (/>'()  partial  solubility 
patameters,  in  whn  h 

y/  72  tl,(l7l.‘>  !•„'  'ifirf-  4  V’>  l‘^l 

where  and  fij,*'  .iie  the  lontiiliutions  to  the  ein  rgv  of 
vapoi i/,iiion  pit  unit  \ohiitie  bv  the  l.ondon  ,uid  |iolar 
loiifv,  ies|H'ciiseK  l  or  the  metals 


^  0.0715(A//,  -  RT,) 

Yl  — 


^.^2/3 


[30] 


where  A//,  is  the  heat  of  sublimation  of  the  metal.  This 
is  the  same  as  Yj  for  all  cubic  metals,  but  not  for  mercury, 
antimony,  zinc,  cadmium,  or  magnesium,  which  require 
special  handling.  Values  for  Yj  are  listed  by  Grosse  (/9) 
for  the  liquid  metals;  these  must  be  multiplied  by  1.14 
to  obtain  Yj  for  the  solid,  to  correct  for  heat  and  volume 
change  of  fusion.  A  very  interesting  but  less  accurate 
procedure  is  based  on  Rabinowicz’  (20)  observation  that 
Yj  is  proportional  to  again  excluding  the  above 

anomalous  metals.  From  this  it  may  be  deduced  that 


A,2^y2  +  ^OOiPj'^-P„2''^) 


[31] 


where  P,n2  is  the  hardness  of  the  metal  measured  while 
it  is  submerged  in  the  liquid.  While  this  is  dimensionally 
unbalanced  it  offers  a  most  convenient  experimental 
method  for  djj,  never  explicitly  published  though  im¬ 
plied  by  Iinaiiaka  (21). 

To  obtain  E,  it  is  also  necessary  to  have  the  entropy 
(dAyJd'l  ).  Foriunately,  this  can  be  obtained  from  the 
relation  cited  by  Duga  (22): 


7V1.2 


[32] 


where  y^  is  the  surface  free  energy  extrapolated  to  0°K. 


n‘>i 

,.,_J 


[28]  and  [27],  dilferentiating  and  substituting  into  [24] 
is  cuiiilicrsomc  but  workable.  The  Justification  for  Equa¬ 
tions  [29],  [30],  and  [32]  is  unpublished  but  has  been 
verified  for  mercury. 

An  alternative  approach  has  been  developed  by  Kara- 
shaev  (23),  who  uses  the  simpler  formula 


B(i  -  1) 
(<  +  2) 


[33] 


where  c  is  the  dielectric  constant  of  the  liquid  and  H  is 
a  constant  depending  only  on  the  electronic  structure  of 
the  metal.  A  method  for  calculating  H  fnim  the  contact 
potential  is  given  by  Zadumkin  (24),  Differentiating 
E<']uation  [33],  simplil  .  ing  by  the  npproximution  <-  ^  < 
and  assuming  that  the  temperature  tle|H.‘ndencc  of  E  is 
entirely  due  to  llieiinal  expansion. 


'iP  «/l'„  \ 

t  VrJ/;) 


[34] 


Equalic..  [33]  has  been  verifietl  by  Karnshaev  for  s(H'- 
rially  purified  gallium  at  2U°G  as  well  as  the  anomalous 
merciiiy.  and  the  lesults  arc  in  an  cird  with  I’hilippnlf's 
predictions  (/.'i).  However,  F.ipiatioii  ('’1)  gives  .1  l.irgei 
value  tif /•.'  than  Rowe's  experiments  fit,  so  ih.il  it  ii.iist 
lie  used  with  c,mlion.  This  is  piobably  due  to  neglei  t  of 
a  tem|ieratme  <le|M  ndeiue  of  //,  whit  h  preMimably  tan 
Ih'  ctiiiectetl 
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.  Computation  of  a  wear  rate  by  Model  IIA  is  somewhat 
complicated  by  the  iterative  procedure  for  p,  and  the 
necessity  to  decide  between  the  two  possible  definitions 
of  o  =  a/i,  or  yafi,  ■  Another  decision  must  be  reached 
on  the  use  of  Equation  (34)  to  replace  the  more  cumber¬ 
some  Equations  [24,  27,  28,  and  32],  or  further  develop¬ 
ment  of  the  empirical  Equation  |3I],  to  obtain  E.  Obvi¬ 
ously,  such  decisions  must  be  made  on  a  performance 
basis,  and  Bayer’s  data  (7)  on  eight  metal  combinations 
and  one  plastic/metal  pair  appears  <^uitable  for  such  test¬ 
ing.  Additional  testing  might  be  based  on  using  Equation 
(23)  in  comparison  with  MacGregor’s  tables  of  F  (4). 
Once  these  decisions  are  made,  there  is  no  shortage  of 
any  of  the  required  data,  listed  in  Table  3,  on  a  wide 
variety  of  surfaces  and  lubricants. 


Table  3-Ii\PUT  D.\TA  FOR  MODEL  IIA 

Eii(;inerring  Data: 

W,  U,  r,.  F\  r,. 

Lubrli anr  Ds'a- 

M,  F*,  h,  y.  T'o  <*. 

Metals  Data: 

y,  B,  MI,\  I'**. 

(  Two  sets  required) 

*  Knlundjni,  for  ahrrnaie  paths. 


MODEL  MB 

Rowe  also  published  (//)  the  appiication  of  Model  IIA 
to  the  lubrication  of  graphite  by  gases  and  tfannrs_  '.'hirh 
falls  outside  the  scope  of  this  paper  and  so  will  not  be 
further  discus.sed. 

ANALYSIS  OF  MODEL  IIC 

Model  IIA  was  designed  only  for  pure  liquids,  and  IIB 
tor  pure  gases.  These  limitations  were  t  oved  by  dcvel- 
ttpment  of  the  idea  of  temporary  re.sidencc  of  both  addi¬ 
tive  and  base  fluid  on  the  ntetal  surface  in  a  dynamic 
(t|uilibrium.  If  o  and  ft  represent  the  fractional  film 
defects  of  the  area  occupied  at  rest  by  the  two  species, 
respectively 

i'  =  -Ka  -  ,35j 


tslieie  rt  is  the  area  fraction  occupied  at  rest  by  additive, 
>'ui«  e  e.ich  of  these  develops  a  film  defect  related  to  its 
lie.u  of  adsorption  as  in  Equation  (18),  this  becomes 


'shetr  Sh  is  tin*  difference  in  heats  of  adsoqrtion,  x  the 
r.iiio  Ilf  intilecular  aieas.  f  ^  the  ratio  of  vibration  times, 
.\  tin-  tiiolei  ulat  iitea  of  tlie  base  fluid,  and  /y  its  \  iSiaiion 
lime 

It  is  nusv  necessary  to  evaluate  and  Rowt  cftiwe  one 
"f  the  several  |Missiblc  wasi.  He  set  up  an  equilibrium 

•4 


constant  for  the  adsorption-desorption  processes  of  both 
species  wh'  ..  will  eventually  result  in  <>  satisfying  the 
Gibbs-Duhem  requirement  of  minimum  free  energy  in 
the  system.  This  results  in  a  special  expression  for  x  =  1 , 
relating  the  wear  rate  to  friction  function  ratio  (f  /yrf) 
for  a  lubricant  with  concentration  c  of  an  additive,  to 
the  same  ratio  (f7y<f)fc  for  the  base  fluid; 

/  j,  \  _  ^-AC/ar.  r /_p  \  _  /  J:\  1  ,  (Z\ 

\yd}^  L\yrf/j  VYd'/J  \yd}^  /'o 

(37) 

where  A5  is  the  entropy  change  of  the  entire  system. 

For  the  case  where  x  values  arc  not  similar,  Rowe  ran 
into  more  difficulty  in  eliminating  p.  This  arose  largely 
because  his  goals  included  evaluating  .v  from  wear  data. 
If  the  approach  shown  under  .Vlodel  IIA  is  taken,  this 
problem  is  minimized  since  x  simply  equals  the  c  ’^c  root 
of  the  ratio  the  molar  volumes  for  additive  and  oase  fluid. 
In  addition,  /'y  becomes  the  ratio  of  liie  (Tifical  tempera¬ 
tures  of  the  two  species.  When,  as  often  happetis,  'I\  for 
the  additive  is  not  known,  it  can  usually  be  estimated  by 
the  Lyderson  method  as  discussed  in  Reid  and  Slienvood 
(25). 

It  is  also  quite  possible  to  avoid  these  difficulties  by 
precalculating  the  composition  of  the  adsorbed  layer  at 
equilibrium.  Everett  (26)  provides  means  for  doing  this 
from  activity  coefficients,  which  can  be  derived  from 
solubility  paiamciers  as  described  by  Frausnitz  (27). 
Another  approach  is  to  apply  the  necess.iry  x  correc  tions 
to  Equation  [37]  and  to  evaluate  the  entropy  change  from 
Equation  (32)  or  (34).  Grozek  (26)  gives  a  procedure  for 
the  direct  determination  of  in  a  flow  microcalorimcter. 

In  any  case,  the  tools  arc  available  to  perform  the 
necessary  manipulation  on  binary  or  muliicompcjnent 
systems.  The  choice  of  going  into  Model  IIC  with  correc¬ 
tions  for*  ^  1  or  into  Model  1I.\  with  a  weighted  aver¬ 
age  value  of  /:  would  be  largely  one  of  convenience. 

One  factor  completely  ignored  in  Model  II  is  the  fact 
that  pure  metals  air  rarely,  if  ever,  encountered  in  lubri¬ 
cation  engineering.  This  is  not  particularly  hard  to  deal 
with  in  the  case  of  solid-solution  alloys  Buckley  has 
shown  that  even  I  atomic  percent  or  less  of  alumiiuim 
alloyed  in  roppc*r  or  iron  will  difluse  to  the  suiface  to 
produce  a  highly  enriched  layer.  This  is  not  any  sfiec  ific 
reaction  of  aluminum  but  merely  the  same  print  iple  that 
causes  soap  to  coat  the  surface  of  water.  It  is  commonly 
believed  that  (iilibs  postulated  altuKtion  of  low  energy 
roinimnents  to  the  surface,  but  this  is  out  the  case.  ’I'his 
prim  iple  ineiely  states  that  w  hen  they  are  brought  to  the 
surface  by  any  process,  ihev  will  tend  to  leinain  there 
bcfause  their  de|KuiuiT  would  cause  an  increase  in  free 
eneigy.  a  .Secoml  loiw  violation.  The  Gibl"  *'uhcm 
etpiation  inciely  pieditis  equilibrium,  and  the  rale  of 
aiuining  it  is  controlled  bv  rcl.iiisely  sluggish  solid 
diffusion. 

'I'his  raises  the  t|ucstion  of  the  lime  scale  of  Model  IIC. 
When  two  as|reritirs  ini[>ari  as  in  Fig  2,  there  ran  be  two 
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REL«Tlve  MOTION 


TIME>0 

Rg-  2— Th*  adhaiiv*  wcpr  precMi— (tort. 


TIME>I9 

Hg.  “5  — Th*  odhttiv*  w*ar  preccst— end. 


different  outcomes.  Either  the  film  survives  by  elastic 
deformation  of  the  asperities  (Fig.  4)  or  welding  takes 
place  (F’ig.  3)  with  the  generation  of  a  hemi-ellipsoidal 
wear  particle  (Fig.  (i).  (The  vertical  scale  in  all  these 
drawings  is  exaggerated  by  50X.  The  times  indicated  may 
be  regarded  as  arbitrarv'.  though  for  the  mechanism  the 
writer  had  in  mind  they  are  approximately  milliseconds.) 
Assuming  that  even  the  survival  case  rccjuires  renewal  of 
the  film  due  to  pai  iial  rupture,  we  may  calculate  certain 
times  on  ar.y  reasonable  set  of  panmeters.  Table  4 
shows  a  typical  case.  The  “wearing”  mcmbeif  were  cal¬ 
culated  on  the  basis  of  time  to  wear  a  depth  of  D,'2 
pcak  to-;.>eak  ('RMS  '/>),  uid  the  “non-wearing”  mem¬ 
bers  on  Nfodel  lA.  It  is  evideni  that  the  time  for  ditiusioii 
of  additive  from  solution  is  ver\-  limited  for  the  “fixed 
spot”  (ball)  member,  and  quite  generous  for  the  “moving 


spot”  (plane)  member,  regardless  of  where  the  wear  is 
taking  place.  Whether  the  time  is  sudicient  for  metallic  ^ 
diffusion  on  the  “moving  spot”  member  remains  to  be 
investigated.  If  it  u,  such  alloys  can  be  handled  simply^ 
by  calculating  y,  by  the  Shain  method  {30).  ‘Hm,  may  ^ 
be  a  metal,  or  sometimes  a  carbide,  etc.  y 

Table  4  represents  an  example  of  the  sort  of  new 
information  that  arises  from  this  type  of  .malysis.  N'eithei 
Model  I  nor  II  alone  would  lead  to  the  valuable  conclu¬ 
sion  that  additives  mu  t  be  tailored  to  fit  the  “moving 


spot”  member  if  they  arc  to  be  effective.  As  far  as  the 
writer  knows,  this  principle  has  never  been  recognized 
before,  but  there  is  every  reason  to  accept  it  as  valiti  based 
both  on  its  derivation  and  on  practical  expcricnee.  On 
this  basis,  JiA’  must  be  ealiulaled  for  the  “moving  spot” 
member,  but  for  th(‘  fia>di>  memlrer.  for  use  in  Equa¬ 
tion  |37].  A  mole  rigorous  solution  in  terms  of  a  anti  n, 
may  be  justifiable  in  the  fulur  '. 

Model  11  may  also  be  useful  in  predicting  the  transition 


ptrint  at  )'  in  Fig.  I,  by  reason  of  the  r.ipid  rise  of  the 


Tww*a 


double  cx|inneiilinl  in  l.(|uaiion  |I8|,  as  pointed  out  by 
Kingsbi  ry  (.//)  whose  e.iily  work  led  to  Model  ll.\.  He 
wTtu  so  far  ns  to  set  up  an  (X|uation  for  this  “rliaracteristie 
tem|xrature.“  7',: 

’■' =  4(' - 


Ng.  3— Tk*  arflwtiv*  w*«f  w*Mtng  tt««ki<ig. 


TMC<g 

Tig  4— fiKtMttwt  b**N<Ty  htk'ltafi**. 


which  drrincs  the  eontlitioiis  for  a  maximum  t/n  rf/‘. 
While  this  model  does  not  iniHiiMe  .t  vertical  litH'  tiansi- 
lion  as  shown  ie  *‘ig,  I,  the  litter  must  Ire  n'gatdrtf  .o 
schematic.  Data  in  the  }'  region  scatter  enough  to  fit  .< 
very  steep  sigmnid  curve  eqnallv  wx  II. 

In  ennrhision.  it  e.in  Ih-  s.iid  that  M<kI'-|  tl  npeio  np 
the  Vvhoir  area  of  revet  >ihle  Miifaer  enercetiis  t*>  ihi 
prrdietion  of  hihiuation  |H‘iloini.inie.  (thvioioU.  the 
genmetiA'  irf  Morkl  I  e.in  Ire  pin  into  ii  vvith  veiv  htilc 
diflUnhy.  i  Ins  vvould  poHliiee  a  rh'sign  progiani  with  ih< 
additional  inputs  show  n  in  T.ifrie  3.  I  h>vv  nu  ompletr  iho 
pirluiT  is  wilt  In  shown  ImIow.  nnihi  .MihIcI  III. 

•S 
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Table  5-ADDmONAL  INPUT  DATA  FOR 
COMPLETE  MODEL  II 

Engineering  Data;  '21  %«  ^i>  ^2<  "v  *2- 

AiMitivc  Data;  M,  F*,  S4,  y,  7",. 

BA$($  Of  MODEL  lit 

Neither  of  the  models  considered  so  far  in  this  study 
have  really  come  to  grips  with  the  irreversible  chemical 
reactions  which  are  known  to  take  place  at  wearing  sur¬ 
faces.  let  alone  those  that  are  merely  suspected.  These  may 
be  put  into  three  groups — corrosion  of  the  metals  by  the 
atmosphere,  reactions  of  the  additives  with  the  metals 
(and  their  corrosion  products)  and  polymerization  of  the 
lubricant,  with  the  metal  as  a  catalyst.  The  general  nature 
of  these  competitive  phenomena  is  illustrated  in  Figure 
6.  Truly,  lubrication  can  be  described  as  “a  very  intense 
process  in  a  very  small  reactor,”  to  quote  Dr.  David 
Tabor. 


OOHrOiCNTS 

TMCBOST-NCMMOeCM 
Ra-  Tkt  a«»S>w««r 

ANALYSIS  OF  MODB.  IIIA 

Only  a  few  papers  seem  to  have  taken  up  the  mathe¬ 
matics  of  corrosive  wear  by  oxygen  and/or  humidity.  Tao 
(32)  has  made  a  fairly  detailed  study,  based  on  the  diffu¬ 
sion  rate  of  o.xygrn  through  h)  drocarbons.  He  started 
with  ;!ie  basic  differential  equations  for  twx>-dimensional 
steady-state  diffusion,  and  some  simplifying  assumptions 
about  the  ball-on-cylinder  apparatus.  These  included 
c^nooih  surfaces,  a  uniform  clearance  ({),  fresh  oil  saiu- 
Mted  with  oxygen  and  humidity  brought  in  by  viscous 
drag,  instantaneous  reaction  of  the  ;  irfacr  followed  by 
immediate  remos'al  of  the  |-'eO(Oif),  and  wear  only  at 
tin*  upper  (ball)  surface.  Frrjm  these  an  eigenvalue  and 
eigenfunrtion  problem  was  built  up  and  solved  as 

^ _ I  r»df 

*  I.W(.‘‘b;p4  •  -  l  08cxp(-3.718A7ycn7) 

|39i 

s%herr  /  is  the  time,  C„  the  initial  oxygen  eoiKeniratkin, 
I ,  tlie  mean  .ul  flow  velocity,  p  the  oil  density,  II  the 
Hertz  cor’.iaei  diameter,  w  the  scar  width  on  the  ball 
inie-iNorctl  iKtrallel  to  tlie  C)linder  axi>),X  the  niau  diffu- 
if.  i'.  f  and  /  *  tlie  sliding  sTlocity. 


Tao  was  able  to  use  this  equation  to  compute  clear¬ 
ances,  and  obtained  quite  plausible  values  based  on  wear 
rates  in  air  by  neglecting  the  exponential  term  and  .sub¬ 
stituting  a  value  of  W  from  a  run  of  known  duration. 
This  was  based  on  the  fact  that  the  exponential  becomes 
negligibly  small  as  the  concentration  of  oxygen  in  the  oil 
leaving  the  action  zone  approaches  zero.  With  this  sim¬ 
plification, 

-  IP  =  0.082p4  CqIU  (40j 

The  value  of  Cq  was  shown  to  be  a  function  of  the 
oxygen  partial  pressure  (p),  the  density  of  the  lubricant 
at  60"F  (p„)  and  other  factors  as  defined  above,  by  the 
writer  (33). 


No  equations  were  set  up  to  handle  humidity  as  it 
appeared  that  air  of  >  50'7r  RH  would  provide  an  ade¬ 
quate  supply  of  H2O  to  keep  the  reaction  supplied.  Tao 
demonstrated  in  fact,  that  F.qu.^tion  [40)  was  only  in 
error  by  about  10*^  regardless  of  the  humidity. 

A  tentative  dependence  of  w  on  the  ‘‘compliance",  a 
complex  term  involving  If  and  D.  was  demonstrated  but 
the  mrlaiioiwliiii »«  Mode!  !B  cannot  b?  ect.nbli'.hcd  v.’ith- 
out  a  good  deal  more  work. 

Obvioiitily,  this  is  only  a  partial  model,  as  Equation 
(40)  would  indicate  that  no  wear  takes  place  at  C„  =  0 
This  is  not  true,  as  shown  in  subsequent  work  by  the  same 
group  (?)  with  argon.  However,  no  attempt  wa.s  made 
then  to  incorporate  this  observation  into  .Model  III.A. 

A  very  recent  paper  by  .Schaizberg  (3-f)  led  to  some 
comments  by  Tao  and  the  wiiter  which  might  be  con¬ 
sidered  to  cunstituie  an  extension  to  inert  atmospheres. 
The  observatioft  vvas  that  the  dry  argon  had  a  low  but 
measurable  wear  rate,  consistent  with  (?)  Schaw'berg 
reduced  .ill  of  his  wear  rates  to  the  forn. 

log  fF  =  m  log  (  -f  C  (42) 

where  C'  is  a  constant  for  anv  given  system.  He  was 
somewhat  surprised  to  find  that  only  four  of  his  sixteen 
cases  fit  tlir  )Kiiirrn  of  Fa)uaii«>n  (16)  wht'  h  Feng  I  ?.i) 
had  shown  to  be  equivalent  to  E()u.xtk>n  |  I2j  sviih 
«  0.2j.  Kight  of  tlie  cav'  shovved  m  =:  0  (t8  to  i>  If) 

four  showed  w  =  0.16  to  O.IH.  and  lour  showed  m  =  O.J  j 
to  0.29  Tltr^’*  results  seem  to  be  strongly  rel.ited  to  those 
of  Bayer  (7).  since  7.  is  pmtMuitonal  to  I.  .Mlnwing  fi»r 
tlie  diireirnec  in  gronutry  due  to  Sshai/lK-fe‘\  use  of  the 
4-batl  ina«  hine  ifor  whk  li  Lqii.riion  |  S<))  h.o  not  vet  Ihtii 
integrated;,  it  vretiis  note-  otihy  that  his  d.tt.i  fall  into 
twx»  modes  ne.»t  llayet"'  pirriirtions  of  m  0  IIR  for 
low-transfer  wear  and  m  »  <)  Jjt)  for  high  ir.msfri  wear 
of  the  ball  Haver  alsvt  irjauied  t*^l  that  a  h  's  eonthina- 
tkmt  srented  to  have  .hi  unstable  mode  wiih  m  taking 
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interniediatc  values.  Schalzberg’s  low  values  were  associ¬ 
ated  with  dry  oxygen  and  argon,  and  his  high  ones  with 
wet  oxygen.  Thus,  the  "Archard  Wear  Law"  of  Equation 
(16]  may  need  to  be  generalized  to  the  form 

„  A-yn  yVd*”* 

V  =  -  [431 


where  m  will  vary  with  the  atmosphere  in  two  or  three 
steps. 

So  far,  all  the  tests  of  Model  lll.\  have  involved  52100 
steel.  The  possibility  of  a  mathem.itieal  basis  for  extend¬ 
ing  this  to  other  metals  and  alloys  is  implied  by  Fuwlc 
(36).  He  studied  the  ratios  of  oxides  formed  on  alloys  and 
found  that  ih"  fraction  coated  with  each  oxide  wz  more 
closely  related  to  the  free  energy  of  formation  of  the  oxide 
than  to  the  comjiosition  of  the  alloy.  This  may  not  be 
helpful,  as  rates  of  oxide  formation  are  not  readily  pre¬ 
dictable  fiitm  these  free  energies  according  to  Kubaschew- 
ski  (37).  .Mso,  in  .ill  of  I'owle's  live  examples,  the  (Jibbs- 
Duhem  principle  ritcil  at  the  etui  of  Model  llCl  would 
have  predictrd  similar  oxide  ratios.  However,  the  possi¬ 
bility  of  extending  Kt|uatiun  (43|  to  other  metals  is  very 
real. 

While  Tao  assumed  for  rntnlel  jnirposes  that  the  corro¬ 
sion  product  was  remosed  as  fast  as  formed,  this  is  not 
rtall)  the  case.  .-Vs  a  rc.aili,  duic  will  be  a  need  lor  surf.ice 
free  energies  of  such  comitounds  in  the  solid  state.  Duga 
(22)  has  made  an  exhaustive  study  of  thc.se  inaieiials  and 
his  irsults  can  be  applied  to  a  hvbrid  Model  IKi  lll.-V 
Uewig  and  /isman  (  W)  have  studied  metal  surfaces  poi- 
ished  under  w.rter  with  what  may  l>e  presumed  to  be  an 
oxide-hydroxide- water  film  left  <»n  them.  The  d.ata  are 
quantttalixely  suitable  for  insertion  into  Equation  (24). 
Stud)  (.i^)  has  studied  the  displacement  of  one  additive 
by  another  this  can  l>e  explainet'  by  Mwlel  IKl,  but 
i  may  involve  Ill.-V  concepts  tucomplite  it  as  .-\nderson  (^t)) 
y shows  similar  dispta  rinent  by  huinidits. 

he  comTpt  of  additives  ioiming  cov.^lent  i-omixmnds 
jwith  the  met.’il  surface  is  at  least  as  old  .ts  that  of  reversi¬ 
ble  adsoqition  ( />.  However,  no  m.iiheinatical  mtKiel  has 
Ux-n  .suggested  up  to  now.  The  ne.uesi  appniach  has  lieen 
that  of  Kreui!  (T/l  who  tletuutistraterl  that  sue!,  hints  can 
Ire  lenwrved  by  the  solseni  action  of  the  base  llnitl,  I  Ins 
bad  alwayv  been  a  stimig  aigttment  lor  tesTrsibIr  adsoi  |>- 
tioh.  as  it  was  hard  to  explain  hosv  a  covalent  lihn  rould 
ever  fail.  llowcsTr.  Kreuz  did  not  varry  his  ssori  to  its 
hrgiral  eomlusion  hs  applying  the  Hildebrand  th<x»rv  of 
regular  M>hitioi«s(  /2i  This  m.iv  have  l»een  h-Tausr  Krruz 
found  quite  sharp  "vihibih/aiion"  trntfietaturrs  for  his 
s«»ap  fihns,  which  svert  esseiitiallv  im!c|ri-iidmt  of  txrn* 
(Totration;  m  fact,  ihrv  rtnild  Ire  relates!  Isr  tire  liansiihm 
})oini  at  I"  in  Eig  1.  Sueh  lielMsnw  ts  not  pnxlietril  by 
the  Hildrbiaiul  th<s  is  isesaiio’  it  tsi.ues  only  irr  t.»iitlon 
forsi-s  ami  the  wraps  ntelude  in.iior  jwdar  and  hvthogen 
Umding  I’redultoo  of  the  Iwh.isior  o(  sovh  'S*nms  n*. 
qsiiirs  the  llan'Cn  nHHhhi.ilion  t/'f).  whtsh  pteiints  a 
i.uh<»  sharp  inivihilns  when  the  molar  vohime  of  ihe 
vs'hil  i\  vsxy  large  ssith  rrsjittt  ts»  lh.it  of  tite  b.tv'  flsiul, 
I  • 


(5gi  -  ««)'  +  0-25(fipi  -  «p2)* 

RT 

+  0.25(6„,-«^)*<^  (441 

where  6,, ,  is  the  hydrogen  bonding  parameter  of  the  soap, 
Sfi2  that  of  the  base  fluid  and  the  other  variables  are  as 
previously  defined  under  Model  IL'\  This  method  will 
require  more  work,  especially  since  6^  and  6^  for  soaps 
are  not  so  easily  calculated  as  for  polymers,  but  it  docs 
oifer  design  possibilities.  It  also  has  the  complication  that 
it  leads  back  to  Model  IIC,  since  the  soaps  are  still 
"ndditives"  and  must  be  treated  according  to  their  £ 
values  on  the  metals. 

ANALYSIS  OF  MODEL  IIIC 

•'(s^  friction  polymer”  concept  is  a  relatively  new  one, 
.111-1  w  still  far  from  computable.  It  was  so  named  in  1958, 
but  I'cin  (43)  was  the  first  to  put  it  into  general  use.  The 
genera!  nature  of  the  reaction  seems  to  be  clear,  it  consists 
of  the  highly  energetic  surface  of  the  metal,  made  even 
more  active  by  rc-cent  abrasion,  acting  as  a  catalyst  in 
the  decomposition  of  the  base  fluid.  The  fragments  then 
recombine,  with  atmospheric  oxygen  if  available  (more 
slowly  without  it)  to  form  a  low  molecular  weight  poly¬ 
mer.  or  oligomer  by  the  usual  standards  cf  polymer 
chcmistiy-.  Chcmic.il  analyses  of  the  product  arc  rare  and 
iiK'umplete,  but  Fein  found  that  it  showed  ketone  and 
ester  groups,  and  even  (when  nitrogen  was  available) 
amides  in  the  infrared  s|xvtra.  This  phenomenon  has  the 
pos.sihility  of  explaining  some  of  the  anomalies  not  han- 
dietl  by  any  other  model  (44)  and  also  may  shed  light 
on  some  morlcs  of  U-aiing  failure,  outside  the  scope  of 
this  study,  due  to  “lacquer"  etc. 

Fein  had  little  to  olTrr  in  the  litw  of  a  iiiaibcmaticai 
model.  Bond  (45.  46).  who  has  done  much  pioneer  vvoik 
in  clarifying  the  rank  of  metals  as  catalysts,  has  not  been 
able  to  provide  quaniitaiivT  prrdiclinn  methods.  This  is 
partly  bvx  ausr  the  iiieials  do  nut  rank  in  quite  the  same 
order  for  different  renciiuns.  For  c.xainple,  in  hydrogen 
ciienuiorptinn  energy  they  rank  as  fu!!'i.'’s‘ 

IF  >  Ta  >  Mo  >  Fe  Ni  Pd  >  Rh 

while  hw  chemisorption  of  eihylrtte  they  rank: 

1*1  Zi  It  >  Pd  >  (Fe,  Ca*.  Ni)  2s  Rh  >  Ro  2:  0»  >  Cu 

Both  ortlen  ftillow  rather  rliMcly  that  of  the  uirface  free 
energies  (2^.  (H>vi.nidy.  the  nature  of  the  base  fluid 
makes  a  great  dilfrtriHr.  aiul  Tabor  has  found  (47)  that 
even  « upper  will  r,ipidlv  ptilynieiire  dittK-ihylMlitrrne. 
Bond  (^^5  ptrrlMied  ih.st  the  mder  of  activity  of  hydro- 
carlMin  would  iriKl  to  be; 

ohfuiv  >  iv<|Mralhiis  >  naphthenes 

>  aouualics  >  |Mtaflim 

'Fhe  work  of  the  latalysi  chemiits  is  based  on  aged 
swtfasrs  and  may  not  irv»-al  ilie  full  extent  of  the  |>rob- 

•T 
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Ta«.e  6- PREDICTIONS  OF  SCORING  RESISTANCE  ON  IRC  JN 


MisciBiir-.  i'r.MPi-PAn m 


At. 

No. 

Rating  feom  LiTEEAnaE 

Equation 

(V  -  r)(48 

Metal  Cevstal 

(54) 

(54*) 

(55) 

(55*) 

deg  K 

Sua;tis 

4 

Be 

HCP 

Very  Poor 

— 

— 

— 

rdO 

+  • 

12 

Mg 

HCP 

Poor 

— 

Good* 

Good/Fair 

108''3 

-*-c 

13 

Al 

FCC 

Poor 

Poor 

— 

£  atr 

2010 

14 

Si 

DIA 

Vfery  Poor 

— 

— 

— 

1976 

- 

20 

Ca 

FCC 

Very  Poor 

— 

— 

— 

24531 

- 

22 

Ti 

HCP 

Very  Poor 

— 

— 

Fair 

1141 

„  t 

24 

Cr 

BCC 

Very  Poor 

— 

Good 

- 

146 

_•  * 

26 

Fe 

BCC 

Very  Poor 

— 

— 

— 

0 

+ 

27 

Co 

HCP 

Very  Poor 

Poor 

— 

— 

139 

+ 

28 

Ni 

FCC 

V'ery  Poor 

— 

Poor 

— 

84 

+ 

29 

Cu 

FCC 

Fair/Poor 

Poor 

— 

Poor 

179 

+ 

30 

Zn 

HCP 

Poor 

Poor 

Fair 

Fair 

7004 

-*-r 

32 

Ge 

DIA 

Good 

— 

— 

— 

4370 

+ 

40 

Zr 

HCP 

Very  Poor 

— 

Poor 

Fair 

1281 

41 

Cb(Nb) 

BCC 

Very  Poor 

— 

— 

— 

216 

+ 

42 

Mo 

BCC 

Very  Poor 

— 

— 

— 

246 

+ 

45 

Rh 

FCC 

Very  Poor 

— 

— 

— 

275 

+ 

46 

Pd 

rCC 

Very  Poor 

— 

— 

— 

450 

+ 

47 

Ag 

FCC 

Good 

— 

Fair 

Gcxxi 

25% 

48 

Cd 

HCP 

Good/Fair 

Good 

Good* 

Fair 

12170 

t 

49 

In 

FCC 

Good 

— 

— 

— 

8246 

+ 

50 

Sn 

DIA 

Good 

_ 

— 

— 

6%7 

+ 

51 

Sb 

LAM 

Good 

— 

— 

— 

10639 

+ 

56 

Ba 

BCC 

Pout 

— 

— 

— 

388% 

- 

58 

Cc 

HCP 

Very  Poof 

— 

— 

— 

IStiilU 

73 

Ta 

BCC 

Ver>  P'.>or 

— 

Fair 

— 

808 

74 

W 

BCC 

Poor 

— 

Fair 

1609 

*1 

77 

Ir 

FCC 

Very  Poor 

— 

— 

- 

950 

+ 

78 

?i 

FCC 

Wry  Poor 

— 

— 

- 

32 

+ 

79 

Au 

FCC 

\try  Poor 

- 

— 

1221 

»• 

8! 

Tl 

HCP 

Good 

— 

— 

— 

14101 

+  ' 

82 

Pb 

FCC 

Good 

— 

Good* 

Good  Fait 

11742 

^  * 

83 

Bi 

LAM 

Good 

— 

Fair 

16901 

»  t 

90 

Th 

FCC 

Very  Poor 

— 

- 

4226 

- 

92 

U 

BOC 

Ner>  Poor 

-■ 

715 

+ 
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Taele  7-PREDlCTION  OF  SCORING  RESISTANCE 


Metal  Pair 

Crystal 

Rating  from  Literature 

MiaciBiLFTY  Temperature 
Equation  (F  —  I)  (48) 

(5'*) 

(55) 

(55*) 

degK 

Success 

Al-Zn 

FCC-HCP 

Poor 

Fair 

1879 

±c 

Co-Cu 

HCP-FCC 

Poor 

— 

622 

AJ-Co 

FCC-HCP 

Poor 

— 

3214 

Cu-Cd 

FCC-HCP 

Poor 

— 

Poor 

9706 

Cu-Zn 

FCC-HCP 

Poor 

— 

Poor/Fair 

4867 

— 

Sb-Zn 

LAM-HCP 

Poor 

_ 

_ 

3 

Al-Cd 

FCC-HCP 

Good 

_ 

_ 

4771 

Bi-Cd 

LAM-HCP 

Good 

_ 

_ 

36 

Cd-Zn 

HCP-HCP 

Good 

_ 

459 

Bi-Zn 

LAM-HCP 

Good 

_ 

_ 

383 

Cu-Ni 

FCC-FCC 

— 

Poor 

498 

+ 

Cr-Mo 

BCC-BCC 

— 

Poor 

— 

825 

+ 

Ti-Zn 

HCP-HCP 

— 

Poor 

— 

0 

Ag-Cu 

FCC-FCC 

— 

Poor 

Fair 

1334 

Cr-Ni 

BCC-FCC 

— 

Poor 

— 

444 

+ 

A1-\V 

FCC-BCC 

- 

Poor 

— 

8526 

- 

Cu-Sn 

FCC-DIA 

— 

Good** 

4583 

+ 

Ni-Sn 

FCC-DIA 

— 

Fair 

— 

8580 

- 

Al-Zn 

FCC-HCP 

Poor 

Fair 

188 

Ag-Ta 

FCC-BCC 

— 

Good 

_ 

7720 

+ 

Cu-Mo 

FCC-BCC 

— 

Very  Good 

— 

903 

Ag-Cr 

FCC-BCC 

— 

Fair/Good 

— 

1487 

- 

Cu-Ta 

FCC-BCC 

— 

Good 

— 

1898 

— 

Al-Pb 

FCC-FCC 

r.  ?- 

A  ait 

4043 

1  tk 

Cu-W 

FCC-BCC 

— 

Good 

2986 

— 

Cu-Pb 

FCC-FCC 

_ 

_ 

Poor 

8527 

— 

Cu-Ti 

FCC-HCP 

_ 

_ 

Poor/Fair 

367 

AlCu 

FCC-FCC 

_ 

Poor 

948 

■f 

Cu-Mg 

FCC-HCP 

Poor 

7933 

-.r 

Pb-Ti 

FCC-HC:P 

_ 

_ _ 

Fair 

6353 

—  e 

Zn-Zr 

HCP  HCP 

_ 

Poor 

3636 

_c 

Ti-Zn 

HCP-HC;P 

_ 

Fair 

3210 

Ag-Zn 

FCC-MC:P 

_ 

Fair 

1401 

Ag-Al 

FCC-FCC 

_ 

Poor 

41 

+ 

Al-Mg 

FCC-HCP 

*•— 

Poor 

3806 

_c 

Al-Ti 

FCC-HCP 

Fair/Poor 

166 

Al-Ni 

FCC-ltIC 

— 

Fair 

2901 

+ 

Ag-Cu 

FC.C-1'C.C 

_ 

Poor 

1335 

+ 

Ag-Ti 

Fcc-H(;;p 

_ 

_ 

Fair 

378 

—  c 

Ag-Zr 

FCC-HCP 

Good 

430 

—  c 

tVI  Mg 

HCP  ik:p 

_ 

1  oor 

84 

+*■ 

Cd-'I'i 

HCPHC;P 

_ _ 

Poor 

7062 

_c 

Ol-Zr 

HCP-HCP 

Good 

8095 

Mg-Ti 

HCP  HI:p 

_ 

Poor 

5897 

Ag-Mg 

.vc:c  ;p.:p 

— 

Fair 

3056 
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lem.  Wear  spots,  such  as  shown  in  Fig.  6,  are  freshly 
abraded  and  emit  electrons  almost  as  if  they  were  radio* 
active.  Polymerization  by  electrical  discharge  has  been 
studied  by  Bradley  (■#P),  who  ranked  some  hydrocarbons; 

styrene  >  naphthalene  >  toluene  >  ethylene 

>  hexamethylbenzene  >  propane 

which  again  point  up  the  uncertainties. 

One  factor  not  previously  considered  is  the  effect  of 
pressure  on  polymerization.  Wcale  (50)  gives  an  equation 
which  is  adaptable  as  follows: 

ln(W  =  |^  (45] 

where  k  and  arc  the  rate  constants  at  pressure  and 
at  one  atmosphere  respectively,  and  is  the  frac¬ 

tional  change  in  volume  on  polymerization.  This  predicts 
fora  typical  reaction  inv\’  ich  =  0.20,  an  accel¬ 

eration  of  approximately  10,000  fold  under  =  20,000 
atmospheres. 

Regardless  of  the  lack  of  consistent  data  and  of  a 
complete  mathematical  model,  this  line  of  attack  looks 
very  promising.  The  term  “surface  resin”  for  the  material 
is  suggested. 

A  very  recent  paper  by  Thompson  (57)  may  be  useful 
as  it  takes  into  account  the  fact  that  the  anti-asperities 
can  carry  part  of  the  load  by  hydrodynamic  lubrication. 
Since  Model  IIIC  involves  partial  filling  of  these  spaces, 
Thompson’s  equation  may  provide  a  means  for  testing 
the  rate  of  surface  resin  formation  by  following  the  time 
required  to  break  in  a  pair  of  surfaces.  This  concept  can 
also  be  used  to  predict  another  “transition  failure”  mech¬ 
anism  via  Equation  [41].  This  would  result  from  the 
layers  of  “surface  resin”  formed  in  the  anti-asperities 
during  break-in,  (Fig.  6)  suddenly  going  into  solution  {52) 
thus  restoring  the  original  surface  roughness.  The  result 
would  be  to  shift  a  good  deal  of  the  load  to  the  asperities, 
which  would  fail  quite  rapidly. 

MODEL  IV 

Investigation  of  the  transition  temperature  has  already 
been  mentioned  under  Models  I!  and  III.  However,  theic 
is  reason  to  believe  that  when  all  else  (including  the  oxide 
film)  fails,  the  nature  of  the  metal-metal  contact  may 
make  the  dilfcrence  between  corrosise  wear  and  scoring 
at  J'or  even  welding  at  Z (Fig.  1 ).  Ernst  (53),  Roach  (5-f), 
and  Coffin  (55)  have  all  made  experiments  along  this  line. 
About  80'^f  of  their  results  can  be  predicted  by  a  simple 
equation  from  Hildebrand  (42)  for  the  miscibility  tem¬ 
perature  of  metals; 

=  |46| 

where  K,  and  Fj  arc  the  molar  volumes  of  the  metals 
while  5,  and  arc  their  .solubility  parameters.  If  T,  is 


less  than  7',,  scoring  will  take  place.  The  results  of  these 
predictions  are  shown  in  Tables  6  and  7.  The  success 
ratings  had  to  allow  for  'he  results  from  the  laboratories 
containing  contradictioi.i,  par  of  which  may  have  been 
caused  by  the  lack  of  a  clear  definition  of  “scoring”.  The 
writer  attempted  to  bring  the  data  to  a  common  basis, 
in  which  “rough”,  “abraded”,  etc.,  were  converted  to 
“poor”  scoring  resistance,  “smooth”  or  low  wear  rates  to 
“good”,  and  intermediate  descriptions  to  “fair”.  Where 
data  in  helium  was  available,  it  was  given  preference  over 
that  in  air  or  in  air-saturated  kerosene.  All  hexagonal 
(HCP)  metals  were  excluded  from  the  ratings;  as  shown 
by  Buckley  (55),  these  follow  different  rules  than  the 
common  cubic  metals.  The  criteria  of  success  were  based 
on  an  arbitrary  choice  of  Poor  =  0-2000  °K,  Fair 
=  2000-4000  °K,  Good  =  4000 -f  °K.  These  unrealistic 
temperatures  result  from  their  being  based  on  complete 
miscibility  at  about  50/50%  by  volume.  Future  programs 
should  be  keyed  to  the  fully  developed  version  of  Equa¬ 
tion  (46),  where  this  can  be  set  at  5/95  or  10/90.  As 
pointed  out  by  Hildebrand  (42),  peH'ect  success  is  unlikely 
due  to  intermetallic  compound  formation  and  also  to 
crystalline  transitions  of  some  metals  at  elevated  temper¬ 
atures. 
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-'^Mechanical  and  lubrication  engineers  picture  solid-solid  contact  as  taking 
place  only  on  the  tips  of  the  higher  asperities  of  each  surface,  so  that  the  contact 
load  is  borne  only  on  a  very  small  portion  of  the  total  surface  area  and  all 
frictional  heat  is  generated  and  absorbed  in  this  small  area.  The  size  and  number 
of  the  Individual  asperity  Junctions  is  in  general  not  known;  but  averages  are  known, 
and  it  appears  that  asperity  heights  and  spatial  frequencies  are  described 
reasonably  well  by  Gaussian  statistics.  This  should  make  it  possible  to  calculate 
frictional  hot-spot  temperatures  and  distributions  on  an  explosive  surface,  giv&n 
some  laboratory  measurements  of  factors  such  as  coefficients  of  friction  frictional 
work  expended  and  surface  characteristics.  With  hot-spot  statistics  available,  the 


probability  of  explosion  can  be  calculated  by  available,  published  methods. 

*  A  coordinated  theoretical  and  empirical  research  program  is  proposed  to 
apply  the  above  engineering  and  mathematical  concepts  to  explosives.  The  results 
should  facilitate  the  identification  of  basic  mechanical  processes  (such  as  metal- 
metal  friction,  or  the  viscous  or  plastic  shear  of  explosives)  which  do  and  do  not 
generate  explosion  hazards.  This  in  turn  will  aid  the  design  of  more  valid 
laboratory  friction  sensitivity  testers  and  the  design  of  safer  explosive  processing 
plants.^^^ 
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